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ABSTRACT 


In  this  Annual  Report  of  Electronics  Research  at  The  University 
of  Texas  at  Austin,  we  describe  recent  findings,  progress,  and  future 
plans  in  the  areas  of  solid  state  electronics,  information  systems,  elec- 
tronic computers,  electronic  controls,  plasma  and  quantum  electronics, 
and  radio  sciences.  Among  the  many  project  reports  in  these  areas  are 
a number  of  significant  findings.  In  the  solid  state  area,  a study  of  the 
transition  layers  between  CVD  VO2  films  and  oxide  substrates  reveals 
complex  substrate- film  interactions  with  significant  implications  for  VO^ 
and  other  material- oxide  substrate  systems. 

The  investigation  of  computer  architecture  for  byte  string  and 
vector  operations,  in  the  information  systems  area,  shows  iterative  pro- 
cessing of  variable  length  operands  to  be  feasible.  In  the  area  of  elec- 
tronic computers,  errors  have  been  reduced  and  reliability  increased  with 
minimal  redundancy,  through  the  use  of  functional  simulation  of  logic  nets. 

A fourth  significant  report  appears  in  the  electronic  controls 
section.  An  easily  implementable , suboptimal  estimator  has  been  devel- 
oped and  applied  to  problems  in  satellite  orbital  perturbations  and  in 
optical  phase  tracking. 

There  were  two  developments  of  interest  in  the  plasma  and  quan- 
tum electronics  area.  First,  digital  bispectral  analysis  has  proven  to  be 
a useful  tool  for  the  interpretation  of  nonlinear  wave  interaction  data. 

In  particular,  this  technique  permits  one  to  discriminate  between  spon- 
taneously excited  and  coupled  modes,  and  to  measure  the  strength  of  the 
interaction  or  coupling  coefficient.  A second  significant  article  reports  a 
general  unifying  theory  for  incoherent  optical  convolvers  and  correlators. 
Improved  flexibility  is  gained  through  relaxed  geometric  constraints.  A 
final  report  in  the  radio  sciences  area  summarizes  improvements  in  tech- 
niques for  on-site  magnetotelluric  data  analysis  and  interpretation. 
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Mirror  Confined  Plasma" 

* E.J.  Powers  and  Y.  C . Kim,  " Bispectral  Wave 
Analysis  of  Nonlinear  Wave-Wave  Interactions 
in  Plasmas" 

University  of  Sterling 
Sterling,  U.K. 

May  2 5 , 1976 

* M.  Fink,  "Relativistic  and  Correlation  Effects 
in  Atoms  Studied  by  Electron  Scattering" 
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University  of  Saint  Andrews 
Saint  Andrews , U,K, 

May  28,  1976 

* M.  Fink,  "Charge  Density  Distributions 
of  Molecules" 

University  of  Manchester 
Institute  of  Science  and  Technology 
Manchester,  U.K. 

June  1,  1976 


* M.  Fink,  "Structure  Determination  by 
Electron  Scattering  and  Counting  Machines" 

University  of  Reading 
Reading,  U.K. 

June  3 , 1976 


* M.  Fink,  "Charge  Density  Distributions 
of  Molecules" 

IX  International  Quantum 
Electronics  Conference 
Amsterdam,  The  Netherlands 
June  14-18 , 1976 

* M.  F. Becker,  “Third  Harmonic  Generation 
at  10 . 6 Microns" 

International  Federation  of 
Automatic  Control  Symposium 
on  Large  Scale  Systems  Theory 
and  Applications 
Udine,  Italy 
June  16-20  , 197 6 

* B.J.  Olufeagba,  R.H. Flake,  and  K.J.  Almquist, 

"A  Multiple  Shooting  and  Sweep  Algorithm  for 
Optimal  Point  Controlled  Distributed  Parameter  Systems" 
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University  of  Bielefeld 
Bielefeld,  W.  Germany 
June  17,  1976 

* M.  Fink,  "Relativistic  and  Correlation  Effects 
in  Atoms  Studied  by  Electron  Scattering" 

University  of  Muenster 
Muenster,  W.  Germany 
June  23,  1976 

* M.  Fink,  "Charge  Density  Distributions  of  Molecules" 

University  of  Kaiserslautern 
Kaiserslautern,  W.  Germany 
June  30,  1976 

* M.  Fink,  "Relativistic  and  Correlation  Effects 
in  Atoms  Studied  by  Electron  Scattering" 

NASA  Seminar/Research  Institute  on 
Differential  and  Algebraic  Geometry 
for  Control  Engineers 
Moffett  Field,  California 
June  30,  1976 

* Steven  I.  Marcus,  "Modeling  and  Stability 

of  Linear  Systems  with  Multiplicative  Poisson 
Noise" 


Gordon  Research  Conference  on 
Lasers  in  Medicine  and  Biology 
Meriden,  N,H. 

July  5-9,  1976 

A,J.  Welch  and  Les  Priebe,  "Prediction  of 
Retinal  Thermal  Injury  from  Laser  Radiation" 

A.J.  Welch  and  Mike  Mandell,  "Receptive 
Field  Measurements  After  Laser  Radiation" 
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University  of  Tuebingen 
Tuebingen,  W.  Germany 
July  20,  1976 

*M.  Fink,  "Structure  Determination  by 
Electron  Scattering  and  Counting  Machines" 

13th  International  Conference  on 
the  Physics  of  Semiconductors, 

Rome,  Italy 
August,  1976 


*R,  M.  Walser  and  R,W.  Bene,  "Solid  Phase 
Reactions  at  Transition  Metal- Silicon  Interfaces" 

Symposium  on  the  Nature  and 
Physical  Properties  of  the 
Earth's  Crust 
Vail,  Colorado 
August  2-6 , 1976 

*F.X.  Bostick  and  H.W.  Smith,  " Magnetotelluric 
and  Dipole-Dipole  Soundings  in  Northern  Wisconsin" 

13th  IEEE  Computer  Society 
International  Conference 
Mayflovrer  Hotel 
Washington,  D.C. 

September  1976 

T.K.M.  Agerwala,  "Designing  with  Microprocessors" 

Division  Meeting  of  Optical  Society 
of  America 
Houston,  Texas 
Septembers,  1976 

*J.W.  Keto,  "Kinetic  Studies  of  Excimer  Lasers 
Using  Single  Photon  Timing" 
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Rice  University 
Houston,  Texas 
September  28  , 1976 

★ Steven  I.  Marcus,  "The  Use  of  Harmonic 
Analysis  in  Suboptimal  Estimator  Design" 

The  1976  International 
Telemetering  Conference 
Los  Angeles,  California 
September  28-30 , 1976 

G.L.  Wise  and  J.B.  Thomas,  "Zero  Memory 
Nonlinear  Transformations  of  Gaussian  Processes" 

Fourteenth  Annual  Allerton  Conference 
on  Circuit  and  System  Theory 
Monticello,  Illinois 
September  29  - October  1,  197  6 

N.C.  Gallagher,  J.W,  Allen,  and  G.L.  Wise, 
"Fourier  Series  Representation  for  Polynomials 
with  Application  to  Nonlinear  Digital  Filtering" 

G.  L.  Wise  and  N.C.  Gallagher,  "A  Represen- 
tation for  Spherically  Invariant  Random  Processes" 

International  Conference  on 
Software  Engineering 
San  Francisco,  California 
October  1976 


E.  W.  Thompson,  "A  Data  Structure  and  Drive 
Mechanism  for  a Table-Driven  Simulation  System 
Employing  Multilevel  Structured  Representations 
of  Digital  Systems" 
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Euromicro  Symposium 
Venice 

October  12,  1976 

G.  J.  Lipovski,  "The  Architecture  of  a Simple 
Effective  Control  Processor" 

EPA  University  Noise  Seminar 
Purdue  University 
October  18-20 , 1976 

E.  L.  Hixson,  "Vehicle  Noise  Source  Level 
Measuring  System" 

ACM  'll 
Houston,  Texas 
October  20 , 197 6 

G.  J.  Lipovski,  "Panel  Session:  Micros- 
Minis-Maxis" 

U.S.  Geological  Survey  Workshop  on 

Electrical  Methods  in  Geothermal  Exploration 
Snowbird,  Utah 
November  4-7,  1976 

*H.W.  Smith,  "Real  Time,  On-Site  Processing 
of  Magnetotelluric  Data" 

F. X.  Bostick,  "A  Simple  Almost  Exact  Method 
of  MT  Analysis" 

F.X.  Bostick,  "The  Effect  of  Lateral  Inhomogeneities 
in  a Magnetotelluric  Survey" 
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Chemical  Engineering  Department 
The  University  of  Texas  at  Austin 
Austin,  Texas 
November  9 , 1976 

* Steven  I.  Marcus,  "The  Use  of  Harmonic 
Analysis  in  Suboptimal  Estimator  Design" 

IEEE  International  Symposium  on 
Pattern  Recognition 
San  Diego,  California 
November  9-11,  1976 

* T.  J.  Wagner,  "A  Distribution-Free 
Performance  Bound  in  Error  Estimation" 

The  American  Physical  Society  Division  of 
Plasma  Physics  Meeting 
San  Francisco,  California 
November  15-19 , 1976 

* Y.  C,  Kim  and  E,  J.  Pov;ers,  "Utilization 

of  Bispectral  Analysis  in  Experimental  Studies 
of  Nonlinear  Wave-Wave  Interactions" 

E.r.  Powers,  et  ai.  , "Fluctuation  Spectrum 
in  the  NASA  Lewis  Bumpy  Torus" 

Surface  Science  Seminar 
Bell  Telephone  Labs 
Murray  Hill,  New  Jersey 
December  1976 

*R.W,  Bene  andR.M.  Walser,  "Schottky  Barriers 
at  Transition-Metal-Silicon  Interfaces" 
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Surface  Physics  Seminar 
Bell  Telephone  Labs 
Murray  Hill,  New  Jersey 
December  1976 

* R.  M.  Walser  and  R.W.  Bene,  "Nucleation  of 
Silicides  at  Transition  Metal-Silicon  Interfaces" 

IEEE  Conference  on  Decision  and  Control 
Clearwater  Beach , Florida 
December  1-3,  1976 

* Steven  I.  Marcus,  'The  Use  of  Spherical 
Harmonics  in  Suboptimal  Estimator  Design" 

Eall  Annual  Meeting 
American  Geophysical  Union 
San  Erancisco,  California 
December  6-10,  1976 

F.X.  Bostick,  H.W.  Smith  and  J,E.  Boehl, 

" Magnetotelluric  Measurement  over  Kilauea 
Iki  Lava  Lake" 

Second  International  Conference  and 
Winter  School  on  Submillimeter  Waves 
and  Their  Applications 
San  Juan,  Puerto  Rico 
December  6~11,  1976 

A.  B.  Buckman,  "Theory  of  an  Efficient  Electronic 
Phase  Shifter  Employing  a Multilayer  Dielectric- 
Waveguide  Structure" 

1977  IEEE  International  Symposium 
on  Circuits  and  Systems 
Phoenix,  Arizona,  1977 

H.  Chang  and  J.  K.  Ag garwal , "Design  and  Simulation 
of  Two-Dimensional  Interpolated  Filter  Systems" 
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10th  International  Hawaii  Conference 
on  System  Sciences 
January  6-7 , 1977 

* T.K.M.  Agerwala  and  B.  Lint,  "A  Comparison  of 
Synchronizing  Primitives" 

University  of  Alberta 
January  23  , 1977 

G,  J.  Lipovski,  "On  Virtual  Memories  and 
Microprocessors" 

4th  Annual  Conference  on 
The  Physics  of  Compound 
Semiconductor  Interfaces 
Princeton,  New  Jersey 
February  1977 

* R.  V\7.  Bene  and  R.M.  Walser,  "The  Effect 
of  a Glassy  Membrane  on  the  Schottky  Barrier 
Between  Silicon  and  Metallic  Silicides" 

Solid  Phase  Epitaxy  Program  Review 
California  Institute  of  Technology 
Pasadena,  California 
March  1977 

R.  M,  Walser,  "Recrystallization  of  Ion-Implanted 
Amorphous  Silicon  Layers" 

R.  W.  Bene,  "Silicide-Silicon  Barriers" 

Intel  Corporation 
March  1977 

R,  M.  Walser,  "Electronic  Conduction  in  Ion- 
Implanted  Poly  Silicon" 
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INTER-NOISE  77 
Zurich,  Switzerland 
March  1977 


E,  L.  Hixson,  "The  Qualifications  for  a 
Noise  Control  Engineer" 

14th  Annual  Electronics  Research  Review 
The  University  of  Texas  at  Austin 
Austin,  Texas 
March  8-9 , 1977 

* H.  W.  Smith,  Overview  of  "Modeling  Subsurface 
Resistivities  with  Magnetotelluric  Data  (MT  Data)" 

B.  F.  Womack,  Overview  of  "Electronic  Control 
Systems" 

* S.  I.  Marcus,  "The  Use  of  Harmonic  Analysis 
in  Sub-Optimal  Estimator  Design" 

G,  L.  Wise,  "An  Overview  of  Some  Recent  Results 
in  Communication  and  Control" 

* J.  K.  Aggarwal,  "Recent  Results  on  Two-Dimensional 
Signal  Processing" 

* T.  J.  Wagner,  "Recent  Results  on  Nonparametric 
Estimation  and  Discrimination" 

* S.A.  Szygenda,  Overview  of  "Electronic  Computer 
Software  System"  and  "Electronic  Computer  Design 
and  Analysis" 

G.  I.  Lipovski,  " Computer  Architecture  and  Micro- 
processor Systems" 

* T.  K.  M.  Agerwala,  "Studies  in  Concurrent  Computer 
Systems" 
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14th  Annual  Electronics  Research  Review 
(continued! 

* E.  W.  Thompson,  "Validation  of  Hardware, 

Software  Systems" 

* R.  T.  Yeh,  "Software  System  Design  Methodology" 

* R.  M,  Walser,  Overview  of  "Basic  Solid  State 
Materials  Research",  and  "Research  on  Solid 
State  Infrared  and  Optical  Devices" 

* R.  W.  Beng,  "Schottlcy  Barriers  in  Transition 
Metal-Silicon  Systems" 

* A.  B.  Buckman,  "Photolysis  in  Pbl^" 

* R.  M.  Walser,  "Recrystallization  of  Ion-Implanted 
Amorphous  Silicon" 

* E.  J.  Powers,  Overview  of  "Atomic  and  Molecular 
Electronic  Processes,"  "Laser,  Holographic,  and 
Optical  Systems,"  "Linear  and  Nonlinear  Wave 
Phenomena  in  Optics  and  Plasmas,"  and  "High 
Energy  Gas  Laser  Systems" 

* M.  Fink,  "Electron  Scattering  from  Atoms  and 
Molecules" 

* M.F,  Becker,  "Infrared  Nonlinear  Optics  in 
Molecular  Gsses" 
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1977  March  Meeting  of  the 
Applied  Physics  Society 
San  Diego,  California 
March  21-24 , 1977 

* M.  Fink,  "Can  Electron  Diffraction  Provide 
More  Accurate  Molecular  Charge  Densities 
Than  Current  Theories" 

Fourth  Symposium  on  Computer  Architecture 
Silver  Spring,  Maryland 
March  23-25,  1977 

T.A.  Welch  and  H.L.  Tredennick,  "High-Speed  Buffering 
for  Variable  Length  Operands" 

Symposium  on  Computer  Architecture 
Washington,  D.C. 

March  24,  1977 

G.J.  Lipovski,  "Panel  Session:  Buzz-words 
in  Computer  Architecture , " and  "On  Virtual 
Memories  and  Micronetworks" 

The  Johns  Hopkins  Conference  on 
Information  Sciences  and  Systems 
Baltimore,  Maryland 
March  30  - April  1,  1977 

* Steven  I.  Marcus,  "Modeling  and  Analysis 
of  Poisson  Driven  Markov  Processes" 

Colloquium  on  Decision  and  Control 
University  of  Texas  at  Austin 
Austin,  Texas 
March  30 , 1977 

B,F.  Womack,  "Design  and  Control  of  Feed- 
back Systems  via  Sensitivity  Theory" 
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The  1977  Conference  on 

Information  Sciences  and  Systems 
The  Johns  Hopkins  University 
Baltimore,  Maryland 
March  30  - April  I,  1977 

G,L.  Wise  and  N.C.  Gallagher,  "A  Novel 
Approach  for  the  Computation  of  Chebyshev 
Polynomial  Expansions" 

University  of  Illinois 
Chicago  Circle  Campus 
April  17,  1977 


R,T.  Yeh,  "An  Outline  of  a Design  Methodology 
for  Data  Base  Systems" 

Joint  Services  Topical  Review  on  Informational 
Aspects  of  Decision  and  Control 
Cambridge,  Massachusetts 
May  4-5 , 1977 

* Steven  I.  Marcus,  "The  Use  of  Spherical 
Harmonics  in  Suboptimal  Esitmator  Design" 


Syracuse,  Ne\w  York 
May  17  , 1977 

G,  J.  Lipovski,  "On  Imaginary  Fields,  Token 
Transfers,  and  Floating  Codes  in  Intelligent 
Secondary  Memories" 


NCC 

Dallas,  Texas 
June  16 , 1977 


G,  J.  Lipovski,  "An  Organization  for  Optical 
Links  Bet\ween  Integrated  Circuits" 
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PUBLISHED  ABSTRACTS 

E.  J.  Powers  and  T,C.  Simonen,  "Microwave  Measurement  of  Plasma 
Density  Fluctuation  Amplitudes , " Bulletin  of  the  American  Physical  Society, 
Vol.  21,  No.  5,  1976,  p.  855. 

* E.  I.  Powers  and  Y.  C.  Kim,  "Application  of  Digital  Spectral  Analysis 
Techniques  to  the  Diagnostic  Problem  of  Fluctuation-Induced  Transport," 
Conference  Record  - The  1976  IEEE  International  Conference  on  Plasma 
Science,  IEEE  Catalog  No.  76CH1083-5-NPS , p.  47. 

E.  J.  Powers,  W.  C.  Turner,  and  T.  C,  Simonen,  "Simultaneous  Enhance- 
ment of  Ion  Velocity  Space  Diffusion  and  Ion  Cyclotron  Oscillations  in  a 
Mirror  Confined  Plasma,"  Conference  Record  - The  1976  IEEE  International 
Conference  on  Plasma  Science,  IEEE  Catalog  No.  76CH1083-5-NPS , p.  77. 

* E.  J.  Powers  and  Y.  C.  Kim,  " Bispectral  Wave  Analysis  of  Nonlinear  Wave- 
Wave  Interactions  in  Plasmas,"  Conference  Record  - The  1976  IEEE  Inter- 
national Conference  on  Plasma  Science,  IEEE  Catalog  No.  76CH1083-5-NPS , 
p.  78. 

E.  J.  Powers,  et  al,  "Fluctuation  Spectra  in  the  NASA  Lewis  Bumpy  Torus"  , 
Bulletin  of  the  American  Physical  Society,  Vol.  21,  1976  , p.  1069. 

* Y.  C.  Kim  and  E.  J.  Powers,  "Utilization  of  Bispectral  Analysis  in  Exper- 
imental Studies  of  Nonlinear  Wave-Wave  Interactions,"  Bulletin  of  th^ 
American  Physical  Society,  Vol.  21,  1976,  p.  1081. 

S.  I.  Marcus,  "Estimation  and  Analysis  of  Nonlinear  Stochastic  Systems," 
IEEE  Transactions  on  Information  Theory,  Vol.  IT-22,  July  1976,  p.  503. 

* A,  B.  Buckman,  "Theory  of  an  Efficient  Electronic  Phase  Shifter  Employing 
a Multilayer  Dielectric-Waveguide  Structure."  IEEE  Conference  Records  - 
Abstracts , Second  International  Conference  and  Winter  School  on  Submilli- 
meter Waves  and  Their  Applications , IEEE  Catalog  No.  76-CH-1152-8-MTT, 
San  Juan,  Puerto  Rico,  December  6-11,  1976. 

* Funded  entirely  or  in  part  by  the  Joint  Services  Electronics  Program. 
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* M.  Fink,  "Can  Electron  Diffraction  Provide  More  Accurate  Charge  Densities 
Than  Current  Theories , " Bulletin  of  The  American  Physical  Society,  Vol.  22, 
1977,  p.  306. 
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CONFERENCE  PROCEEDINGS 

* R.  T.  YehandC.  Reynolds,  " Verification  of  Nondeterministic  Programs  , " 
Proceedings  of  the  11th  Asimolar  Conference  on  Circuits,  Systems,  and 
Computers , November  1975. 

* T,  A.  Welch,  "An  Investigation  of  Descriptor  Oriented  Architecture , " 
Proceedings  of  the  Third  Symposium  on  Computer  Architecture , Clearwater, 
Florida,  January  19-21 , 1976,  Cat.  No.  76CH1043-5C,  pp.  141-146. 

* T.  K.M.  Agerwala,  "Timing  and  Priority  Considerations  in  Concurrent 
Computer  Systems , ' Proceedings  of  the  Conference  on  Information  Sciences 
and  Systems . The  Johns  Hopkins  University,  Baltimore,  Maryland,  March 
31  - April  2,  1976. 

G.  L.  Wise,  "Spectral  Densities  of  Markov  Processes."  Proceedings  of 
the  1976  Conference  on  Information  Sciences  and  Systems,  The  Johns 
Hopkins  University,  March  31  - April  2,  1976,  p.  161. 

* R.  T.  Yeh,  " Prospect! ves  on  the  Nature  of  Software  Problems , " Proceedings 
of  the  IFIP  Working  Conference,  Tokyo,  Japan,  April,  1976. 

* S.  I.  Marcus,  "Analysis  of  Poisson  Driven  Bilinear  Systems , " Proceedings 
of  the  1976  Conference  on  Information  Sciences  and  Systems,"  The  Johns 
Hopkins  University,  April  1976,  pp.  411-415. 

* A.  B.  Buckman,ed.,  "1976  Region  V IEEE  Conference  Digest,"  IEEE  Catalog 
No.  76CH1068-6-REG5,  Austin,  Texas,  April  14-16,  1976. 

R.  H.  Flake,  "Qualitative  Theory  in  Estuarine  Water  Quality  Modelling," 
Proceedings  of  the  IEEE  Region  V Conference,  Austin,  Texas,  April  14-16, 
1976. 

* R.  T.  Yeh,  "An  Approached  Program  Verification,"  Proceedings  of  the 
Design  Automation  Workshop,  San  Francisco,  California,  June,  1976. 


* Funded  entirely  or  in  part  by  the  Joint  Services  Electronics  Program, 
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* S.  I.  Marcus,  "Modeling  and  Analysis  of  Linear  Systems  Driven  By  Multi- 
plicative Poisson  White  Noise,"  Proceedings  of  the  NASA  Seminar/Res^rch 
Institute  on  Differential  and  Algebraic  Geometry  for  Control  Engineers, 

June  197 6 . 

* B.  J.  Olufeagba,  R.  H.  Flake,  and  K,  J.  Almquist,  "A  Multiple  Shooting 
and  Sweep  Algorithm  for  Optimal  Point  Controlled  Distributed  Parameter 
Systems  , " Proceedings  of  the  International  Federation  of  Automatic  Control 
Symposium  on  Large  Scale  Systems  Theory  and  Applications,  ed.,  G.  Guarda- 
bassi  and  A.  Locatelli,  Udine,  Italy,  June  16-20,  1976. 

* R.  M.  Walser  and  R.  W.  Bene,  "Solid  Phase  Reactions  at  Transition 
Metal-Silicon  Interfaces."  Proceedings  of  the  International  Conference 
on  Physics  of  Semiconductors,  Rome,  Italy,  August  1976, 

G.  L.  Wise,  "A  Comment  on  the  Estimation  of  Probability  Density  Functions," 
to  appear  in  Proceedings  of  the  Workshop  on  Decision  Informatiori  f^r 
Tactical  Command  and  Control,  Airlie,  Virginia,  September  22-25,  1976, 

G.  L,  Wise  and  J,  B,  Thomas,  "Zero  Memory  Nonlinear  Transformations 
of  Gaussian  Processes,"  (invited  paper)  , Proceedings  of  the  197  6 Interj 
national  Telemetering  Conference,  Los  Angeles,  California,  September 
28-30,  1976,  pp.  136-142. 

N.  C.  Gallagher,  T.  W.  Allen,  and  G.  L.  Wise,  "Fourier  Series  Represen- 
tation for  Polynomials  with  Application  to  Nonlinear  Digital  Filtering, 
Proceedings  of  the  Fourteenth  Annual  Allerton  Conference  on  Circuit  and 
System  Theory,  Monticello,  Illinois,  September  29-October  1,  1976  , 

pp.  211-218. 

G.  L.  Wise  and  N,  C.  Gallagher,  "A  Representation  for  Spherically 
Invariant  Random  Processes."  Proceedings  of  the  Fourteenth  Annual  Allerton_ 
Conference  on  Circuit  and  System  Theory,  Monticello,  Illinois,  September 
29  - October  1,  1976,  pp.  460-469, 

* S.  A.  Szygenda,  N.  Billawala,  and  E.  W.  Thompson,  "A  Data  Structure 
and  Drive  Mechanism  for  a Table-Driven  Simulation  System  Employing 
Multilevel  Structural  Representations  of  Digital  Systems,"  Proceedings  of 
the  International  Conference  on  Software  Engineering,  October  1976. 
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T.  I.  Wagner,  "A  Distribution- Free  Performance  Bound  in  Error  Estimation," 
Proceedings  of  the  IEEE  International  Symposium  on  Pattern  Recognition , 
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A.  INVESTIGATION  OE  TRANSITION  LAYERS  BETWEEN 

CVD-VOg  THIN  FILMS  AND  OXIDE  SUBSTRATES  t 
Professor  R.  M.  Walser,  Professor  R.  W.  Bene,  and  Ron  Remke 


Our  research  on  VO^  thin  films  was  initiated  as  a model 
system  in  which  to  address  fundamental  questions  about  the  behavior  of 
electronic  devices  operating  at  steady  states  far  from  thermodynamic 
equilibrium.  In  VOg  thin  film  resistors  the  dissipative  structure  is  a 
high  current  filament  with  an  abrupt  metal-semiconductor  phase  boundary. 
Our  previous  investigations  of  the  switching  [1]  and  bolometrlc  response 
[2]  of  these  filaments  revealed  several  anomalies  indicating  that  the 
phase  boundary  was  pinned  for  small  thermal  perturbations.  Phase  bound- 
ary pinning  is  expected  to  significantly  alter  the  interpretation  of  the 
many  previous  studies  of  electronic  and  optical  properties  of  VO^  thin 
films  and  devices  [3]. 

Our  research  in  the  past  year  has  been  directed  toward  deter- 
mining the  origin  of  this  pinning.  The  results  obtained  conclusively  show 
that  the  observed  pinning  is  due  to  an  interfacial  layer  between  the  chem- 
ically deposited  VO  films  and  the  oxide  substrates,  i.e.  , sapphire 
(AlgOa)  quartz  (SiOg).  Our  findings  and  their  interpretations  are 

summarized  in  a forthcoming  journal  article  [4]. 

Concentration  profiles  obtained  from  AES  analysis  show  that 
substrate  metal  atoms  (A1  from  sapphire  and  Si  from  quartzl  penetrate 
the  VOg  film  a distance  lOOOA  from  the  interface.  The  large  depth 
of  inter-diffusion  is  surprising  considering  the  stability  of  the  bulk  metal- 
oxide  bonds  of  these  systems  and  the  low  (500°G)  reaction  temperature. 
Similar  low  temperature  interfacial  reactions  are,  however,  commonly 
observed  at  metal-semiconductor  interfaces.  Because  the  68  C metal- 
semiconductor  phase  transition  is  sensitive  to  small  concentrations  of 
metal  impurity  atoms,  but  relatively  Insensitive  to  stress,  one  can  de- 
termine the  electronic  properties  of  both  transition  regions  by  comparing 
the  measured  temperature  dependence  of  the  VO^  film  resistance  as  a 
function  of  film  thickness. 


+Thls  research  was  supported  entirely  by  the  Joint  Services 
Electronics  Program  under  Contract  F44620-76-C-0089 . 
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In  summary,  the  AES  data  show  that  while  both  systems  have 
chemical  transition  regions,  the  interface  of  VOg-AlgQj  is  electrically 
inactive  compared  with  that  of  VOg-SlO^.  The  data  shows  that  the  transi- 
tion temperature  of  the  first  system  is  identical  to  bulk  VO^  and  inde- 
pendent of  film  thickness  while  that  of  the  latter  varies  as  shown  in  Fig. 

1.  The  electrical  inactivity  of  the  VO2-A1^03  implies  that  AI  must  be 
segregated  during  the  crystallization  of  VOg  and  that  the  transition  re- 
gion has  a much  lower  conductance  than  that  of  semiconducting  VOg  . 

The  opposite  is  true  for  the  VO^-SiOg  interface;  Si  is  evidently  incor- 
porated in  the  VO^  crystalline  grains  and  the  interfacial  region  has  a 
conductance  comparable  with  VO^  and  a higher  transition  temperature. 

The  existence  of  the  transition  region  has  serious  implications 
for  understanding  switching  and  phase  boundary  dynamics  in  VOg  struc- 
tures with  film  thicknesses  comparable  with  the  transition  region.  For 
example,  in  VOg-SiO^,  the  dominant  conduction  path  before  thermal 
switching  wilt  be  along  the  transition  region,  shifting  to  the  surface  after 
the  switch . occurs . Thus,  the  nucleation  of  the  phase  boundary  will  be 
thermally  driven  by  the  subsurface  region,  but  the  easy  growth  path  of  the 
filament  will  be  along  the  surface.  Clearly  the  switching  dynamics  of  the 
phase  boundary  in  this  system  will  be  quite  different  than  those  of  the 
VOg-Alg03  system  where  filament  nucleation  and  growth  should  both  occur 
on  the  surface. 

To  further  test  our  hypothesis  on  impurity  segregation  in  these 
systems,  thin  deposited  films  of  Al  and  Si  were  diffused  into  the  as- 
deposited  amorphous  VOx  film  and  into  polycrystalline  Vg  Os  after 
oxidation-crystallization  of  the  a-VO,  . From  the  shift  in  T^  relative  to 
that  of  the  bulk  crystal,  we  conclude  that  neither  Al  or  Si  can  diffuse 
into  crystalline  VgOs  or  VOg  at  500‘’C,  but  both  are  incorporated  into 
the  crystalline  phase  if  present  in  the  amorphous  phase.  This  implies 
that  the  Al  in  the  AES  profile  of  VOg-AlgO^  was  not  in  the  metallic  state 
during  the  reaction  cycle. 

In  principle,  the  existence  of  a transition  layer  provides  an 
explanation  for  the  anomalies  we  observed  earlier.  These  are  discussed 
in  more  detail  in  our  forthcoming  paper. 
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In  summary,  we  expect  to  detect  a transition  from  interface- 
dominated  to  surface-dominated  effects  in  relatively  thick  (»500  A)  VO^ 
films  on  SiO,  substrates  at  high  bias  levels.  This  transition  is  shown 
in  Fig,  2 for  the  bias-dependent  black-body  radiance  profiles  of  the  fil- 
ament and  in  Fig,  3 for  the  bias  dependence  of  the  spatial  profiles  of  the 
bolometric  response  of  the  filament. 
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Infrared 

Radiance 


Figure  2,  Spatial  Infrared  Radiance  Profiles  Of  *[he  Filament  As  A Function 
Of  Bias . 
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Figure 
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3.  Spatial  Profiles  Of  Bolometric  Response  Of  The  Filament  As  A 
Function  Of  Bias . 
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B.  SOLID-SOLID  INTERFACE  REACTIONS  t 

Professor  R.  W.  Ben6,  Professor  R.  M,  Walser,  W.  J,  Schaffer, 
lames  Hu,  Mark  Chonko,  and  Dan  Richter 


The  purpose  of  this  work  is  to  determine  the  relationship  of 
thermodynamic  equilibrium  parameters  to  interphase  reactions  in  thin 
films  at  low  temperatures,  where  it  appears  that  the  new  phases  formed 
are  determined  kinetically. 

Based  upon  experiments  on  compound  nucleation,  surface  recon- 
struction, Schottky  barrier  heights  , recrystallization  , etc.,  we  have 
postulated  a model  [1,  2,  3]  for  the  dynamics  of  phase  formation  in  Si- 
metal  and  Ge-metal  systems.  A prominent  feature  of  this  model  is  a 
predicted  glassy  membrane  region  between  crystalline  phases.  These 
glass-like  interphases  are  expected  to  significantly  affect  the  kinetics 
of  thin  film  phase  formation,  diffusion,  and  the  static  properties,  such 
as  Schottky  barrier  heights,  in  metal-semiconductor  systems.  In  an 
attempt  to  verify  the  existence  of  these  membrane  regions  and  also  to 
study  the  phase  nucleation  in  thin  films  in  more  detail  concerning  param- 
eter variations,  we  have  begun  measurements  on  systems  of  thin  metals 
(5-1000  k)  sputtered  onto  single  crystal  silicon  substrates.  The  measure- 
ments we  are  doing  at  present  are  x-ray  diffraction  and  transmission 
electron  diffraction  (TED)  , and  will  be  expanded  to  include  Auger  spec- 
troscopy and  Schottky  barrier  measurements  in  the  near  future. 

Transmission  electron  diffraction  measurements  are  presently 
underway  on  thin  metal  films  of  Pd,  Pt,  Ni  and  Co  on  (100)  Si.  We  are 
doing  two  types  of  experiments;  one  is  to  deposit  films  ranging  from 
=“8  ito  100  A on  Si  substrates  and  then  jet-thin  through  the  Si  wafer 
to  prepare  the  TED  samples,  and  the  other  is  to  prepare  relatively  thick 
films  (~  1000  k)  and  then  do  a double  jet-thinning  operation  from  each 
side  to  the  interphase  region  (near  the  Si-silicide  interphase)  . Prelim- 
inary measurements  on  the  Pd-Si  system  have  indicated  that  we  get  an 
amorphous  layer  for  Pd  films  less  than  about  15  k and  then  polycrystalline 
PdaSi  (288  H)  formation  for  films  thicker  than  this.  Figures  1 and  2 show 


t This  research  was  supported  in  part  by  the  Joint  Services 
Electronics  Program  under  Contract  F44620-76-C-0089 , and  by  the 
Office  of  Naval  Research, 
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typical  TED  patterns  for  25  A and  12  I sputtered  films  of  Pd  on  (100) 
Si.  Both  samples  were  unannealed.  More  extensive  measurements  may 
indicate  that  we  are  indeed  seeing  the  membrane  region  for  the  thinner 
films . 


The  Pt,  Ni,  and  Co  systems  also  form  silicides  upon  deposition 
(without  annealing)  as  indicated  by  TED  and/or  x-ray  diffraction.  These 
systems  are  more  complicated  than  the  Pd  case  because  they  each  form 
more  than  one  silicide  in  the  low  temperature  regime  and  so  we  are  not 
as  far  along  in  our  investigations  for  these  systems.  In  addition,  the 
double  etching  process  is  difficult  in  the  Pt  system  because  the  etches 
tend  to  undercut  the  Pt  by  etching  the  silicides  faster  than  Pt.  In  the 
Ni-Si  system,  we  have  observed  formation  of  the  NiaSi(6)  phase  upon 
deposition  using  TED  (not  visible  in  x-ray)  and  also  NiSi  using  x-ray 
diffraction.  We  have  not  as  yet  sorted  out  the  morphology  of  these  poly- 
crystalline regions,  nor  have  we  looked  for  amorphous  regions. 

One  of  the  things  that  needs  to  be  studied  in  more  detail  is  the 
role  of  interfacial  oxygen  on  the  interphase  region.  In  this  regard.  Figs. 
3,4,  and  5 show  TED's  of  a 700  A film  of  Co  on  Si  annealed  at  450°  C 
for  30  minutes,  which  has  subsequently  been  double  jet-thinned  to  expose 
the  interphase  region.  Figure  3 shows  polycrystalline  CogSi  formation 
and  faint  rings  of  CogSiO^  . Figure  4 shows  another  region  in  the  inter- 
phase where  the  CosSi04  rings  are  dominant.  Different  regions  in  the 
interphase  are  brought  out  in  the  TED  due  to  unevenness  of  the  etching 
process.  Finally,  in  Fig.  5 we  show  another  region  in  the  interphase 
where  a very  amorphous  structure  is  obse'^ved. 

We  plan  to  continue  our  efforts  to  characterize  the  interphase 
region  in  these  systems  more  fully  than  has  been  done  in  the  past,  and 
in  addition,  study  the  effects  of  various  parameters  such  as  substrate 
orientation,  substrate  preparation,  sputtering  parameters , substrate 
bias  and  annealing,  on  the  interphase  structure  and  growth  kinetics.  In 
this  regard  we  expect  to  shortly  begin  Auger  and  Schottky  barrier  measure- 
ments to  supplement  the  structural  information  we  are  now  obtaining. 
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Figure  I.  TED  of  25A  Pd  on  (100)  Si 
Showing  Pd^  Si  Formation. 


Figure  2,  TED  of  12  A Pd  on  (100)  Si 
Showing  Amorphous  Rings. 
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Figure  3.  TED  of  Interphase  of  700l  Co  on  (100)  SI 
Showing  Co^  Si  Formation. 


Figure  5.  Same  as  Figure  3 Showing 
Highly  Amorphous  Ring 
Structure . 


Figure  4.  Same  as  Figure  3 Showing 
Co.  Si  Rings . 
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C.  INTEGRATED  OPTICAL  BEAM  ADDRESSER  t 

Professor  A.  B.  Buckmdn 

0 

The  very  large  range  of  refractive  index  (2. 2-2. 7 at  6238  A) 
accessible  by  photolysis  of  Pbl^  films  at  ~ 165“C  raises  the  possibility 
of  multiple  exposure  of  a large  number  of  elementary  holograms  in  guided 
wave  structures  containing  such  films.  As  with  an  ordinary  thick  phase 
hologram,  the  number  of  elementary  holograms  which  can  be  multiply  ex- 
posed increases  with  refractive  index  range.  In  a guided  wave  geometry 
a device  containing  many  superimposed  Bragg  gratings  could  function  as 
a beam  addresser  with  a large  number  of  resolvable  spots,  if  the  mode 
propagation  constant  were  controlled  electrooptica Ily. 

In  the  JSEP-supported  part  of  this  work,  we  have  been  concerned 
primarily  with  design  of  the  experiment  for  exposing  elementary  holograms 
in  PbIg  films,  allowing  measurement  of  diffraction  efficiency.  In  a re- 
lated project  sponsored  by  the  Office  of  Naval  Research,  multiple  layer 
waveguiding  structures  are  being  analyzed  theoretically,  with  special 
attention  given  to  structures  where  small  refractive  index  variations  in 
one  of  the  media  produce  larger  variations  in  mode  effective  index  [1-3]. 

To  record  a hologram  in  hot  Pblg  , random  phase  changes  in  the 
interfering  beams  due  to  air  convection  near  the  heated  film  must  be  elim- 
inated. To  accomplish  this,  a prism  of  triangular  cross  section  was  placed 
in  optical  contact  with  the  glass  substrate.  A diverging  beam  from  an  argon 
laser,  incident  on  the  opposite  edge  of  the  triangular  prism,  was  ob^served 
to  produce  good  quality  interference  fringes,  of  period  ~4000-5000  A 
(depending  on  prism  geometry  and  refractive  index!  in  photoresist  coated 
on  the  glass  substrate.  The  glass  substrate  forms  part  of  an  enclosed  box, 
which  is  held  at  a temperature  of  165°C  by  means  of  resistance  heating. 

The  box  is  heavily  Insulated,  to  avoid  producing  convection  currents  in 
front  of  the  prism  which  would  destroy  the  interference  pattern.  The  first 
experiments  on  Pbig  films  will  begin  shortly. 

By  measuring  the  diffracted  intensity  of  the  phase  gratings  we 
produce,  we  can  evaluate  diffraction  efficiency.  This  will  be  done  for 
both  guided  and  ungulded  modes,  but  emphasis  will  be  on  the  gulded- 


+ This  research  was  supported  entirely  by  the  Joint  Services 
Electronics  Program  under  Contract  F44620-76-C-0089. 
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wave  geometry.  Results  for  the  guided-wave  case  will  be  interpreted  using 
the  theory  of  guided  waves  in  multilayer  waveguides , including  recent 
contributions  to  the  theory  of  periodic,  stratified  media  [4], 


2. 


3. 
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D.  MECHANISM  FOR  LARGE  REFRACTIVE  INDEX 

CHANGE  IN  PHOTOLOYZED  Pbl^  t 
Professor  A,  B,  Buckman 


The  very  large  (~0.5)  refractive  index  change  in  Pblg  that 
occurs  upon  photolysis  at  150-180°C,  observed  by  us  over  two  years  ago, 
has  been  the  subject  of  a variety  of  experiments  designed  to  determine 
the  mechanism  behind  this  striking  effect,  Ellipsometry  allowed  us  to 
attribute  the  effect  to  dispersion  arising  from  a broad,  photolysis-induced 
band  near  580  nm,  which  shifted  to  shorter  wavelengths  as  photolysis 
proceeded  [1],  However,  the  origin  of  this  band  could  not  be  determined 
from  optical  data  alone. 

The  experiments  used  to  shed  light  on  the  mechanism  behind  the 
large,  photolysis-induced  refractive  index  change  were:  (1)  spectral  de- 

pendence of  photoconductivity,  and  (2)  x-ray  diffraction  on  films  photo- 
lyzed  under  various  ambients  and  passivating  overlayers.  The  experimen- 
tal details  are  reported  in  Ref.  2. 

The  spectral  dependence  of  the  photoconductivity  shows  that 
the  photolysis-induced  band  arises  from  electron  transitions  which  pro- 
duce free  carriers.  The  x-ray  diffractometry  traces  from  films  photolyzed 
in  vacuum,  air,  Og  and  CO^,  and  while  covered  with  a transparent,  pass- 
ivating layer,  show  that  the  large  refractive  index  change  depends  on  the 
following  conditions;  first,  gaseous  iodine  must  leave  the  film,  via  the 
surface  and/or  grain  boundaries;  second,  oxygen  or  CO^  must  enter 
from  the  atmosphere  and  bond  with  the  Pb  left  behind,  inhibiting  agglom- 
eration of  neutral  Pb,  The  defects  thus  formed  must  give  rise  to  the  ob- 
served absorption  band.  These  results  are  described  in  detail  in  Ref.  2. 

V/e  intend  to  investigate  the  spatial  resolution  of  this  effect 
by  forming  phase  holograms  in  the  Pbig  films.  Since  the  above  exper- 
iments clearly  demonstrate  the  photochemical  origin  of  the  large  refrac- 
tive index  change,  and  since  amplitude  hologram  recording  in  PbIg  has 
already  been  demonstrated  [3],  we  expect  the  resolution  to  be  sufficient 
for  phase  holography. 

t This  research  was  supported  entirely  by  the  Joint  Services 
Electronics  Program  under  Contract  F44620-76-C-0089 . 
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A,  DESIGN  OF  TWO-DIMENSIONAL  RECURSIVE 

FILTERS  BY  INTERPOLATION  t 
Professor  J.  K.  Aggarwal  and  H.  Chang 


1,  INTRODUCTION 

Recursive  digital  filtering  of  two-  and  multi-dimensional  signals 
is  an  important  technique  in  the  processing  of  two-  and  multi-dimensional 
data.  However,  the  design  of  two-  and  multi-dimensional  recursive  filters 
is  difficult  due  to  the  fact  that  polynomials  in  two  or  more  variables  may 
not,  in  general,  be  factored  into  lower  order  polynomials.  Specifically, 
the  stability  test  procedure  for  a two-dimensional  filter  is  lengthy  and 
tedious.  Further,  if  the  filter  turns  out  to  be  unstable,  there  is  no  guar- 
anteed procedure  which  may  stabilize  the  unstable  filter.  These  difficul- 
ties can  be  partly  circumvented  by  designing  two-  and  multi-dimensional 
systems  as  one  dimensional  systems.  To  do  this,  we  have  developed  a 
technique  for  rotating  the  frequency  response  of  a separable  filter.  In  the 
present  technique,  transfer  functions  having  rational  powers  of  z are 
introduced  and  realized  by  input/output  signal  array  interpolations. 


2-  MAIN  RESULTS 

The  frequency  transformation  technique  is  one  of  the  traditional 
approaches  in  deriving  desired  frequency  selective  filters  from  simple  pro- 
totype filters.  Extending  the  notion  of  frequency  transformation  to  include 
rational  powers  of  Zi  and  Zs  , one  can  rotate  the  frequency  response  with 
band  contraction.  To  begin  with,  consider  the  transformation 


Zl 


- -4 

= ZiZg 


where  a and  ^ are  integers.  The  transformation  has  the  following 
properties: 


tThis  research  was  supported  in  part  by  the  Joint  Services 
Electronics  Program  under  Contract  F44620-76-C-0089 , and  in  part  by 
NSF,  ENG  74-04986. 
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(1)  It  may  be  thought  of  as  a mapping  from  the  unit  circles  of  the 
(z^  , Zg) -domain  into  the  unit  circles  of  the  (zj^  , z^) -domain. 

(2)  When  p/a^O,  the  transformation  is  causal,  otherwise  it  is 
noncausal. 

(3)  In  the  spatial  domain  it  is  equivalent  to  the  rotation  of  recur- 
sion direction  with  a new  sample  interval. 

(4)  In  the  frequency  domain  the  effect  of  the  transformation  is  a 
rotation  of  the  frequency  response  with  a contraction  of  the 
bandwidth. 

(5)  It  is  clear  that  the  stability  of  the  resulting  filter  is  not  affected 
since  the  transformation  can  be  viewed  as  a transformation  of  the 
delay  unit  which  may  be  specified  by  both  recursion  direction 
and  sampling  interval. 

Since  we  allow  noninteger  delays,  the  signal  values  for  the  trans- 
formed filter  are  not  defined  at  the  new  grid  points,  so  the  transformation 
is  not  readily  applicable  to  standard  rectangular  arrays.  If  one  uses  a 
suitable  interpolation  function  that  can  reconstruct  a continuous  signal 
from  a discrete  signal,  it  is  possible  to  generate  an  interpolated  array 
where  signal  values  are  defined  on  the  new  grid  points.  In  general,  there 
are  many  functions  available  for  the  interpolation.  Specifically,  the  band- 
limited  interpolation  is  employed  in  our  design.  The  band-limited  inter- 
polation function  usually  stretches  from  -»  to  +®.  To  obtain  an  interpo- 
lation function  with  finite  terms,  we  apply  windowing  to  it.  With  all  of 
these  features,  it  is  called  the  interpolated  filter  system  which  consists 
of  input  interpolator,  (z^, Zg)~domain  rotated  filter,  and  output  interpolator. 

The  overall  frequency  characteristic  of  the  interpolated  filter  system  is 
determined  by  the  frequency  characteristics  of  input/output  interpolators 
and  (z^ , Zg) -domain  rotated  filters.  Several  useful  filters  were  actually 
designed  such  as  a circularly  symmetric  low-pass  filter,  nonseparable 
band-pass  filter,  and  a ring-shaped  band-pass  filter  to  prove  the  validity 
of  our  design  technique.  In  the  design  of  these  filters,  we  assumed  that 
the  input  signals  are  given  as  finite  area  arrays  so  that  noncausal  filters 
can  be  realized  through  the  array  reorientations.  Also,  we  employed  the 
zero-phase  response  since  the  linear  phase  is  very  important  in  the  image 
processing.  Actual  filter  characteristics  turned  out  to  be  very  smooth  and 
to  have  excellent  circular  symmetry.  A primary  advantage  of  employing 
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separable  filters  in  the  interpolated  filter  system  is  that  one  may  perform 
two-  or  multi-dimensional  manipulations  by  means  of  a series  of  one- 
dimensional manipulations.  With  this  scheme  the  same  interpolated  fil- 
ter system  may  be  used  iteratively,  together  with  proper  array  reorienta- 
tions, resulting  in  an  appreciable  saving  in  hardware. 

Since  the  theory  and  design  techniques  for  two-dimensional 
digital  filters  are  largely  incomplete,  in  the  next  period  we  will  investi- 
gate these  problems  further  until  we  have  fairly  general  theory  and  design 
techniques  for  those  filters.  Then,  we  hope  to  generalize  the  theory  and 
design  techniques  available  for  one-  and  two-dimensional  digital  filters 
to  cover  the  time-varying  digital  filters. 
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NONPARAMETRIC  DISCRIMINATION  AND  ESTIMATION  t 
Professor  C.  S.  Penrod,  Professor  T.  I.  Wagner, 
and  L,  P.  Devroye 


In  the  nonparametric  discrimination  problem,  a statistician  ob- 
serves a random  vector  X taking  values  in  ]R“,  and  wishes  to  estimate 
the  value  of  an  associated  random  variable  0 which  takes  values  in 
{1,  . . . ,M1  . The  random  variable  9 is  usually  referred  to  as  the  state  of 
the  observation  X,  The  statistician  assumes  no  knowledge  of  the  prob- 
ability distribution  of  (X,0)  other  than  what  may  be  inferred  from  his 
data,  a sequence  (X;^  , 9i)  , . . . , (Xo  , 0j  of  independent  identically  dis- 
tributed random  vectors  with  the  same  distribution  as  (X,0)  , given  by 

(a)  P{  0 = j]  , and 

(b)  given  that  0 = j,  X has  probability  measure 

on  the  Borel  sets  of  IR”*  . 

It  is  also  assumed  that  (X,e)  is  independent  of  the  data, 

A discrimination  rule  is  simply  a specific  method  of  choosing 
an  estimate  of  9 from  the  data  and  the  observed  value  of  X,  If  we  allow 
0 to  denote  the  rule's  estimate  of  0,  we  can  define  the  conditional  prob- 
ability of  error  for  the  rule  with  a given  data  set  by 

L„  = P{0  e 1 (X,,9i) (X„  ,0^)1  . (2) 

The  strong  law  of  large  numbers  assures  that  the  relative  frequency  of 
errors  will  converge  to  La  if  we  use  the  rule  and  the  given  data  to  esti- 
mate the  values  of  9 for  a large  number  of  independent  observations. 
Hence,  La  is  a measure  of  the  performance  of  the  rule  with  a given  data 
set. 


t This  research  was  supported  in  part  by  the  Joint  Services 
Electronics  Program  under  Contract  F44620-76-C-0089  and  in  part  by  the 
Air  Force  Office  of  Scientific  Research  Grant  72-2371, 
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The  nonparametric  estimation  problem  is  quite  similar  to  the 
discrimination  problem  in  its  formulation.  The  principle  difference  is 
that  9 (which  is  now  called  the  parameter  associated  with  X)  is  allowed 
to  take  values  in  IR*"  rather  than  [ 1 , . . . , M} , and  we  now  assume  that 
(X,  0)  is  a random  vector  with  joint  distribution  F(x,  9)  . We  also  assume 
that  a loss  function  L is  g^iven,  defined  on  IR'’  x IR*"  so  that  L(0,  6 ) 
is  the  loss  incurred  when  9 is  the  rule's  estimate  of  0.  We  now  define 
the  conditional  expected  loss  of  a rule  with  a given  data  set  as 

L„  = E[L(9,9)  I (Xi,9i) (X„,0J]  . (3) 

As  in  the  discrimination  problem,  Ln  measures  the  performance  of  the 
rule  with  a given  data  set  since  it  is  the  average  loss  which  will  be  in- 
curred if  the  values  of  9 are  estimated  for  a large  number  of  independent 
observations . 

In  both  the  discrimination  and  estimation  problems,  if  the  under- 
lying distribution  of  (X,  9)  is  known,  it  is  usually  possible  to  construct 
an  optimal  estimate  9*(X)  of  9 (e.g,,  in  the  estimation  problem  one  may 

have  to  be  satisfied  with  an  ®-optimal  estimate) . Such  an  estimate  does 
not  depend  on  the  data  and  is  optimal  in  the  sense  that  its  probability  of 
error  or  expected  loss  is  less  than  or  equal  to  that  of  any  other  rule . The 
optimal  rules  are  known  as  Bayes  rules,  and  their  expected  loss  (or  prob- 
ability of  error)  is  denoted  R*.  R*  then  provides  a standard  of  excellence 
against  which  we  can  compare  the  performance  of  all  nonparametric  rules. 

One  of  the  basic  approaches  to  the  nonparametric  discrimination 
problem  is  suggested  by  the  manner  in  which  an  optimal  Bayes  rule  is  con- 
structed. Let 

Pj  (X)  = P{9  =il  X]  , 1 sj  sM. 

Then  it  is  easily  shown  that  a Bayes  optimal  rule  results  if,  given  X,  we 
choose  9*(X)  = j if 

Pj  (X)  > Pi  (X)  for  all  i j . 
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(Ties  can  be  broken  by  employing  some  appropriate  form  of  randomization.) 
Of  course,  in  our  nonparametric  case,  tfie  functions  Pj  (x)  are  unknown 
since  they  depend  on  the  distribution  (1)  , For  example,  if  the  measures 
(ij  are  absolutely  continuous  with  respect  to  Lebesgue  measure,  then 
there  are  probability  density  functions  fj  , 1 ^ j sM  corresponding  to 
the  M.J  . These  densities  can  be  used  to  define  versions  of  Pj(x)  by 

/ttj  fj  (x)/f(x)  , if  f(x)  > 0 

' jo  , if  fM  =0, 

M 

where  f(x)  fj  (x)  . This  suggests  that  the  nonparametric  problem  could 

1 

be  approached  by  using  the  data  to  estimate  the  density  functions,  and  then 
using  the  density  estimates  in  place  of  the  true  densities  in  the  manner 
described  in  (4).  Such  rules  are  known  as  "two  step"  rules.  In  Devroye 
and  Wagner  [1]  the  results  of  an  extensive  study  of  density  estimation  and 
two  step  rules  are  presented.  Conditions  are  shown  for  which  consistent 
density  estimates  yield  rules  for  which  Lc  - R*  in  probability  or  with 
probability  one.  In  addition,  a detailed  study  is  made  of  the  Parzen-Ros- 
enblatt  [2,  3]  and  Loftsgaarden-Quesenberry  [4]  density  estimates,  ex- 
tending many  of  the  previously  known  results. 

Another  approach  to  the  nonparametric  discrimination  problem 
is  the  familiar  nearest  neighbor  rule.  This  rule  simply  chooses  0 - 9i 
if  Xi  is  the  closest  observation  to  X from  the  data^.  Various  generali- 
zations include  the  k nearest  neighbor  rule  where  9 is  chosen  to  be 
the  state  which  occurs  most  frequently  among  the  k nearest  neighbors 
to  X.  Among  the  well-known  propetues  of  nearest  neighbor  rules  is  the 
fact  that  if  k is  allowed  to  depend  on  n in  such  a way  that  k ^ ® and 
k/n  ^ 0,  then  it  is  possible  to  have  L„  R*  (Fix  and  Hodges  [5])  . In 
(1)  this  result  is  extended  to  a more  general  class  of  nearest  neighbor 
rules  which  weight  the  influence  of  the  closer  observations  more  heavily. 

Nearest  neighbor  rules  have  also  been  used  as  solutions  to  the 
nonparametric  estimation  problem  (Cover  [6]) . The  k nearest  neighbor 
rule's  estimate  of  9 is  obtained  by  averaging  the  parameters  associated 
with  the  k closest  observations  to  X from  the  data.  Cover  shows  con- 
ditions for  which 

ELb  -(1  + l/k)R* 
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for  k nearest  neighbor  rules  with  squared-error  loss  functions.  In  Penrod 
and  Wagner  [7] , these  results  are  extended  to  show  that 

Lj  - (1  + l/k)R*  in  probability. 

The  significance  of  this  result  is  that  it  concerns  the  convergence  of  L, 
for  a single  long  data  sequence  as  opposed  to  the  convergence  of  the  aver- 
age value  of  La  averaged  over  a large  number  of  long  data  sequences. 

In  addition,  conditions  are  given  for  which  Ln  -R*  jin  probability  for  k 
nearest  neighbor  rules  which  have  k ^ and  k/n  0. 

The  results  mentioned  thus  far  have  all  been  asymptotic  results 
concerning  the  performance  of  various  rules  with  large  data  sets.  For 
obvious  reasons,  however,  the  designer  of  a nonparametric  discrimination 
or  estimation  rule  is  always  interested  in  knowing  the  value  of  for 
this  rule  with  the  data  available  to  him.  But,  due  to  the  fact  that  it 
cannot  be  computed  without  knowledge  of  the  underlying  probability  dis- 
tribution of  (X,9  ) , which  is  not  available,  he  must  be  content  with  an 
estimate  of  the  true  value  of  Lg.  In  (1)  the  properties  of  three  different 
estimates  of  Lg  for  discrimination  rules  are  studied.  These  estimates 
are  known  in  the  literature  (see  Toussaint  [8])  as  the  deleted  estimate, 
the  holdout  estimate,  and  the  resubstitution  estimate.  For  linear  discrim- 
ination rules,  local  rules,  and  two  step  rules,  conditions  are  shown  for 
which  P(  iLg  - Lg  1 ^e),  where  Lg  is  one  of  the  above  estimates  of  Lg, 
can  be  bounded  independently  of  the  distribution  of  (X,  0) . These  results 
allow  the  designer  of  a rule  to  construct  confidence  intervals  for  the  per- 
formance of  his  rule.  These  results  are  extended  to  nonparametric  esti- 
mation rules  in  [7].  The  report  [7]  also  contains  the  results  of  an  exten- 
sive simulation  study  of  the  performance  of  deleted  and  holdout  estimates 
of  Lg  with  nearest  neighbor  discrimination  and  estimation  rules. 


References 

1.  L.  P.  Devroye  and  T,  J.  Wagner,  "Nonparametric  Discrimination 

and  Density  Estimation,"  Technical  Report  No,  183,  Electronics 
Research  Center,  University  of  Texas  at  Austin,  1976. 


I 


(II.  INFORMATION  SYSTEMS) 


E,  Parzen,  "On  the  Estimation  of  a Probability  Density  Function 
and  the  Mode,"  Annals  of  Mathematical  Statistics,  Vol.  33, 
1962,  pp.  1065-1076. 

M,  Rosenblatt,  "Remarks  on  Some  Nonparametric  Estimates  of 
a Density  Function,"  Annals  of  Mathematical  Statistics,  Vol. 

27,  1957,  pp.  832-837. 

D.  O.  Loftsgaarden  and  C.  P.  Quesenberry,  "A  Nonparametric 
Estimate  of  a Multivariate  Density  Function,"  Annals  of  Math- 
ematical Statistics,  Vol.  36,  1965,  pp.  1049-1051. 

E.  Fix  and  J.  L.  Hodges,  " Discriminatory  Analysis . Nonpara- 
metric Discrimination:  Consistency  Properties,"  Report  4, 
Project  No.  21-49-004,  United  States  Air  Force  School  of  Avia- 
tion Medicine,  Randolph  Field,  Texas,  1951. 

T.  M.  Cover,  "Estimation  by  the  Nearest  Neighbor  Rule,"  IEEE 
Transactions  on  Information  Theory,  Vol.  IT-14,  January  1968, 
pp.  50-55. 

C.  S.  Penrod  and  T.  J.  Wagner,  "Nonparametric  Estimation  with 
Local  Rules,"  Technical  Report  #182,  Electronics  Research  Cen- 
ter, University  of  Texas  at  Austin,  1976. 

G.  T.  Toussaint,  "Bibliography  on  Estimation  of  Misclassifi- 
cation,"  IEEE  Transactions  on  Information  Theory,  Vol.  IT-20, 
July  1974,  pp.  472-479. 


(II  . INFORMATION  SYSTEMS) 


C,  A COMPUTER  ARCHITECTURE  FOR  BYTE  STRING 

AND  VECTOR  OPERATIONS  t 
Professor  T.  A.  Welch  and  H.  L.  Tredennick 


1.  INTRODUCTION 

This  report  presents  the  findings  of  a computer  design  exper- 
iment carried  out  at  The  University  of  Texas  at  Austin.  A computer  arch- 
itecture for  variable  word  length  processing  is  proposed  to  support  high- 
level  languages.  The  proposed  architecture  uses  byte  iterative  processing 
and  specialized  hardware.  The  specialized  hardware,  shown  in  Fig.  1, 
consists  of  a byte  serial  ALU  and  a separate  address  generator  and  mem- 
ory buffer  for  each  operand  byte  stream. 

Variable  word  length  processing  has  been  investigated  to  de- 
termine which  subtasks  are  difficult  and  which  are  easy  and  to  determine 
the  feasibility  of  variable  word  lenth  processing.  Additionally,  require- 
ments and  implementation  strategies  for  the  ALU  , operand  address  gener- 
ation, and  memory  buffering  have  been  analyzed. 

Use  is  made  of  cheaper  and  more  complex  logic  and  memory 
componenets  to  effect  a more  complex  design  for  simpler  programming. 

The  processor  is  designed  to  manipulate  more  complex  operands  rather 
than  provide  more  complex  operations.  The  processor  is  designed  for 
processing  operands  which  might  be  unequal  length  elements,  arrays, 
or  arbitrary  data  structures.  Design  centers  on  processing  for  character 
strings  and  variable-length  arithmetic,  but  scientific  processing  is  not 
ignored. 


2.  VARIABLE  FIELD  LENGTH  PROCESSING 

The  following  problems  are  associated  with  variable  field 
length  processing: 

(1)  It  is  impractical  to  store  either  complete  variable-length 
operands  or  partial  results  in  registers. 

tThis  research  was  supported  in  part  by  the  Joint  Services 
Electronics  Program  under  Contract  F44620-76-C-0089,  and  NSF  Grant 
GI-42514. 
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(2)  Variable-length  remove  and  replace  are  difficult  in  a fixed 
\word  length  memory  with  random  operand  placement. 

(3)  Operand  alignment  may  be  required  for  unequal  length 
operands . 

(4)  Simultaneous  operand  fetching  is  desirable  in  an  iterative 
solution  which  uses  multiple  operand  streams. 

(5)  Sign  and  exponent  placement  are  significant  problems  due 
to  the  serial  scan  nature  of  processing. 

(6)  Optional  scan  direction  is  desirable. 

(7)  Sign,  zero,  or  blank  extension  is  desirable  for  processing 
unequal  length  operands. 


3.  VARIABLE  LENGTH  ARITHMETIC 

Iterative  processing  is  preferred  over  length  limitation  because 
limiting  operand  lengths  to  the  maximum  register  or  accumulator  length 
was  felt  to  be  too  restrictive.  Byte  serial  processing  was  selected  over 
bit  and  nibble  processing  because  the  narrower  data  paths  for  bit  and 
nibble  transfers  slow  processing.  Data  lengths  greater  than  a byte  rein- 
troduce many  problems  (such  as  shifting , masking  , and  alignment)  and 
ALU  costs  increase  rapidly  with  width. 

The  two  viable  choices  for  number  representation  and  scan 
direction  are:  (D  a right-to-left  scan  with  a two's  complement  number 
representation  or,  (2)  a left-to-right  scan  with  a sign  and  magnitude 
number  representation.  The  two's  complement  representation  and  right- 
to-left  scan  were  chosen  because  of  greater  ease  and  efficiency  of  the 
addition  scan. 

Character  and  integer  instruction  implementations  are  simple 
and  multiply  and  divide  instruction  implementations  are  difficult.  The 
character  instructions,  such  as  move,  compare,  shift,  and  Boolean, 
and  the  integer  instructions,  such  as  add,  subtract,  move,  compare,  and 
replace  all  require  only  a simple  single  pass  scan.  The  multiply  and 
divide  instructions  require  multiple  scans  and  complex  scanning  sequences 
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4.  ALU  IMPLEMENTATION 

Implementation  of  a byte  serial  ALU  capable  of  supporting  the 
basic  character  and  integer  operations  would  require  about  52  IC's.  The 
incremental  increase  in  ALU  hardware  required  to  support  floating-point 
operations  is  about  50%.  An  increase  of  about  50%  over  hardware  for  the 
basic  ALU  will  support  multiply  operations.  Additional  hardware  to  sup- 
port the  divide  instruction  is  about  an  80%  increase  over  hardware  to 
support  the  multiply.  A complete  implementation  requires  about  160  IC's. 

Time  to  perform  a 32-bit  byte  serial  add  is  about  480  gsec  com- 
pared to  about  160  ^pec  for  a 32-bit  parallel  add.  Performance  of  the  add 
in  a byte  serial  manner  represents  about  a 20%  degradation  in  the  total 
instruction  execution  time.  Execution  of  a 32-bit  multiply  in  byte  serial 
fashion  requires  about  a sixteenth  of  the  hardware  for  a parallel  multiply 
and  is  about  sixteen  times  slower.  A byte  serial  multiplier  is  about  equal 
in  speed  and  cost  to  a bit  serial  implementation.  A byte  serial  divide 
implementation  is  about  the  same  cost  as  a fixed  length  divide,  but  is 
about  half  as  fast. 


5.  OPERAND  ADDRESSING 

All  of  the  problems  associated  with  variable-length  processing 
affect  the  design  of  the  operand  address  generators.  A three  address 
instruction  format  of  the  form  A op  B -•  C was  selected  for  the  pro- 
cessor design  and  a separate  address  generator  was  provided  for  each 
operand  stream. 

The  operand  address  generator  provides  the  following  element 
addressing  features; 

(1)  Optional  element  scan  direction 

(21  Terminate  on  short  element  (Truncation) 

(3)  Repeat  short  element 

(4)  Extension  of  leading/trailing  byte/zeroes 

(5)  Element  alignment 
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The  operand  address  generator  provides  the  following  vector 
addressing  features; 

(1)  Skip  internal  fields 

(2)  Vector  chaining 

(3)  Repeat  vector 

(4)  Repeat  last  element 

(5)  Terminate  on  short  vector 

(6)  Optional  vector  scan  direction 

The  Operand  Address  Generators  can  provide  byte  addressing 
for  variable  word  length  operands,  addressing  for  simultaneous  byte 
streams,  operand  alignment,  and  operand  extension.  In  addition,  the 
address  generators  can  provide  complex  addressing  sequences  required 
to  implement  complicated  iterative  processes,  and  support  array  operands 
with  interspersed  data  or  concatenated  fields.  Implementation  of  the 
address  generators  to  support  these  operands  and  operations  requires 
complex  control,  many  data  paths,  and  much  interprocess  communication. 
Each  Operand  Address  Generator  requires  about  90  IC's  for  implementation 
(using  a mix  of  standard  and  Schottky  TTL'  . 


6.  VARIABLE  FIELD  LENGTH  BUFFERING 

Problems  confronted  in  design  of  the  memory  buffer  include  the 
lack  of  variable-length  registers,  accessing  of  variable-length  operands, 
operand  alignment,  and  provision  of  simultaneous  data  streams. 

The  following  observations  can  be  made  about  the  requirements 
for  implementing  the  buffer:  (1)  It  must  be  fast  enough  to  compensate  for 

the  lack  of  variable-length  registers;  (2)  it  must  be  large  enough  to  hold 
several  operands  of  typical  size;  (3)  automatic  buffer  allocation  is  re- 
quired because  allocation  is  too  difficult  and  dynamic  for  the  programmer 
to  manage;  (4)  dynamic  byte  sequence  control  is  desirable  for  efficiency; 
(5)  the  ALU  requires  a byte  interface  and  the  memory,  a word  or  block 
interface;  (6)  there  should  be  a short  initiation  latency  (time  from  termi- 
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nation  of  one  process  to  the  beginning  of  the  next  should  be  short);  and 
(7)  multiple  access  paths  are  necessary  for  efficient  processing  using 
simultaneous  data  streams  in  an  iterative  process. 

Four  schemes  were  considered  for  implementation;  a single 
word  buffer,  a single  cache  buffer,  a separate  queued  buffer  for  each 
operand  stream,  and  a separate  cache  buffer  for  each  operand  stream. 

Single  Word  Buffer.  The  single  word  buffer  has  the  advantages 
of  simplicity,  limited  data  currency  problems,  and  short  start-up  time. 

It  has  the  disadvantages  of  being  slow,  not  allowing  block  transfers,  and 
holding  only  a single  word  of  previous  results  (for  retrieval)  . 

Single  Cache.  The  single  cache  buffer  has  the  advantages  of 
allowing  main  memory  block  transfers,  providing  some  register-like  cap- 
abilities, being  faster  than  the  one  word  buffer,  and  saving  main  memory 
cycles.  The  single  cache  suffers  from  access  priority  selection  problems 
and  requires  sequential  access.  Conventional  cache  design  improvements 
are  of  minimal  avail  because  they  are  designed  to  reduce  access  conflicts 
for  space.  In  an  iterative  processor  employing  simultaneous  data  streams, 
access  conflicts  in  time  are  much  more  important  than  conflicts  for  space 
due  to  the  parallel  nature  of  the  data  streams. 

Queued  Buffer.  Separate  queues  for  each  operand  stream  provide 
simultaneous  operand  streams,  simplify  buffer  allocation,  and  can  be  fast 
enough  to  support  a byte  serial  ALU.  The  queued  buffering  scheme  has 
many  disadvantages:  there  may  be  multiple  copies  of  data,  queue  searches 
are  difficult,  queue  control  is  complex,  the  queues  require  much  data  and 
address  storage,  queues  are  inflexible,  and  initialization  and  emptying 
are  slow  processes. 

Separate  Cache  Buffers.  The  separate  cache  buffering  scheme, 
shown  in  Fig.  2,  provides  the  following  advantages:  virtual  variable- 
length  registers , hardware  controlled  allocation,  byte  addressing  for  the 
ALU  byte  streams  and  block  addressing  for  main  memory,  and  simultaneous 
data  paths  for  the  separate  data  streams.  It  has  the  disadvantages  that 
there  may  be  multiple  copies  of  the  data,  many  address  paths  and  com- 
parisons are  required,  complex  control  functions  are  required,  and  per- 
verse address  sequences  exist  which  defeat  the  allocation  scheme 
(direct  mapping) . 
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7.  CONCLUSIONS 

The  proposed  architecture  has  been  evaluated  and  shows  that 
byte  serial  iterative  processing  of  variable-length  operands  is  feasible. 
Variable  word  length  operations  are  significantly  faster  on  the  proposed 
machine  than  on  any  fixed  word  length  machine  in  its  class(minicomputers'  . 
Fixed  word  length  operations  are  faster  than  their  software  counterparts 
on  a fixed  length  machine,  about  as  fast  as  a microprogrammed  implemen- 
tation, and  slower  than  specialized  hardware. 

The  MSI  and  LSI  technologies  allow  sophisticated  ALU  design. 
The  byte  serial  ALU  is  not  the  processing  bottleneck. 

Operand  byte  access  is  the  limiting  factor.  Stores  are  slower 
than  reads  and  multiple  stores  may  be  required. 
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Figure  1.  Variable  Length  Data  Processors. 
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Figure  2.  Separate  Cache  Memory  Scheme. 
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A.  SIMULATION  AND  RELIABILITY  OF  FUNCTIONAL  DEVICES  t 
Professor  S.  A.  Szygenda,  Professor  E.  W,  Thompson,  and  A.  K.  Bose 


Simulation  of  logic  nets  has  been  used  extensively  for  the  veri- 
fication of  logic  design  and  to  study  the  reliability  of  logic  systems.  The 
level  of  simulation  however,  has  been  restricted  to  simple  Boolean  gates 
and  flip-flops.  This  restriction  has  not  posed  any  severe  problems  to 
users  in  the  past,  since  the  size  of  nets  being  used  has  been  small  and 
the  devices  simple.  With  the  current  trend  towards  design  of  large 
systems  with  the  use  of  complex  devices , logic  simulation  has  become 
inadequate.  Simulation  at  a functional  level  provides  a possible  solu- 
tion to  this  problem. 

This  report  discusses  several  aspects  of  this  trend.  Major 
problems  that  remain  unsolved  in  the  area  of  functional  level  simulation 
are  discussed  and  solutions  offered  for  some  of  them.  Two  major  func- 
tional devices  are  examined  and  ways  to  improve  their  reliability  are 
suggested . 

The  lack  of  understanding  that  exists  for  functional  devices 
indicates  that  functional  level  simulation  has  to  evolve  from  gate  level 
simulation.  Moreover,  since  the  functional  description  of  a system  is 
not  adequate  to  perform  an  in-depth  analysis,  such  as  prediction  of  haz- 
ards and  races,  an  interchangeable  functional  and  gate  level  description 
of  systems  is  suggested  [1].  Such  a capability  would  provide  for  efficient 
macroscopic  simulation  by  using  the  functional  description,  or  accurate 
timing  and  failure  analysis  by  using  the  gate  level  description. 

Under  certain  situations,  it  may  be  possible  to  simulate  a func- 
tional device  as  accurately  as  a gate.  This  can  be  realized  by  using  a 
data  structure  of  the  type  shown  in  Table  1 to  represent  multivalued  sig- 
nals at  the  input  of  a system  [2]. 


tThis  research  was  supported  in  part  by  the  Joint  Services 
Electronics  Program  under  Contract  F44620-76-C-0089,  and  in  part  by 
Naval  Surface  Weapons  Center  N609-76-C-A168. 
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Table  1.  Signal  Representation. 


Signal  Description 


Current 

Value 


Future 

Value 


Logical  1 


Logical  0 


Upward  Transition 


Downward  Transition 


Unknown  Signal 


Using  this  data  structure  and  the  known  values  of  the  input 
signals,  the  following  binary  vectors  may  be  defined; 

(1)  The  current  input  vector: 

X = <X.  , X„  , . . . X > 

c Ic  2c  nc 

where  x.  is  the  current  value  of  the  input  signal  x, . 

1C  1 

(2)  The  future  input  vector: 


Xp  = <x^p,  X2p,  ...  x^p> 

where  x^^  is  the  future  value  of  the  input  signal  x. . 

(3)  Intermediate  input  vectors: 

''2k'  ••• 


nk  '^ik'-^^ic' 


where  the  vectors  are  generated  by  taking  all  possible 
different  combinations  of  current  values  and  future  values 
of  the  inputs. 
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These  binary  vectors  can  be  used  to  evaluate  a binary  descrip- 
tion of  the  functional  device,  and  the  outputs  obtained  from  the  current 
and  future  input  vectors  (y^,  and  y^  ) can  be  combined  using  Table  1 to 
yield  a multivalued  output.  The  outputs  obtained  from  the  intermediate 
inputs  (yj)  can  be  used  to  predict  hazards  by  executing  the  following 
steps: 

Step  1:  Using  the  intermediate  input  vectors,  generate  inter- 
mediate outputs,  y^j^  = 

Step  2:  If  Y = y„  and  y 7^  y^  = y„  , then  a static  hazard  exists. 

C r 11  C r 

step  3:  If  y^^  Yp  » then  if  all  y^^  are  identical,  no  dynamic 
hazard  exists.  If  y^^.  ^ y^j^  for  any  j and  k,  then  a 
dynamic  hazard  is  present. 

The  analysis  of  faulty  behavior  of  functional  devices,  especially 
those  of  a sequential  nature,  poses  severe  problems.  A single  fault  in 
a serial  shift  register,  a functional  device  that  has  wide  applications 
in  computer  systems,  has  a global  corrupting  influence  on  the  data  prop- 
agating through  it.  The  correction  of  errors  introduced  in  the  data  is 
impossible,  but  the  detection  of  errors  can  be  realized  through  a nominal 
amount  of  redundancy  [3] . 

An  analysis  of  various  shift  register  configurations  shows  that 
the  effect  of  a single  fault  is  always  restricted  to  a single  type  of  data 
bit  (0  or  1).  This  result  indicates  that  errors  in  a shift  register  can  be 
detected  by  counting  the  number  of  ones  entering  the  shift  register  and 
comparing  this  with  the  number  of  ones  coming  out  of  it.  A hardware 
scheme  to  do  this  can  be  realized  by  using  an  up/down  counter.  Since 
a counter  requires  around  loggN  flip-flops  to  count  up  to  N,  while  a 
N bit  shift  register  requires  N flip-flops,  the  redundancy  involved  is  of 
the  order  of  log^N . 

A computer  Input/output  controller  is  a functional  unit  whose 
complexity  is  very  high  (300-400  gates).  Since  these  controllers  are 
used  to  interface  with  electromechanical  Input/output  devices , their 
failure  could  be  catastrophic.  To  prevent  such  a situation,  these  units 
are  often  duplicated,  A detailed  analysis  of  its  organization  reveals 
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that  it  can  be  made  highly  reliable  by  introducing  only  30%  additional 
hardware  [4] . 

These  results  indicate  that  great  cost  savings  can  result  from  j 

a close  analysis  of  functional  devices.  The  often-used  approach  of  i 

duplicating  systems  is  uneconomic  and  unnecessary.  The  only  feasible  | 

way  to  do  the  analysis  is  by  simulating  the  system.  As  simulation  is  an  | 

expensive  process  and  is  required  at  various  levels  of  the  design  phase,  i 

the  capability  to  simulate  at  a functional  level  is  essential.  | 

Since  functional  level  simulation  is  a fairly  new  concept,  the  i 

current  state  of  the  art  is  far  from  satisfactory.  There  are  problems 

associated  with  the  description  of  functional  devices,  simulation  proce-  ! 

dures,  simulation  accuracy  and  fault  handling  strategies  that  need  to  be 

solved  before  functional  level  simulation  is  meaningfully  realized.  | 
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B.  development  of  methodologies  and  models  for  the  design  of 

DIGITAL  SIMULATION  SYSTEMS  WITH  VERY  HIGH  PERFORMANCEt 
Professor  S.  A.  Szygenda  and  Yu-Huei  Jea 

I . INTRODUCTION 

Digital  simulation  is  a mandatory  tool  in  most  of  today's  design 
automation  systems.  An  ideal  simulator  should  satisfy  three  major 
criteria,  namely,  accuracy,  cost-effectiveness,  and  ease  of  use  [1]. 
This  report  outlines  the  methodologies  and  models  developed  to  fit  these 
requirements.  Detailed  descriptions  can  be  found  in  [2]. 


2 . BOOLEAN  VECTOR  E-ALGEBRAS 

Due  to  the  advances  of  technology,  there  is  an  increasing  inter- 
est in  multivalued  logic  systems  where  more  than  two  logical  values  are 
used.  For  example,  the  following  topics  which  are  important  in  the  com- 
puter area,  are  related  to  multivalued  logic:  multivalued  logic  circuits, 
digital  simulation,  test  generation,  multibus  digital  circuits,  etc. 

In  the  application  of  multivalued  logic,  each  logic  value  could 
be  represented  by  a Boolean  vector,  i.e.  , a vector  with  binary  compo- 
nents (0  or  1).  Therefore  it  is  quite  important  to  have  a thorough  under- 
standing of  properties  embedded  in  the  algebraic  structures,  using  Boolean 
vectors  as  basic  operands.  Here  we  briefly  introduce  a new  mode  of 
complement  and  a new  type  of  algebra. 

A P mode  complement  is  a permutation  which  is  not  an  identity 
mapping.  Only  when  a permutation  is  equal  to  the  product  of  disjoint 
transpositions  can  a P mode  complement  have  the  order  of  2. 

An  E- algebra  {s,l,  + ,'}  is  a set  S with  two  binary  operations 
. and  + and  one  unary  operation  ' such  that  for  all  x,y,z  e S we 
have 


El.  x*x=x 

x+x=x 

(Idempotent) 

E2.  x*y=yx 

x+y=y+x 

(Commutative) 

■^This  research  was  supported  in  part  by  the  Joint  Services 
Electronics  Program  under  Contract  F44620-76-C-0089,  and  in  part  by 
the  Naval  Surface  Weapons  Laboratory  N 60921-76-C-A168. 
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E3.  X*  (y  z)=(x- y)  • z x+ (y+z)=(x+y)+z  (Associative) 

E4.  There  exists  0,1  e S such  that 
x-0=0  x+0=x 

X*  l=x  x+l  = l (Universal  Bounds) 

E5.  For  every  x e s,  there  exists  a unique  x'  e S such  that 

x"=x  (Involution) 

E6.  (x-y)'=x'+y' 

(x+y)'=x'»y'  (de  Morgan) 

A Boolean  vector  E-algebra  is  an  E-algebra  {S, . , + ,'1  where  S 
is  a Boolean  vector  set,  . and  + are  two  Boolean  vector  binary  opera- 
tions, and  ' is  a Boolean  vector  unary  operation. 

One  application  of  Boolean  vector  E-algebras  has  been  suggested 
in  the  digital  simulation  area.  It  provides  one  way  to  define  the  informa- 
tion representation  for  signal  models  in  a parallel  simulation  environment. 


3.  VARIOUS  APPROACHES  TO  REPRESENTING  SIGNAL  MODELS  AND 
RELATED  LOGICAL  OPERATIONS 

In  the  design  of  a digital  simulator,  one  of  the  first  decisions 
to  be  made  is  to  determine  how  many  signal  states  we  want  to  propagate. 
Then  we  must  decide  the  information  representation  of  these  signal  values. 
Also  the  implementation  of  related  logical  operations  has  to  be  settled. 
Usually  the  good  representation  of  operations  depends  on  the  proper 
representation  of  operands,  and  vice  versa. 

Four  major  approaches  can  be  used  in  various  situations: 

(1)  Table  look-up  method 

(2)  Switching  theory  approach 

(3)  Relational  approach 

(4)  Boolean  vector  E- algebra  approach 
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Although  global  optimization  of  information  representation 
might  be  quite  difficult  or  even  impossible,  near  optimal  representations 
can  be  achieved  through  the  optimal  implementation  of  those  types  of 
elements  which  occur  most  frequently,  in  most  digital  networks. 


4.  GATE  MODELING  AND  TIMING  ANALYSIS  ALGORITHM 

i Gates  are  the  basic  elements  of  digital  circuits.  They  may  have 

I several  inputs  but  only  one  output.  Two  operations  are  normally  per- 

j formed  by  them: 

(1)  V-operation  (or  value  operation) 

(2)  D-operation  (or  delay  operation) 

The  V-operation  of  any  gate  is  usually  deterministic  in  nature  and 
clearly  defined  by  the  corresponding  operation  table.  The  complexity 
of  the  D-operation,  however,  depends  on  the  delay  model  adopted.  It 
i can  be  as  simple  as  an  identity  mapping  in  a zero  delay  model,  or  as 

complicated  as  a composition  of  the  following  four  mappings  in  a precise 
j delay  model:  P-mapping  , S-mapping,  A-mapping,  and  T-algorithm. 

1 T-algorlthm  (timing  analysis  algorithm)  is  not  simple  in  the 

t three- valued  model  [3,4]  and  even  more  complicated  in  the  five- valued 

[ model.  In  this  research,  we  not  only  find  a five-valued  T-algorithm, 

[ but  also  prove  its  accuracy  and  functional  completeness. 

) 

f There  are  four  types  of  delay  parameters  used  in  the  precise 

delay  model:  djj^  (min.  rise  delay),  d„|,  (max.  rise  delay),  d„p  (min. 
fall  delay),  and  d,^-  (max.  fall  delay).  Their  usage  in  our  model,  which 
is  different  from  [5],  can  be  summarized  as  follows: 

• Add  d_.  when  the  input  signal  changes  to  U. 

• Add  d„.  when  the  input  signal  changes  to  1 . 

nn 

! • Add  d.  r when  the  input  signal  changes  to  D. 

I ■ ^ 

I • Add  d^p  when  the  input  signal  changes  to  0. 

• Add  min  (d  » , d -)  when  the  input  signal  changes  to  E. 

I ffl  R ' B r 


j 
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5.  HIGH-LEVEL  INPUT  LANGUAGES  FOR  DIGITAL  NETAVORKS 

A system  which  is  difficult  to  use  usually  finds  little  use.  This 
statement  gives  rise  to  the  important  question  of  how  to  represent  the 
digital  network  such  that  the  user  can  easily  transcribe  it  into  a ready- 
for- input  format. 

Postfix  notation,  also  known  as  Polish  notation,  is  a very  con- 
venient method  of  incorporating  both  the  simplicity  of  a functional  descrip- 
tion and  the  preciseness  of  a structured  description.  For  example, 

ABC3+DEF4 

denotes  a circuit  with  an  OR  gate  of  three  inputs:  A,  B,  and  C,  and  an 
AND  gate  of  four  inputs:  one  is  the  output  of  the  above  OR  gate,  and 
the  others  are  D,  E,  and  F. 

Some  modifications  can  be  made  so  that  we  can  represent  non- 
tree type  circuits,  such  as  circuits  with  feedback,  reconvergent  fanout, 
etc.  Also  functional  elements  can  be  easily  represented  in  the  modified 
postfix  notation. 

In  order  to  accomplish  further  improvements  over  the  coding 
efficiency,  some  element  type  hiding  techniques  can  be  used.  Thus  we 
include  a declaration  paragraph  before  the  postfix  expression.  For 
instance,  the  above  example  can  be  rewritten  to  be 

BEGINDECLARE  3+  : Gl;  4.  : G2;ENDDECLARE 

r ABC  +G1  D E F #G2 


6.  CONCLUSION 

The  models  and  methodologies  given  above  present  a practical 
approach  to  the  design  of  a digital  simulation  system  which  is  cost- 
effective,  accurate,  and  easy-to-use.  The  effort  is  directed  at  both 
gate-level  and  functional  simulation.  The  implementation  problem  has 
also  been  considered.  These  have  definitely  pushed  the  state  of  the 
art  in  the  technology  of  digital  simulation. 
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C.  COMPUTER  SOFTWARE  SYSTEM  RESEARCH  IN 

THE  AREA  OF  MODULE  INTERFACE  t 
Professor  Ed  Thompson  and  Ali  Mohseni 


The  overall  objective  of  this  research  is  to  reduce  the  high  costs 
of  software  systems,  which  is  the  major  cost  of  practically  any  computer 
application  today.  In  order  to  achieve  this,  two  approaches  have  been 
selected . 


The  first  approach  in  this  research  is  the  development  of  pro- 
cedures and  languages  for  the  formal  specification  and  verification  of  the 
data  interfaces  between  program  modules.  The  research  on  data  interfaces 
between  program  modules  is  based  upon  the  development  of  a data  interface 
specification  language  and  its  Implementation. 

The  second  approach  is  focussed  on  modular  programming  in 
COBOL.  The  type  of  the  module  and  the  way  of  passing  parameters  has 
some  effect  on  costs  associated  with  execution  time  and  storage  spaces. 
For  a particular  problem,  good  questions  to  be  answered  are:  "What  type 

of  module?"  and  "How  are  parameters  passed?" 

A preliminary  module  interface  specification  language,  MIS-1, 
has  been  designed,  and  a precompliler  named  BRIDGES,  which  supports 
the  application  of  MIS-1,  has  been  implemented  and  debugged.  .A  MIS-1/ 
BRIDGES  users  manual  has  been  written,  and  experience  in  using  MIS-1 
and  BRIDGES  has  been  gained  in  a variety  of  contexts.  An  evaluation  of 
that  experience  has  been  completed.  Four  general  observations  can  be 
made: 


(1)  MIS-1  allows  powerful  validation  algorithms  to  be  executed 
by  BRIDGES,  but  makes  more  difficult  the  application  of  the 
technique  of  data  hiding. 

(2)  The  reaction  of  programmers  to  using  MIS-1  is  to  partition  a 
program  into  fewer  modules  and  to  use  single  modules  to 
implement  data  abstractions. 


t This  research  was  supported  entirely  by  the  Joint  Services 
Electronics  Program  under  Contract  F44620-76-C-0089. 
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(3)  BRIDGES  spends  most  of  its  time  doing  I/O  and  scanning 
Fortran  statements  rather  than  performing  its  validation  rou- 
tines . 

(4)  The  validation  algorithm  in  BRIDGES  is  heuristic  rather  than 
rigorous . 

Based  on  observations  1 and  2 above,  a new  specification  language,  MIS- 
2,  has  been  designed  to  be  more  concise  and  to  incorporate  data  hiding. 
Observation  3 implies  BRIDGES  could  be  integrated  into  a compiler  with 
minimum  overhead,  and  optimization  algorithms  could  be  used,  i.e.  , 
incorporate  compliler  flow  analysis  algorithms  as  part  of  its  validation 
technique. 

In  the  area  of  modular  programming  in  COBOL,  a series  of  small 
programs  were  designed  in  COBOL  to  answer  the  following  questions: 

(1)  Effect  of  type  on  parameter  costs 

(2)  Parameter  costs  vs.  the  number  of  parameters 

(3)  Implementation  costs  of  a module 

(4)  Data  reference  costs 

(5)  Method  of  parameter  passing 

(6)  Segmentation  and  its  limits 

The  results  from  their  compilation  and  execution  on  four  differ- 
ent computers,  the  Xerox  SIGMA-5,  the  CDC-6600,  the  DEC  PDP-10,  and 
the  IBM  360  are  being  analyzed. 

From  these  experiments,  we  will  obtain  some  information  which 
will  give  us  a better  understanding  of  the  trade-offs  Involved  in  using  sec- 
I tions,  as  opposed  to  subprograms,  so  that  intelligent  choices  can  be  made 

for  system  implementation.  Also,  this  information  will  be  useful  in  making 
recommendations  for  the  implementation  of  future  COBOL  compilers. 
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D.  FAULT  TOLERANCE  OF  CONTROL  CIRCUITS 

DESIGNED  USING  PETRI  NETS  t 
Professor  Tilak  Agerwala,  Professor  Ed  Thompson, 
and  Ernesto  Pacas-Skewes 


► 


I 


Petri  Nets  [5]  are  a useful  tool  for  the  design  and  description 
of  systems.  They  can  be  used  to  represent  concurrent  operation  and  ar- 
bitrariness, and  can  describe  systems  at  various  levels  of  detail.  They 
are  especially  suited  for  control  circuits.  Their  use  in  the  design  of  cir- 
cuits that  perform  algorithmic  procedures,  like  adders,  yields  circuits 
that  are  more  complex  than  the  ones  produced  by  other  methods.  Once  a 
Petri  net  for  a circuit  is  obtained,  it  can  be  directly  translated  into  an 
implementation  using  circuits  that  have  been  proposed  for  net  substructures. 

Three  dynamic  properties  of  Petri  Nets  have  direct  implications 
on  the  physical  implementation  of  the  net.  These  properties  are  safeness, 
liveness,  and  conflict.  The  implementation  of  a non-safe  net  will  produce 
a circuit  requiring  counters  to  keep  track  of  the  number  of  signals  arriving 
to  the  non-safe  places.  The  implementation  of  a live  net  is  guaranteed 
to  be  fully  utilized,  i.e,,  every  part  of  the  circuit  will  be  used.  The  im- 
plementation of  a conflict-free  net  will  be  a circuit  with  no  races. 

In  the  area  of  fault  tolerance,  research  has  been  done  towards 
the  modeling  of  fault  phenomena  in  control  circuits  by  Petri  Nets  [2], 
Structures  that  can  model  loss,  diversion,  active  deadlock,  skew,  mul- 
tiple instances  and  tardiness  of  control  have  been  developed.  In  this 
approach,  the  faults  that  can  occur  in  any  given  control  circuit  are  class- 
ified according  to  their  effects  on  the  circuit  behavior. 

The  objective  of  this  research  was  to  determine  the  fault  tolerance 
capabilities  of  circuits  designed  using  Petri  Nets.  In  order  to  accomplish 
this,  a circuit  implementation  of  an  element  used  to  implement  Petri  Net 
transitions  was  analyzed  for  fault  response.  The  effects  of  these  responses 
on  the  circuit  containing  the  faulty  element  were  then  considered. 


t This  research  was  supported  entirely  by  the  Joint  Services 
Electronics  Program  under  Contract  F44620-76-C-0089. 
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The  element  analyzed  is  the  c-circuit  proposed  by  David  Muller 
[4].  The  c-circuit  has  been  widely  used  in  the  circuit  implementation  of 
Petri  Nets.  This  element  has  been  mostly  used  in  circuits  utilizing  simple 
signalling  (where  tokens  are  represented  by  high  to  low  as  well  as  low  to 
high  signal  transitions)  , and  David  Misunas  f3]  has  suggested  using  the 
c-circuit  in  systems  utilizing  reset  signalling  (where  tokens  are  represented 
only  by  low  to  high  signal  transitions)  . 

The  c-circuit  shown  in  Fig.  1,  was  analyzed  for  stuck-at  type 
faults.  This  circuit  can  present  seventeen  distinguishable  faults.  The 
effects  of  six  of  these  faults  would  produce  unpredictable  operation  of  a 
control  circuit.  The  unpredictable  operation  results  from  the  fact  that  in 
the  presence  of  one  of  these  faults,  the  c-circuit  will  operate  properly 
for  some  set  of  input  and  output  values,  but  under  certain  input  and  out- 
put conditions,  the  output  of  the  circuit  will  follow  one  of  the  inputs.  For 
example,  if  the  second  input  to  gate  C is  stuck  at  zero,  the  circuit  oper- 
ates properlv'  for  a signal  value  of  zero  at  input  A.  For  a one  at  input  A 
though,  the  output  of  the  circuit  will  follow  input  B. 

The  other  eleven  faults  will  produce  an  easily  detectable  mal- 
function of  the  control  circuit.  This  is  because  under  the  presence  of 
any  one  of  these  faults,  the  c-circuit  will  undergo  at  most  one  transition. 
For  example,  If  input  A is  stuck  at  one,  the  output  of  the  circuit  will  at 
most  undergo  the  transition  from  low  to  high.  If  the  control  circuit  is  the 
implementation  of  a live,  safe  and  conflict-free  Petri  Net  that  operates  in 
a cyclic  form,  a fault  of  this  type  will  cause  the  operation  of  the  circuit 
to  stop  after,  at  most,  one  cycle.  Thus,  detection  of  this  type  of  fault 
is  easily  accomplished.  Fault  isolation  to  the  c-element  level  could  also 
be  done  if  there  are  means  of  determining  at  which  point  in  the  cycle  the 
circuit  stopped.  If  the  c-element  is  implemented  on  an  IC  chip,  then  all 
the  pin  faults  will  cause  the  operation  of  the  circuit  to  stop.  This  is  an 
attractive  feature,  pin  faults  being  the  most  common  faults  in  digital  cir- 
cuits . 


It  was  found  that  the  use  of  the  c-circuit  for  simple  signalling. 
Imposes  a restriction  on  the  structure  of  the  Petri  Nets  that  can  be  imple- 
mented. The  restriction  has  not  been  graph  theoretically  approached,  but 
the  Petri  Nets  shown  in  Fig,  2,  if  Implemented  using  c-circuits  and  sim- 
ple signalling,  will  not  operate  properly,  although  the  net  is  live,  safe 
and  conflict-free  for  the  marking  shown.  Tokens  (represented  by  signal 
transitions  that  are  recorded  by  both  c-clrcuits  implementing  the  transi- 
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tlons  sharing  an  Input  place)  are  not  actually  removed  wrhen  one  of  the 
transitions  fires.  In  Figure  2,  ^when  tj^  fires,  a signal  transition  will 
occur  at  the  shared  input  of  the  circuit  implementing  tg  . When  ta 
fires,  a signal  transition  is  received  at  the  other  input  of  tg  , so,  tg 
will  fire  without  waiting  for  tj^  to  fire  again. 

Structures  more  complicated  than  the  one  shown  in  Figure  2 
can  produce  the  same  malfunction.  Thus,  one  structural  restriction  is 
that  input  places  should  not  be  shared,  unless  the  Petri  Net  is  a Free 
Choice  Petri  Net  [4],  This  restriction,  however,  greatly  reduces  the 
synchronization  capabilities  of  Petri  Nets, 

For  a circuit  that  would  allow  place  sharing  in  simple  signalling, 
some  kind  of  interaction  between  circuits  with  a common  input  would  have 
to  exist,  and  some  kind  of  memory  elements  would  have  to  be  included  at 
the  inputs.  This  would  complicate  the  circuits  and  could  give  rise  to 
races  and  hazards.  Since  our  interest  is  in  fault  tolerance,  it  was  decided 
that  simple  signalling  should  not  be  used,  and  a circuit  for  reset  signalling 
was  designed.  The  circuit  shown  in  Figure  3,  implements  a transition  with 
two  input  places  and  two  output  places. 

The  circuit  in  Fig.  3 requires  more  logic  than  a c-circult  and 
thus,  can  have  more  faults.  The  fault  analysis  of  the  circuit  showed  that 
regarding  most  pin  faults,  the  circuit  has  the  same  fault  detection  charac- 
teristics of  the  c-clrcuit.  I.E.  , most  pin  faults  would  cause  the  system 
operation  to  stop.  Fourteen  out  of  the  sixteen  distinguishable  pin  faults 
will  exhibit  this  characteristic.  The  two  pin  faults  that  could  produce 
unpredictable  behavior  of  the  system  are  Rj  stuck-at-one  and  CL  stuck- 
at-one.  The  stuck-at-one  fault  in  Rj  will  produce  the  desired  effect  if 
a restriction  is  Imposed  on  the  Petri  Net  structure.  The  restriction  is  that 
no  transition  should  share  all  its  input  places.  This  requirement  is  not 
as  restrictive  as  the  one  believed  to  be  imposed  by  c-elements  in  simple 
signalling.  A stuck-at-one  in  CL  could  result  in  the  circuit  being  wrongly 
initialized,  producing  a system  malfunction.  For  a cyclic  control  circuit 
implementing  a live,  safe  and  conflict-free  Petri  Net,  this  fault  can  be 
detected  and  perhaps  even  marked  out.  If  such  a system  is  forced  to  under- 
go it's  cycle  once,  the  faulty  element  will  either  stop  the  cycle,  or  will 
be  properly  initialized  when  the  cycle  is  over. 
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Petfi  Nets  structures  that  will  model  each  of  the  faults  that 
can  be  present  in  the  circuit  (except  for  four) , have  also  been  developed 
as  part  of  the  research. 

Control  circuits  designed  using  Petri  Nets  require  more  hardware 
than  circuits  designed  by  other  techniques,  but  they  seem  to  have  more 
attractive  fault  tolerant  capabilities.  The  advantage  of  Petri  Nets  is  that 
the  fault  tolerance  features  of  the  designed  circuits  are  intrinsic  to  the 
design  method,  that  is,  the  designer  does  not  have  to  include  them  in  the 
design. 


Another  approach  that  has  been  taken  to  Implement  Petri  Nets 
uses  an  array  of  diodes  [1].  These  circuits  seem  to  be  more  suited  for 
place  sharing,  and  apparently  the  method  imposes  less  restrictions  on 
the  nets  that  can  be  implemented.  Research  on  the  fault  tolerant  capa- 
bilities of  these  circuits  will  be  done. 

Another  area  where  research  will  be  done  is  simulation  of  func- 
tional devices.  The  possibility  of  using  Petri  Nets  as  a description  tool 
in  the  simulation  of  functional  devices  will  be  studied. 
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E.  CONCURRENT  SYSTEM  SYNTHESIS  t 

Professor  Tiiak  Agerwala  and  Yong-chai  Choed-amphai 

Concurrent  systems,  characterized  by  multiple  loci  of  control, 
are  extremely  complex  and  difficult  to  analyze.  The  primary  objective  of 
this  research  is  to  formulate  procedures  for  the  systematic  synthesis  of 
systems  with  known  properties  so  that  no  subsequent  analysis  is  required, 
or  of  systems  whose  behavior  is  easily  verifiable.  The  synthesis  proce- 
dures should  allow  the  use  of  hierarchic  integration  (bottom-up)  tech- 
niques as  well  as  structural  decomposition  (top-down)  techniques. 

The  basic  approach  was  to  use  the  Petri  net  model  [1,2].  This 
model  has  proven  to  be  a convenient  tool  for  the  representation  and  study 
of  concurrent  systems.  Techniques  were  considered  for  synthesizing 
Petri  nets  with  known  properties. 

It  is  assumed  here  that  the  reader  is  familiar  with  the  Petri  net 
model.  Petri  nets,  though  known  to  be  good  models  for  concurrent  i 

systems,  are  difficult  to  analyze.  The  conventional  practice  is  to  design  ] 

a net  and  then  analyze  its  properties.  A number  of  iterations  of  design 
and  analysis  may  be  required  before  a "correct"  net  is  obtained.  This 
process  can  be  quite  expensive.  In  [3]  a Petri  net  is  analyzed  by  obtain- 
ing an  incidence  matrix  and  solving  a system  of  linear  equations  to  obtain 
invariants.  An  invariant  is  a minimum  set  of  places  such  that  the  sum  of 
tokens  in  the  places  in  the  set  remains  constant.  These  invariants  are 
then  used  to  prove  or  assist  in  proving  the  correctness  of  Petri  nets. 

In  contrast,  we  have  taken  a synthesis  approach.  A net  is 
constructed  from  a basic  building  block  according  to  a construction  rule. 

The  final  net  is  guaranteed  to  have  the  important  property  of  boundedness. 

(Boundedness  means  that  each  place  will  always  have  less  than  N tokens 
where  N is  an  a priori  known  value.)  This  property  is  essential  for  imple- 
mentation. In  addition,  all  invariants  of  the  net  are  obtained  as  by- 
products from  the  synthesis  process.  The  results  are  summarized  below. 

Proofs  of  theorems  have  been  omitted. 


+ This  research  was  supported  entirely  by  the  Joint  Services 
Electronics  Program  under  Contract  F44620- 76- C-0089 . 
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Definition  1:  Given  a net  = (P^  , T^,  A^)  (p^  is  the  set  of  places 

T is  the  set  of  transitions,  A^  is  the  set  of  arcs)  where  p , p^  are 
any  two  places  of  N^.  If  is  modified  to  N^  = (P^  , T^  , A^)  such 
that: 

r 1 

= T^  and 

A^  = :'_A^  - ’^(p,,t)  I i = 1,  2}  - ■i'(t,  p.)  ( i = 1,2}  ] 

Ul(P,t)  I (p.,  t)  e A , i = 1,2"^  U ^(t,  P)  I (t,Pi)  e A , 
'll  1 ^ 


then  p,  and  p,  are  said  to  have  been  merged  to  p to  form  a new 
net  Ng  , 

Synthesis  basics.  Let  Pj^  and  p^  be  any  two  places  of  a Petri  net  N^^  = 
(P^  , Tj^  , A^)  such  that: 

1)  each  of  them  is  in  some  invariant  of  N\, 

2)  for  all  t e T^  , (Pj  , t)  e A^  =»  (p^  , t)  / A^  and 

(t,p^)  e Aj  ^ (t,  Pg)  / Aj 

Let  111  , Iig  , . . . Ii„  and  Ig^  , Igg  , ...  be  invariants  such  that 

Pi  e Ii  j , 1 < j < n and  p^  ® Ig  1, . 1 < k < m . 

If  p^  and  p,  are  merged  to  Pg  , and  the  new  net  is  , then 

Theorem  1:  If  I is  an  invariant  of  N^  containing  Pi  and/or  p^  , 

then  I is  not  an  invariant  of  . 

Theorem  2:  The  set  of  all  invariants  of  not  containing  p^  is 

identical  to  the  set  of  invariants  of  not  containing  p^  and/or  Pg  . 


60 


^ A>.  ■ 


(III.  ELECTRONIC  COMPUTERS) 


Theorem  3:  All  invariants  of  Nj,  containing  Pg  are  obtained  as 
follows: 

For  every  pair  of  invariants  I^j  and  1^,  ^ such  that 

’ll  Li  = “ Li  ■ 

the  set  Ij  j U I^  j U f P3  1 “ [Pi  ' Pa]  an  invariant  of 
provided  that  no  other  invariant  of  N^  is  its  proper  subset. 

The  Synthesis  Procedure.  The  synthesis  process  is  as  follows: 

1.  Begin  with  a building  block. 

2.  Bring  in  another  building  block  or  merge  two  places.  The 
merging  can  be  performed  only  if  every  place  will  still  be 
contained  in  at  least  one  invariant  after  the  merging. 

3.  If  the  required  net  is  completed  then  stop;  else  go  to  step  2, 

Checking  for  the  merging  requirement  in  step  2 can  he  done  using 
theorem  1 to  3.  A building  block  is  any  Petri  net  all  of  whose  invariants 
are  known  and  every  place  is  in  some  invariant. 


Theorem  4:  The  Petri  nets  built  according  to  the  synthesis  proce- 

dure are  bounded. 

Theorem  5:  The  class  of  nets  synthesized  is  equivalent  to  the  class 

of  bounded  nets.  (The  definition  of  equivalence  is  taken  from  [1]  ). 

Example 

A very  simple  example  is  presented  to  illustrate  the  procedure. 
Two  "consumer"  processes  are  to  be  modeled.  The  consumers  should 
not  consume  at  the  same  time. 

The  consumer  processes  are  first  synthesized  as  shown  in 
Fig.  la.  t^  and  tt^  represent  the  process  of  consuming.  This  synthesis 
can  be  done  using  some  basic  building  blocks.  The  invariants  for  the 
net  in  Fig.  la  are: 


1 

i 

I 


i 
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'i  ' ■'sl 

12  " {P0'  P7l 

13  * ft.,) 

14  - {P2I 

and  Pg  are  then  merged  to  get  the  net  in  Fig.  lb.  The  invariants  are; 

I,  = {P4.  P5} 

h “ fP6'  P?) 

'3  * tPsJ 

Since  each  place  is  in  some  invariant,  the  net  is  bounded  for  any  initial 
marking.  Initially,  assume  a single  token  is  placed  in  each  of  the 
places  p,^  , Pg  and  P3  . Since  t^  and  t,^  both  have  an  input  place  in  the 
same  invariant  ({Pjjjf  and  since  the  sum  of  the  tokens  in  the  places  of 
this  invariant  is  1 , t^  and  t,^  fire  in  a mutually  exclusive  manner.  There- 
fore, the  two  consumers  cannot  consume  at  the  same  time.  A knowledge 
of  the  three  invariants  can  also  be  used  to  prove  the  absences  of  dead- 
lock. 

Plans  for  the  Next  Period 


1.  Further  rules  for  bottom-up  synthesis  will  be  considered.  The 
top-down  approach  will  also  be  considered,  i.e.  , how  can  transi- 
tions be  expanded  into  nets  so  that  certain  properties  are  preserved. 

2.  The  limitations  of  the  use  of  invariants  in  proving  the  correctness 
of  nets  will  be  investigated. 

3.  Using  the  formal  rules,  an  attempt  will  be  made  to  provide  guide- 
lines for  synthesizing 

• parallel  programs  (operating  systems) 

• data  flow  programs 

• software  for  distributed  processor  systems. 
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The  overall  aim  is  to  preserve  certain  properties  during  synthesis 
so  that  system  behavior  is  easily  verifiable. 
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COMMUNICATION  IN  CONCURRENT  SYSTEMS  + 
Professor  Tllak  Agerwala 


Well  developed  theories  of  computation  and  numberous  system 
analysis  and  design  techniques  have  allowed  the  construction  of  high 
performance,  essentially  serial  computing  devices . Speed,  efficiency, 
and  fault  tolerance  requirements,  however,  indicated  that  the  notion  of 
concurrency  or  parallelism  be  exploited  in  computing  systems.  Rapidly 
dropping  hardware  costs  made  it  feasible  to  build  systems  consisting  of 
interconnections  of  numerous  autonomous  computing  elements.  Such 
systems,  in  general  have  not  performed  even  close  to  original  expecta- 
tions. One  reason  is  insufficient  emphasis  on  the  important  role  of  com- 
munication in  such  systems  and  its  relationship  to  computation.  The  prime 
objective  of  this  research  is  to  provide  a framework  in  which  communica- 
tion aspects  can  be  formally  studied.  Another  objective  was  to  investigate 
techniques  for  the  synchronization  of  processes  in  a multiprogramm.ed  or 
multiprocessor  system. 

The  approach  to  studying  communication  aspects  was  basically 
to  "start  small."  A simplistic  model  of  comm  jnication  was  developed  and 
various  situations  were  analyzed  using  this  model,  to  provide  some  in- 
sight into  the  role  of  communication  in  concurrent  sys'ems  and  its  rela- 
tionship to  computation. 

A simple  model  of  communication  and  a measure  of  communica- 
tion complexity  was  proposed.  This  model  consists  of  a graph  where  each 
node  consists  of  a processor  with  some  local  memory.  The  processor  is 
capable  of  performing  simple  arithmetic/  logical  operations  according  to 
a program  stored  in  memory  or  wired  in.  The  arcs  of  the  graph  represent 
links  along  which  data,  signals,  etc.  are  transmuted  'rom  one  node  to 
another.  Everything  occurring  inside  the  node  is  termed  computation. 
Activity  on  the  arcs  is  termed  comma. dcation . The  total  number  of  prim- 
itive operations  (summed  over  all  node'  is  defined  as  "computation  com- 
plexity" and  the  total  number  of  arcs  Is  defined  as  "communication  com- 
plexity, " 


^Thls  research  was  supported  entirely  by  the  loint  Services 
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It  is  assumed  that  each  algorithm  statement  contains  sufficient 
information  to  deterraine  the  kind  of  system  graph  it  will  execute  fastest 
on.  The  complexity  of  the  graph  then  gives  the  complexity  of  the  algo- 
rithm. The  communication  provided  by  a particular  hardware  configuration 
can  be  determined  by  modeling  it  as  a graph. 

The  problems  listed  below  were  analyzed  using  this  simple 
model  for  communication: 

(1)  Finding  the  maximum  of  a sequence  of  N numbers 

(2)  Computing  the  Discrete  Fourier  Transform 

(3)  Permuting  data  in  a SIMD  machine 

(4)  Sorting  a sequence  of  N numbers 

(5)  Communicating  in  a computer  network 

(6)  Determining  the  minimum  spanning  trace  of  a network 

This  exercise  illustrated  the  following: 

(1)  Communication  plays  a critical  role  in  distributed  computing 
systems 

(2)  Systems  with  different  amounts  of  computation  and  commun- 
ication can  be  designed  for  performing  the  same  overall  task 

(31  Computation/communication  complexity  trade-offs  are  poss- 
ible 

(41  It  may  be  possible  to  increase  speed  of  execution  by  going 
to  higher  communication  complexity 

(51  Given  current  technology,  minimizing  execution  time  on  a 
sequential  machine  or  minimizing  hardware  may  not  be 
appropriate . 
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for  details,  the  reader  is  referred  to  [I].  The  work  on  imple- 
mentation was  primarily  a continuation  of  research  conducted  during 
1974-197  5 (supported  by  U.S.  Atomic  Energy  Commission  under  contract 
AT(ll-l-3288) ) . New  synchronization  primitives  were  proposed  and  shown 
to  be  complete:  they  can  represent  any  desired  interaction  between  pro- 
cesses without  the  use  of  conditionals.  The  usefulness  of  the  primitives 
was  Illustrated  and  disadvantages  of  existing  synchronizing  mechanisms 
pointed  out,  A straightforward  but  efficient  implementation  of  the  prim- 
itives was  also  obtained. 

For  details,  the  reader  is  referred  to  [2], 

The  following  will  be  undertaken; 

(1)  Development  of  a theoretical  foundation  for  studying  com- 
munication in  parallel  systems, 

(a)  The  model  obtained  thus  far  is  quite  simplistic.  It  does- 
not  adequately  cover  all  important  aspects,  (Ex,  pro- 
cessor memory  communication  within  a node  has  not  been 
considered)  . A more  sophisticated  model  will  be  devel- 
oped. 

(b)  The  measure  of  communication  complexity  is  quite  prim- 
itive, Other  measures  such  as  link  capacities,  number 
of  transfers,  fan-in/fan-out,  etc.  will  be  considered 
and  a more  realistic  model  developed, 

(2)  Application  of  theoretical  framework  to  the  design  and  anal- 
ysis of  concurrent  systems. 

(a)  Selected  architectures  and  algorithms  will  be  analyzed 
using  the  theoretical  framework. 

(b)  Existing  building  blocks  such  as  microprocessors  will 
be  evaluated  to  determine  the  complexity  of  communica- 
tion they  can  realistically  be  expected  to  handle.  Arch- 
itectural features  to  enhance  communication  will  be 
studied . 
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(c)  New  approaches  for  algorithm  analysis  will  be  considered. 
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G. 


SOFTWARE  VERIFICATION  t 
Professor  R.  T.  Yeh 


Our  primary  activity  during  the  past  twelve  months  was  con- 
centrated on  programming  methodology,  especially  in  the  area  of  program 
verification.  Two  distinct  parts  of  our  investigations  are  reported  here. 

(1)  Verification  of  non-deterministic  programs  by  predicate 

transformation.  In  this  worlc,  we  present  a methodology  for 
the  verification  of  both  deterministic  and  non-deterministic 
programs.  Our  method  is  based  on  the  concept  of  predicate 
transformer  introduced  by  Dijkstra.  The  predicate  transformer 
is  used  to  define  formal  semantics  of  programs  by  a mapping 
which  transforms  a set  of  states  after  the  execution  of  a 
program,  to  the  set  of  all  possible  states  before  the  execu- 
tion of  the  same  program.  Thus,  the  difference  between  the 
concept  of  determinism  and  non-determinism  loses  its  signif- 
icance in  this  semantic  concept.  The  focus  is  on  the  nature 
of  the  computation;  and  hence,  the  concept  of  iteration, 
not  how  the  program  iterates,  is  of  concern  in  the  current 
context . 

By  the  very  nature  of  the  semantic  definition  using  predicate 
transformers,  it  is  implicitly  assumed  in  our  method  that  the 
termination  property  of  a program  is  an  inherent  property  of 
the  algorithm  that  realizes  the  given  program.  Furthermore, 
termination  and  consistency  can  be  handled  by  the  same 
modus  operand!.  These  points  are  illustrated  by  a number 
of  examples. 

In  this  work,  the  main  results  were  the  method  for  generating 
the  weakest  preconditions.  We  have  Investigated  two 
methods;  by  generation  of  a seguence  of  approximations  of 
the  final  predicate  and,  by  solving  a solution  for  recursive 
equation.  These  two  techniques  indeed  can  be  used  to  comp- 


This  research  was  supported  in  part  by  the  Joint  services 
Electronics  Program  under  Contract  F44620-76-C-0089  and  in  part  by 
the  National  Science  Foundation. 
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liment  each  other  in  locating  the  loop  invariant.  The  gen- 
eration of  successive  turns  will  provide  a guideline  in  ob- 
taining a candidate, and  the  recursive  equation  provides  a 
test  for  checking  the  adequacy  of  the  candidate.  Therefore, 
these  two  techniques  are  useful  additions  to  any  verification 
system  based  on  inductive  insertion  method.  The  impact  of 
the  work  here  is  that  we  treat  total  correctness  instead  of 
partial  correctness.  On  the  other  hand,  from  the  practical 
viewpoint,  a basic  limitation  is  imposed  by  our  method  by 
the  exploding  complexity  for  finding  the  weakest  precondition. 

(2)  In  this  part,  our  investigation  is  centered  in  the  study  of  the 
nature  of  verification  processes.  We  have  considered  the 
inductive  mechanisms  for  five  techniques  of  verifying  inter- 
active/recursive program  structures:  inductive  assertion, 

predicate  transformers,  subgoal  induction,  computation  in- 
duction and  structure  induction.  We  have  discovered  that 
all  five  techniques  can  be  justified  by  a single  theorem  of 
inductive  proof  technique.  We  have  also  shown  that  all 
five  techniques  face  exactly  the  same  problem  of  finding 
properties  that  will  carry  an  induction.  Such  properties,  I 
refer  to  as  "Inductive  set".  We  have  proven  that  in  the 
inductive  set,  all  five  techniques  are  easily  related  to  one 
another  and  that  a program  proof  of  any  one  technique 
can  be  easily  converted  to  a proof  of  any  one  of  the  other 
techniques.  We  were  able  to  then  conclude  the  computer 
programs  simply  are  inductive  definitions  of  the  functions 
they  compute.  Induction,  therefore,  is  the  only  method  by 
which  they  can  be  proved  and  hence,  the  problem  of  induction 
is  therefore  unavoidable.  It  seems  that  the  same  method  is 
applicable  to  cahracterize  termination  proofs  and  we  are 
working  on  that  prospect. 
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A.  ELECTRONIC  CONTROL  SYSTEMS--RESEARCH  IN  ESTIMATION 
AND  ANALYSIS  OE  NONLINEAR  STOCHASTIC  SYSTEMS  t 
Professor  Steven  I.  Marcus 


The  objective  of  this  research  is  the  development  of  a theory 
for  the  design  and  analysis  of  nonlinear  stochastic  dynamical  systems. 
These  problems  are  of  considerable  importance  in  many  areas  of  commun- 
ications, control,  guidance,  and  tracking.  Optimal  estimators  previously 
have  been  derived  for  very  general  classes  of  nonlinear  systems.  Hov/- 
ever,  the  optimal  estimator  requires,  in  general,  an  infinite  dimensional 
computation  to  generate  the  conditional  mean  of  the  system  state  given 
the  past  noisy  observations.  This  computation  involves  either  the  solu- 
tion of  a stochastic  partial  differential  equation  for  the  conditional  density 
or  an  infinite  dimensional  system  of  coupled  ordinary  stochastic  differen- 
tial equations  for  the  moments  of  this  density.  Hence,  approximations 
must  be  made  for  practical  implementation. 

One  of  the  major  research  activities  undertaken  in  the  past 
year  under  the  partial  support  of  this  grant  is  an  investigation  of  the  use 
of  the  inherent  structure  of  particular  subclasses  of  nonlinear  stochastic 
systems  in  order  to  design  high-performance,  easily  implementable  sub- 
optimal  estimators  for  these  systems.  In  particular,  the  tools  of  harmonic 
analysis  on  Lie  groups  and  homogeneous  spaces  have  been  used  to  design 
suboptimal  estimators  for  nonlinear  systems  whose  state  evolves  on  a 
compact  Lie  group  or  homogeneous  space,  such  as  the  n-sphere  [2],  [4], 
[7]-[9].  These  models  have  applications  in  phase  tracking  and  demodu- 
lation, satellite  tracking,  and  rigid  body  orientation  estimation. 

The  approach  taken  in  this  research  consists  of  the  truncation 
of  the  infinite  set  of  coupled  equations  for  the  moments  of  the  conditional 
density  of  the  state,  by  means  of  an  assumed  density  approximation.  That 
is,  the  conditional  density  of  the  state  is  expanded  in  the  appropriate  har- 
monic expansion  (e.g,  , Fourier  series  on  the  circle  or  spherical  harmonics 


t This  research  was  supported  in  approximately  equal  parts  by 
the  Joint  Services  Electronics  Program,  National  Science  Foundation,  and 
Bureau  of  Engineering  Research  at  the  University  of  Texas  at  Austin.  Some 
of  the  work  reported  here  is  a continuation  of  research  begun  at  the  Mass- 
achusetts Institute  of  Technology  under  the  support  of  AFOSR,  NASA,  and 
NSF. 
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on  the  sphere)  . However,  the  higher  order  coefficients  are  not  merely 
set  equal  to  zero,  for  it  has  been  shown  that  this  does  not  result  in  good 
performance.  Instead,  it  is  assumed  that  the  conditional  density  is  a so- 
called  "normal"  density  on  the  particular  Lie  group  or  homogeneous  space; 
these  densities  have  many  of  the  same  properties  as  the  Gaussian  density 
on  1R“  , and  depend  on  a small  finite  number  of  parameters.  This  assump- 
tion effectively  truncates  the  infinite  dimensional  optimal  estimator  by 
expressing  higher-order  moments  in  terms  of  lower-order  ones,  without 
eliminating  the  higher-order  moments  (e.g.  , the  high-frequency  components 
in  a Fourier  series  expansion)  . This  approach  has  been  highly  successful, 
as  the  following  results  indicate. 

In  [2],  estimation  on  the  2-sphere  is  considered.  This  prob- 
lem is  motivated  by  the  desire  to  traci;  a satellite  in  a spherical  orbit  about 
a celestial  body,  where  it  is  assumed  that  the  orbit  is  perturbed  by  random 
effects  which  can  be  modeled  as  Gaussian  white  noise.  In  this  case,  a 
reasonable  model  for  the  evolution  of  the  position  x(t)  of  the  satellite  at 
time  t is  that  of  a stochastic  bilinear  differential  equation  driven  by 
Gaussian  white  noise;  the  observations  are  assumed  to  be  linear  obser- 
vations of  x(t)  corrupted  by  additive  Gaussian  white  noise.  The  optimal 
estimator  and  the  suboptimal  estimator  resulting  from  the  normal  assumed 
density  approximation  on  are  derived.  Monte  Carlo  computer  simula- 
tions, reported  in  [9],  show  that  this  suboptimal  estimator  performs  at 
least  as  well,  under  a variety  of  parametric  conditions,  as  the  extended 
Kalman  filter  and  the  Gustafson-Speyer  state-dependent  noise  filter,  and 
with  comparable  or  fewer  computations.  In  fact,  the  particular  estimator 
simulated  in  [9]  is  computationally  quite  simple,  having  been  truncated  . 
after  the  first-order  moments  of  the  conditional  density  (this  is  called  the 
first  order  filter)  . 

In  [4],  the  corresponding  results  for  bilinear  systems  evolving 
on  compact  Lie  groups  are  presented;  the  concepts  of  group  representation 
theory  are  utilized  in  the  analysis  and  in  the  derivation  of  the  normal 
assumed  density  approximation  on  a compact  Lie  group.  Both  optimal  and 
suboptimal  estimators  are  studied.  This  class  of  problems  is  motivated 
by  the  problem  of  estimating  the  orientation  of  a rotating  rigid  body;  the 
state  of  the  system  is  a direction  cosine  matrix,  which  evolves  on  the 
special  orthogonal  Lie  group  SO(3)  . 
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The  methods  discussed  here  are  successfully  applied  to  an 
optical  phase  tracking  problem  in  [7].  In  this  problem,  the  state  of  the 
system  is  the  phase  9{t)  , and  sin0(t)  is  used  to  modulate  the  intensity  of 
a coherent  light  source,  such  as  a laser,  which  is  pointed  at  a remote 
receiver.  As  opposed  to  the  previous  problems  in  which  the  observations 
consisted  of  nonlinear  functions  of  the  state  corrupted  by  additive  white 
noise,  the  observation  in  this  problem  consists  of  the  output  of  a photo- 
detector--a  doubly  stochastic  Poisson  process  with  rate  modulated  by 
sin9(t)  . The  objective  is  to  estimate  9(t)  given  the  past  observations  of 
the  photodetector.  The  state  (i.e.,  the  phase  9(t)  ) evolves  on  the  circle, 
so  Fourier  series  and  the  folded  normal  density  are  the  appropriate  tools 
for  designing  the  normal  assumed  density  suboptimal  estimator.  Monte 
Carlo  computer  simulations  have  shown  that  the  first  order  filter  designed 
by  this  technique  performs  better  at  low  frequencies  and  comparably  or 
slightly  worse  at  high  frequencies,  than  the  previously  used  "quasi- 
optimum" estimator  of  Snyder  for  this  problem.  The  phase  estimation 
problem  in  which  both  jump  process  observations  and  nonlinear  observa- 
tions with  additive  white  noise  are  available,  is  considered  in  [8]. 

A major  effort  has  been  undertaken  during  the  past  year  to  in- 
vestigate the  properties  of  nonlinear  stochastic  systems  driven  by  multi- 
plicative Poisson  white  noise,  with  the  eventual  goal  of  solving  estimation 
and  control  problems.  Such  systems  have  many  applications,  including  that 
of  modeling  a satellite  subject  to  random  micrometeorite  collisions.  In 
[1]  and  [5],  the  analysis  of  linear  systems  with  multiplicative  Poisson  white 
noise  is  considered.  The  stochastic  calculus  of  McShane  is  extended  to 
this  case  by  means  of  an  injection  procedure  which  is  similar  to  that  of 
McKean.  Other  properties  of  such  systems  are  studied,  and  moment 
equations  and  necessary  and  sufficient  conditions  for  stochastic  stability 
are  derived.  These  results  are  extended  considerably  in  [3]  and  [6]  to 
general  nonlinear  systems  with  Poisson  white  noise.  It  is  proved  that 
this  extension  possesses  the  same  desirable  properties  as  McShane's 
calculus.  This  work  represents  a considerable  step  in  the  direction  of 
providing  a unified  calculus  which  applies  to  a large  class  of  noise 
processes . 
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B.  RESEARCH  IN  LARGE  SCALE  SYSTEMS  AND 

DESIGN  TECHNIQUES  t 
Professor  Baxter  F.  Womack  and  Associates f 


Three  problems  have  been  investigated  in  the  area  of  large 
scale  systems  theory.  In  [1],  three  different  configurations  of  dynamic 
multiple-connected  subsystems  are  considered.  The  general  state  variable 
equations  for  circular,  polar,  and  circular-polar  configurations  of  n sub- 
systems are  given.  For  a specific  class  of  subsystems  and  n=  4,  com- 
plete input-output  relations  are  developed.  System  stability  and  control- 
ability  are  investigated.  In  [2],  properties  of  the  matrix  Lyapunov  equation 
are  developed  using  spectral  norms  to  establish  upper  and  lower  bounds 
of  eigenvalues  that  guarantee  a positive  definite  solution  of  the  Lyapunov 
equation.  This  work  is  continuing  and  will  be  useful  in  designing  large 
scale  systems  created  by  connected  subsystems.  In  [3],  a new  numerical 
calculation  method  is  treated  for  solving  simultaneous  ordinary  differential 
equations  of  higher  order  as  initial  value  problems.  The  concept  of  a diag- 
onalized modal  property  is  introduced  to  achieve  necessary  simplification. 

New  design  techniques  for  control  systems  are  being  investi- 
gated. In  [4],  it  is  demonstrated  that  a system  response  can  be  improved 
by  means  of  a multi-valued  switching  circuit  that  is  controlled  by  measuring 
on-line  system  variables  and  analogous  on-line  Nyquist  points.  Other 
research  in  progress  in  the  design  area  includes  (a)  nonlinear  synthesis 
procedures  in  the  frequency  domain  when  time  domain  tolerances  are 
specified,  and  (b)  nonrecursive  digital  filter  design  specifications,  arch- 
itecture, and  limitations. 

Sensitivity  theory  is  being  developed  and  utilized  in  several 
ways  in  continuing  research  projects.  In  one  project,  time  domain  sens- 
itivity functions  are  being  used  to  optimize  a large  dimensional,  non- 
linear dynamic  system.  In  two  other  projects,  new  and  more  useful  form- 
ulations of  sensitivity  functions  are  being  sought  which  can  be  used  on 
(a)  systems  with  multiple  eigenvalues,  and  (b)  adaptive  systems  with 
multiple  inputs  and  outputs. 


t This  research  was  supported  in  part  by  the  Joint  Services 
Electronics  Program  under  Contract  F44620-76-C-0089 , the  International 
Research  and  Exchanges  Board,  the  Japanese  Ministry  of  Education,  and 
Japanese  Scientific  Research  Grants. 
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New  and  novel  methods  of  characterizing  and  solving  control 
problems  are  continually  sought.  One  effort  in  this  direction  is  indicated 
in  [5]  and  (6],  where  coordinate  transformations  of  parameterized  arcs 
have  enabled  the  invariative  structure  inherent  in  the  Hamilton  and  La- 
grange equations  to  be  exposed  by  means  of  tensor  analysis.  A continuing 
goal  of  this  research  is  to  utilize  tensors  and  extensors  to  characterize 
problems  having  multiple  performance  indices.  Another  effort  focuses  on 
understanding  human  tasks  as  an  intermediate  step  in  formulation  of  new 
control  strategies  for  broader  classes  of  control  problems.  In  [7],  a model 
is  postulated  and  tested  to  predict  major  shifts  in  the  search  strategy  of 
human  subjects  analyzing  multimodal  surfaces.  In  [8],  a detailed  investi- 
gation of  human  search  experiments  reveals  that  the  human  concept  forma- 
tion of  the  multimodal  optimum  point  has  five  subconcepts,  namely  a point- 
segment  concept,  line-segment  concept,  plane-segment  concept,  local- 
surface  concept,  and  global-surface  concept. 

f Associates:  V.M.  Levykin,  J.J.  Montemayor,  SadaoAzuma, 

Y.  Tjandra,  D.O,  Nilsson,  C.R.  Brownina,  T.  Nagooka,  and  M . Oda. 
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C.  RECENT  RESULTS  IN  lARGE  SCALE  APPROXIMATIONS 

or  POINT  CONTROLLED  DISTRIBUTED  SYSTEMS  t 
Professor  R.  H.  Flake 


A new  numerical  method  for  optimal  point  control  of  elliptic  dis- 
tributed parameter  systems  has  been  developed  under  the  Joint  Services 
Electronics  Program  support  and  applied  in  several  studies  [1,6].  In  these 
studies,  a new  efficient  large  scale  system  numerical  solution  technique 
based  on  multiple  shooting  [3,  4,  9],  called  extended  multiple  shooting 
(EMS)  , was  developed  for  solving  mixed  boundary  value  problems  for  linear 
and  semi-linear  elliptic  systems  defined  on  IR^  and  irregular  regions  in 
]R^.  This  new  numerical  approach  may  be  employed  with  large  scale  finite 
difference  or  finite  element  approximations  for  the  operator.  The  method 
entails  an  order  reducing  feature  in  which  the  original  point  control  prob- 
lem is  transformed  into  an  equivalent  problem  possessing  much  fewer  un- 
knowns. This  feature  has  been  demonstrated  for  systems  involving  dis- 
cretization of  up  to  2,000  grid  points  which  were  reduced  to  equivalent 
problems  containing  around  fifty  variables  [6].  The  resultant  convergence 
rate  for  the  boundary  value  problem  solution  was  almost  three  times  as 
fast  as  the  successive  over-relaxation  technique  in  a large  estuarine 
modelling  problem  [8].  The  optimal  parametric  point  controls  were  obtained 
by  applying  the  Goldfarb's  search  algorithm.  These  earlier  applied  studies 
under  the  National  Science  Eoundation  sponsorship  involved  the  determin- 
ation of  optimal  waste  treatment  plant  capacities  for  Jamaica  Bay,  a small 
estuary  in  New  York  [5]  and  for  Corpus  Christ!  Bay,  a large  estuary  on  the 
> Gulf  Coast  of  Texas  [8],  Control  and  state  variable  inequality  constraints 

I were  employed  to  impose  realistic  waste  treatment  plant  characteristics 

' and  water  quality  standards  in  the  analyses.  It  was  apparent,  in  both 

studies,  that  minimum  waste  treatment  capacities  for  satisfying  water 
quality  standards  in  the  estuaries  could  be  achieved  at  about  one-third 
the  capital  and  operating  costs  required  by  the  uniform  plant  treatment 
! policies  usually  espoused  by  the  regulatory  agencies. 

A brief  problem  statement  for  optimal  point  control  of  elliptic 
systems  and  the  method  of  solution  are  outlined  below. 


This  research  was  supported  by  the  Joint  Services  Electronics 
Program  under  Contract  F44620- 76- C-0089 . Earlier  studies  applied  in 
this  topic  were  supported  by  the  National  Science  Foundation  under  Grant 
No.  GI-34870X. 

81 


(IV.  ELECTRONIC  CONTROLS) 


Optimal  Point  Control  for  Elliptic  Systems 

The  solutions  of  optimization  problems  have  been  obtained  for 
several  elliptic  distributed  parameter  systems  which  have  controls  applied 
pointwise  in  the  problem  domain.  A brief  problem  statement  follows: 

Find  a*eS  such  that 


N 

a 

inf  J(a)  - inf  ^ f(aj)  (1) 

aeS  a j = 1 

subject  to  state  equations,  which  for  a linear  second  order  elliptic  system 
are: 


(uC)  - (vC)  - K(x,y)  C 
6 X ' by 


h (a.)  6 (x  - X . , y - y ,) 


with  boundary  conditions 


(x,y)  e nc]R' 


{?.) 


YC(x,y)  + p, 


bC(x,y) 
b n 


0 


(x , y)  e b 0 
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where 

IyI  + Ul=  1 

and  state  variable  inequality  constraints 

C(x,y)  < g(x,y)  , (x,y)eQ  (3) 

where  the  set  S is  the  admissible  control  parameter  set  which  is  a 
N 

convex  subset  of  IR 

f (•  ) e ; g(  • , •),  h(  • , •)  € ^ E^(  • , , E^(  • , , 

u(  • , , v(  •,  •)  e 


6(x,y)  is  the  delta  function  at  x=0,  y = 0, 


E^  and  E^  > 0 , V(x,  y)  g q . 


The  above  problem  is  reformulated  in  terms  of  a per^v^bation  parameter 
vector  B = {a  ~ clq)  and  the  sensitivity  function  = S(x,y)  where 

h is  assumed  to  be  a linear  function  of  a . 
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Solution  Method 

The  optimization  problem  defined  by  (1)  , (2)  , and  (3)  , has  a 
unique  solution  and  is  equivalent  to  the  problem. 

Find  0 e /3  which  minimizes 


N 


a 


1(0)  - 


f(0 .) 


E . 

j = 1 


subject  to 


S (x.y)  0 < g(x,y)  - C(x,y) 
a 


(x,y)  efi 


where,  B e j?,  is  the  constraint  set  and  C (x,y)  is  the  solution  of  (2) 
with  a = ttg 

where  S (x,y)  is  the  sensitivity  function  of  the  state  to  the  control 
a 

parameters  a . 

This  sensitivity  function  satisfies  the  boundary  value  problem 


57  (n^) 


ds  MuSj 


S(vS) 


a 


ft  X 


ft  y 


- K (x , y)  S 


N 


a 


E 

i = 1 


ft  h . , . 

a(x  - X, , y - y.) 

6 a.  ] ] 


(x , y)  e n c IR 
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with  boundary  conditions 

yS  (x , y)  + p,  =0  , (x , y)  e ft  n . 

ftn 

The  proof  of  the  equivalence  of  this  alternative  formulation 
with  the  original  problem  is  given  in  [7],  This  result  means  that  the 
original  formulation  is  equivalent  to  a constrained  minimization  problem 
which  may  be  solved  for  discretized  problems  by  non-linear  programming 
methods.  This  approach  requires  both  the  solution  of  the  state  and  the 
state  sensitivity  equations. 

The  EMS  method  is  based  on  the  equivalence  between  the  solu- 
tion of  the  boundary  value  problem  and  an  initial  value  problem.  The 
equivalence  is  expressed  by  the  following  statement: 

The  unique  solution  to  the  boundary  value  problem 
(2)  is  equivalent  to  the  solution  of  a Cauchy  prob- 
lem defined  by  the  same  differential  operator,  with 
boundary  values  and  derivatives  on  a part  of  the 
boundary,  ftOp  , such  that  the  domain  of  dependence 
of  the  initial  values  contains  Q and  all  boundary 
conditions  are  satisfied  on 
Jifl,  = I (x,y)  e ftn/  ftQp]  . 

This  statement  may  be  verified  by  an  application  of  the  properties  of  the 
solutions  to  (1)  . The  extended  shooting  technique  is  a method  which  can 
be  used  to  calculate  the  unknown  initial  value  for  the  equivalent  Cauchy 
problem.  The  resulting  Cauchy  problem  is  not  well-posed,  in  general, 
as  can  be  verified  from  Hadamaard's  well  posedness  condition  [2],  In 


practical  problems,  this  condition  manifests  itself  by  divergence  of  the 
algorithms  [3]  due  to  numerical  ill-conditioning.  This  problem  is  allev- 
iated by  the  EMS  method  which  decomposes  the  original  domain  into 
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where 

^ 0 if  ^ n ^ 0 0 

and  is  based  on  the  following  result. 

The  unique  solution  to  the  boundary  value 
problem  (2)  is  equivalent  to  the  solution 
of  a multi-boundary  value  problem  defined 
on  Qj  , L j = 1, 2 , . , , ,N^^]  satisfying  the 

same  operator  of  (2)  , the  original  boundary 
conditions  and  additional  conditions 
C(x,y)  e wjk)  . 

The  proof  of  this  statement  follows  from  the  properties  of  the  solutions  to 
the  system  state  equations  and  from  the  previous  statement  on  equivalence. 
The  multiple  boundary  value  problems  defined  with  respect  to  the  subregions 
indicated  above  may  be  solved  by  shooting.  The  resulting  technique  is 
referred  to  as  the  EMS  method.  EMS  numerical  algorithms  have  been  em- 
ployed in  the  modelling,  analysis  and  optimal  control  of  elliptic  distrib- 
uted parameter  systems.  The  discrete  implementation  of  the  EMS  tech- 
nique involves  integrating  the  operator  equations  by  means  of  numerical 
integration  and  assuming  arbitrary  values  for  unspecified  initial  values 
at  the  subregion  boundaries.  A Newton  type  iteration  procedure  is  then 
employed  to  compute  a solution  which  satisfies  the  boundary  conditions 
at  the  final  points.  The  domain  of  dependence  of  each  initial  value  prob- 
lem was  selected  to  minimize  the  ill-conditioning  caused  by  growth  of  the 
trial  solutions.  Large  scale  approximations  to  the  elliptic  problems  of  the 
class  considered  have  been  solved  by  this  method. 

In  [5],  we  used  the  EMS  technique  to  solve  the  multi-point 
boundary  value  problems  satisfied  by  several  water  quality  variables  for 
Jamaica  Bay,  New  York,  A numerical  integrator  for  ordinary  differential 
equations  was  used  for  integrating  the  operators.  The  problem,  which 
had  four  point  controllers,  was  decomposed  into  eighteen  spatial  sub- 
regions.  The  multiple  shooting  algorithm  reduced  the  error  in  the  Newton 
iteration  for  the  correct  initial  values  to  10  in  two  iterations. 
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In  [6],  a multiple  shooting  technique  was  used  to  integrate  the 
dispersion-advection  equations  for  Corpus  Christi  Bay,  a large,  two- 
dimension,  estuary  in  Texas.  In  contrast  to  the  previous  problem,  a two- 
dimensional  integrator  was  derived  from  the  five  point  finite  difference 
approximations  [10]  for  the  operator.  Altogether,  269  mesh  points  were 
involved  in  the  problem  along  with  twenty-three  controllable  inputs.  The 
subregions  employed  in  the  shooting  algorithm  required  thirty-two  initial 
values  to  define  an  EMS  formulation  consistent  with  the  original  boundary 
value  problem.  The  Newton  iterations  for  satisfying  the  boundary  condi- 
tions converged  in  about  three  iterations.  The  method  was  found  to  be 
almost  three  times  as  fast  as  an  earlier  solution  technique  applied  to  this 
problem  that  was  equivalent  to  the  standard  successive  over-relaxation 
method . 


In  [1],  the  EMS  method  was  employed  to  solve  a discrete-time 
continuous  space  approximation  for  parabolic  distributed  systems.  High 
accuracy  spatial  approximations  were  achieved  with  algorithms  requiring 
less  than  twenty  initial  values.  These  studies  indicate  that  the  EMS 
method  has  excellent  convergence  characteristics  and  permits  the  compu- 
tation of  realistic  large  scale  approximations  for  elliptic  and  parabolic 
distributed  systems. 

Several  of  the  studies  mentioned  above  also  involved  the  compu- 
tation of  optimal  point  controls  for  large  scale  elliptic  problems.  In  each 
case,  the  sensitivity  equations  for  the  problem  relative  to  the  controls  were 
computed  by  the  EMS  method.  The  equivalent  formulation  then  reduces  to 
mathematical  programming,  which  is  defined  by  a non-linear  cost  function 
subject  to  linear  inequality  constraints  and  bounds. 

These  recent  results  demonstrate  that  the  EMS  approach  is  a 
powerful  method  for  solving  large  scale  approximations  to  elliptic  dis- 
tributed parameter  systems.  In  future  investigations  under  JSEP  sponsor- 
ship, we  propose  detailed,  theoretical  studies  of  solutions  for  parabolic 
problems  as  well  as  the  computation  of  open  loop  and  closed  loop  controls 
for  parabolic  problems  with  point  controls.  The  EMS  technique  should 
provide  a powerful  tool  for  solving  the  optimality  equations  for  parabolic 
systems.  It  is  also  important  to  note  that  the  numerical  state  decoupling 
and  order  reduction  characteristics  of  the  EMS  technique  is  of  a different 
nature  from  those  used  in  the  aggregation  and  other  contemporary  large 
scale  methods.  Various  current  techniques  of  decentralized  control  and 
other  suboptimal  policies  could  also  be  combined  with  the  resulting  high 
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resolution  EMS  models  to  study  very  large  problems. 
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A.  BISPECTRAL  ANALYSIS  OF  NONLINEAR  WAVE  INTERACTION  DATA  t 
Professor  Edward  J.  Powers  and  Y.  C.  Kim 


1.  INTRODUCTION 

In  the  past  we  have  developed,  partially  under  JSEP  sponsor- 
ship, techniques  which  have  proven  useful  in  experimentally  analyzing 
and  interpreting  fluctuation  data  associated  with  waves  and  instabilities 
[1,2].  The  approach  is  to  digitize  the  fluctuation  data  of  interest  and 
then  from  the  corresponding  raw  time  series  data  compute  various  spec- 
tral and  statistical  functions  of  interest.  Of  particular  interest  is  the 
fact  that  digital  spectral  analysis  based  on  EFT  techniques  is  capable 
of  resolving  a fluctuation  spectrum  into  a linear  superposition  of  the 
various  component  waves.  In  particular,  one  can  simultaneously  deter- 
mine the  amplitude,  frequency,  wavenumber,  and  coherence  of  each  wave 
present  in  a fluctuation  spectrum.  This  phase  of  the  work  is  now  sup- 
ported by  the  National  Aeronautics  and  Space  Administration. 

In  analyzing  fluctuation  data  characterized  by  nonlinear  or 
parametric  wave-wave  interactions,  "linear"  spectral  analysis  techniques 
are  of  obviously  limited  value.  Thus,  during  the  past  year  we  have  turned 
our  attention  to  the  problem  of  developing  procedures  that  are  suitable 
for  analyzing  "nonlinear  data.  " In  particular,  higher-order  spectral  tech- 
niques are  required  to  analyze  such  data.  For  example,  for  a system 
which  is  quadratically  nonlinear,  the  two-dimensional  frequency  spectrum, 
the  bispectrum,  is  the  appropriate  spectral  quantity  to  consider. 


2.  THE  BISPECTRUM 

The  bispectrum  is  best  introduced  by  way  of  application  to  a 
specific  problem.  Suppose  that  some  fluctuating  signal  has  been  digi- 
tized and  that  classical  "linear"  spectral  analysis  indicates  waves 
present  at  frequencies  , tUg  , and  ± ou^  with  wavenumbers  k^  , kg  , 
and  k^  ± ' respectively.  The  fact  that  waves  exist  at  sum  and  differ- 

ence frequencies  suggests  that  a nonlinear  interaction  has  taken  place. 
However,  if  all  the  waves  are  normal  modes  of  a system,  there  exists  the 


tThis  research  was  supported  principally  by  the  Joint  Services 
Electronics  Program  under  Contract  F44620- 76- C-0089 , and  in  part  by 
the  Texas  Atomic  Energy  Research  Foundation. 
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possibility  that  the  sum  and  difference  frequency  waves  were  spontane- 
ously excited.  Thus,  the  first  question  we  considered  was  how  does  one 
discriminate  between  spontaneously  excited  waves  and  those  which  are 
coupled  through  some  nonlinear  or  parametric  interaction?  The  second 
question  concerns  the  possibility  of  experimentally  determining  the  com- 
plex coupling  or  interaction  coefficient  from  the  experimental  data. 

We  consider  the  discrimination  problem  first  by  introducing  the 
bispectrum  [3]  and  reviewing  some  of  its  properties.  Mathematically, 
the  bispectrum  B(u)j  corresponds  to  two-dimensional  Fourier  trans- 
form of  a second-order  correlation  c(t^  ,Tg)  = E[x(t)x(t  + T^)x(t  + t^)] 
where  E denotes  an  expected  value.  The  fluctuating  quantity  x(t)  is 
assumed  to  be  statistically  stationary.  Since  the  second-order  correla- 
tion is  defined  in  terms  of  a third-order  product  moment,  the  bispectrum 
is  a mean  cube  spectral  density.  That  is,  B(uj^  , )dou^ dou^^  is  the  contri- 
bution to  the  mean  cube  E[x^  of  those  spectral  components  lying  in  the 
to  + dti)^  and  to  + dtu^  intervals  in  the  - plane.  Although 
the  bispectrum  is  equivalent  to  a two-dimensional  Fourier  transform  of 
c(Tj^  slso  be  written  as 

B(co^,0)2)  = yE  rX(aJ^)X(u)2)X*(cu)l  ; o)  = (1) 

where  X(<u)  is  the  Fourier  transform  of  x(t) , T is  the  time  duration  of  the 
x(t)  signal  and  the  asterisk  represents  a complex  conjugate  [4], 

Examination  of  Eq.  (1)  clearly  indicates  that  the  bispectrum  will 
be  zero  unless  the  following  two  conditions  are  met: 

1)  Waves  must  be  present  at  the  frequencies  , and 

2)  A phase  coherence,  or  phase  consistency,  must  be  present 
between  the  waves  mentioned  in  1). 

In  more  physical  terms,  we  note  that  if  waves  are  present  at  , and 

± u-'g  , then  condition  1)  is  satisfied.  Suppose,  however,  that  these 
waves  are  spontaneously  excited  normal  modes  of  the  system.  In  this 
case  each  wave  will  be  characterized  by  statistically  independent  random 
phases,  and,  thus,  when  the  statistical  averaging  indicated  by  the  ex- 
pectation operator  is  carried  out,  the  bispectrum  will  ideally  take  a zero 
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value.  On  the  other  hand,  if  the  sum  and  difference  frequency  waves 
were  generated  through  some  nonlinear  or  parametric  interaction,  then 
the  respective  phases  of  the  waves  would  no  longer  be  statistically  inde- 
pendent but  rather  a phase  coherence  would  exist.  As  a result  of  this 
phase  coherence,  the  statistical  averaging  will  not  lead  to  a zero  value 
of  the  bispectrum.  Thus,  our  approach  in  attacking  the  problem  of  dis- 
criminating between  coupled  waves  and  those  which  are  spontaneously 
excited  is  based  on  the  facts  that  coupled  waves  are  phase  coherent  and 
that  the  bispectrum  may  be  used  to  measure  the  degree  of  phase  coherence, 
Strictly  speaking,  a normalized  bispectrum,  known  as  the  bicoherence 
spectrum,  is  used  to  measure  the  degree  of  phase  coherence. 


3.  THE  BICOHERENCE  SPECTRUM 


The  bicoherence  spectrum  is  defined  as 


T%  [p(u)^)P(u)2)PH]  ^ 


UD  = tju  ± 

1 2 


(2) 


where  P(u^)  = ^ E jx(u))x*(u))J  is  the  auto-power  spectrum.  The  statistical 

properties  of  the  bicoherence  spectrum  may  be  interpreted  as  follows. 
Suppose  that  between  random  variables  X(cUg)  and  X(uj)  there  exists 


a relation 


X(u))  = A(cu^  ,cu2)X(uo^)X(uJ2) 


U)  = ± 


(3) 


where  physically  X(cu^)  and  X(tu^)  represent  the  amplitudes  of  the  primary 
interacting  waves  and  X(ou)  denotes  the  amplitude  of  the  resulting  sum  or 
difference  frequency  wave.  The  quantity  A(cu^  , is  the  coupling  coeffi- 
cient. The  goodness  of  the  quadratic  (i.e. , double  product)  relation  in 
Eq.  (3)  can  be  measured  by  a quadratic  correlation  coefficient  [5], 


p(oD^  ,tu„)  = 


E [x(tu^)X(u)2)X*(tu)] 


E [lX(a,^)X(co2)l^]^  E [lXH|^] 


U)  = . (4) 
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The  quantity  ,0)^)  basically  measures  the  degree  of  quadratic  correla- 
tion on  a spectral  basis  and  is  analogous  to  the  coherency  spectrum 
1^12  I which  is  a measure  of  the  degree  of  linear  correlation  between 
two  signals  on  a spectral  basis  [6].  Since  the  waves  at  and  are 
statistically  independent  of  each  other,  |p(<^i  to  Eq.  (2). 

With  the  aid  of  the  Schwarz'  inequality  one  can  show  that  | p(u)i  *'“2)  | * 
and  thus  b(tu^  ,cUg)  is  bounded  by  0 s b ^ 1.  Physically,  b(ou^  ,tu^)  will 
take  on  a value  close  to  unity  when  a nonlinear  (strictly  speaking, 
quadratic)  interaction  has  taken  place.  On  the  other  hand,  a value  of 
b(tu^,tu^)  near  zero  suggests  an  absence  of  quadratic  nonlinearity  and 
thus  suggests  that  any  waves  present  at  and  % ± are  probably 

spontaneously  excited  modes  of  the  system  rather  than  (quadratically) 
coupled  modes.  Application  to  experimental  data  is  discussed  in  a later 
section. 


4.  DETERMINATION  OF  COUPLING  COEFFICIENT 

In  a linear  approximation,  one  represents  a random  fluctuation 
in  terms  of  a linear  superposition  of  independent  waves.  In  a real  phys- 
ical system,  however,  the  spectral  components  of  a fluctuation  do  not 
always  appear  to  be  statistically  independent.  Rather,  some  spectral 
components  often  show  a high  degree  of  coherence  to  other  spectral  com- 
ponents which  can  be  measured  by  the  bicoherence  spectrum  as  discussed 
above.  The  high  degree  cl  coherence  between  spectral  components  sug- 
gests the  toUowing  "generalized  Fourier  series"  expansion  of  a random 
fluctuation,  i.e., 

x(t)  =2^x“’  (•»)■*“*  + X)  (».)x‘“(»  )e  ' ' ^*  + ...(5) 

where  A(uu^  ,^^)  is  a quadratic  coupling  coefficient.  After  Fourier  trans- 
forming each  side,  Eq.  (5)  becomes 

X(tu)  = X^^V)  + X^^^tu)  + ...  (6) 
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where 

. (7) 

(t'l  + uj^  = tu 

Clearly,  X^^\u))'s  represent  the  complex  amplitudes  of  the  statistically 
independent  "linear  waves"  and  is  the  amplitude  of  the  sum  and 

difference  frequency  waves  due  to  a quadratic  nonlinear  interaction  of 
X^^^{tu)'s.  As  a special  case,  suppose  that  there  exists  only  one  set  of 
and  tu^  which  satisfies  the  constraint  ou  = ± tD^  . In  other  words, 

only  two  waves  at  and  oj^  interact  and  generate  a third  and  a fourth 
wave  at  the  sum  and  difference  frequencies.  Then  Eq.  (7)  can  be 
written  as 

X{o))  - A(tu^,u32)X(uj^)X(iU2)  (8) 

which  is  equivalent  to  Eq.  (3).  Next  we  multiply  each  side  of  Eq.  (8)  by 
X*(cu^  )X*(ujg)  and  take  an  expectation  value  to  get 

E [x*(cu^)X*(u)2)X(cu)]  = A(ou^,«^2)E  |^|X(aj^)|^|X(uJ2)I^]  • 0) 

Making  use  of  the  fact  that  the  waves  at  and  at  are  statistically 
independent  of  each  other,  we  finally  obtain  the  following  expression 
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electromagnetic,  seismic,  and  acoustic  waves.  To  demonstrate  the 
utility  of  these  techniques  we  now  apply  them  to  wave  phenomena  observed 
in  plasmas.  Acoustic-type  waves  propagating  in  plasmas  give  rise  to  a 
space-time  modulation  of  the  plasma  density.  These  fluctuations  in  den- 
sity are  monitored  with  Langmuir  probes,  digitized  with  the  aid  of  Bioma- 
tion  8100  waveform  recorders,  and  transferred  to  magnetic  tape  under 
the  control  of  an  LSI- 1 1 microcomputer.  The  magnetic  tapes  are  then 
read  at  a later  time  and  the  data  analyzed  on  a CDC  6600  computer.  Al- 
though the  bispectrum  is  defined  over  the  entire  plane, 

we  make  use  of  its  symmetry  conditions  and  plot  it  only  in  the  first  quad- 
rant. Furthermore,  since  the  data  is  digitized,  the  highest  spectral  fre- 
quency plotted  is  given  by  the  Nyquist  frequency  f,^ . Consequently,  in 
the  figures  that  follow  the  bicoherence  spectrum  is  plotted  in  a triangular 
region  of  the  first  quadrant  defined  by  0 < fj^  s f^j  , fj  s f^  , f^  + f^  < f^  . 

Identification  of  Nonlinearly  Coupled  Waves 

In  Fig.  1 is  shown  a self-excited  fluctuation  spectrum  which  was 
observed  during  a study  of  the  evolution  of  drift-wave  turbulence  [7].  In 
the  middle  portion  of  the  figure  is  shown  a three-dimensional  plot  of  the 
bicoherence  spectrum,  and  at  the  bottom  is  shown  a contour  plot.  In 
particular,  we  wish  to  consider  the  waves  oscillating  at  frequencies 
f^  = 28,  fj^  = 53,  fj,  = 81,  f^  = 134,  and  f^  = 162  kFIz.  Each  wave 
was  identified  as  a drift  wave  propagating  mainly  in  the  azimuthal  direc- 
tion. The  phase  spectrum  (not  shown)  indicates  that  these  waves  are  not 
dispersive,  and  thus  the  nonlinear  wave  interaction  selection  rules  for 
wavenumber  will  be  automatically  satisfied,  if  the  frequency  selection 
rules  are  satisfied. 

Note  that  the  bicoherence  spectrum  In  Fig.  1 exhibits  several 
peaks  located  at  (f^  /f^),  (f^  ,f^j),  and  (f^  ,f^).  Let  us  first  consider  the 
peak  at  (f^,  ,f^).  Since  f^  = 2f^  , harmonic  generation  is  suggested  and 
in  fact  the  peak  at  (f^  ,f^)  indicates  that  the  wave  at  f^  is  the  genuine 
second  harmonic  of  the  wave  at  f^.  . Next,  the  peaks  at  (f^  ,f^)  and 
(f^^  ,f^)  also  indicate  phase  coherence  between  the  waves  f^  and  f , and 
the  difference  frequency  wave  f»  = fe  - h and  the  sum  frequency  wave 
• Thus  the  bicoherence  spectrum  indicates  that  the  waves 
present  at  f,^  / f[,  / f^  / f^  / f^  are  coupled,  rather  than  independently 
excited . 
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In  Fig.  2 is  shown  a turbulent  power  spectrum  and  below  It  the 
corresponding  bicoherence  plots.  In  the  turbulent  regime  we  might  synthe- 
size the  power  spectrum  in  terms  of  many  independent  oscillators  oper- 
ating at  different  frequencies.  These  oscillators  are,  of  course, 
characterized  by  statistically  independent  phases.  Under  such  condi- 
tions we  would  expect  the  bicoherence  spectrum  to  be  near  zero  and, 
indeed,  the  computer- generated  bicoherence  plot  is  very  near  zero  for  all 
combinations  of  f^  and  . 

Experimental  Determination  of  Coupling  Coefficients 

Equation  (10)  indicates  that  the  coupling  coefficient  A((u^ 
may  be  obtained  from  measurements  of  the  bispectrum  and  the  auto-power 
spectrum.  Our  work  [8]  on  this  phase  of  the  problem  is  fairly  recent  and, 
therefore,  the  results  are  somewhat  preliminary.  To  test  the  feasibility 
of  experimentally  determining  coupling  coefficients  we  have  applied  the 
approach  to  ion-acoustic  wave  harmonic  generation  in  plasmas.  The 
power  and  phase  spectra  are  shown  in  the  upper  half  of  Fig,  3,  Harmonics 
of  the  fundamental  (16  kHz)  are  clearly  evident  at  32  and  48  kHz.  The 
respective  azimuthal  mode  numbers  are  determined  from  the  phase  spec- 
trum as  m = 1,  2,  and  3,  Note  that  the  bicoherence  spectrum  shown  in 
the  lower  half  of  Fig.  3 exhibits  strong  peaks  at  (f;^  ,fj)  and  (f;^ 
corresponding  to  generation  of  the  first  and  second  harmonic. 

The  harmonic  generation  coupling  coefficient  A(f;^  ,f^)  relating 
the  amplitude  of  the  first  harmonic  (2f^)  to  the  fundamental  (f^)  was 
determined  in  the  following  way.  The  fluctuating  waveform  was  digi- 
tized and  Fourier  transformed  on  the  computer  using  the  FFT  algorithm. 
Next,  the  auto-power  spectrum  and  the  bispectrum  were  computed  from  the 
Fourier  transforms  and  the  coupling  coefficient  computed  using  Eg.  (10). 

The  results  are  given  in  Table  1 for  three  different  observations 
(l.e.,  observations  made  under  different  experimental  conditions)  of  ion 
acoustic  harmonic  generation.  In  Table  1 the  symbol  V is  used  to  indicate 
a dimensionless  coupling  coefficient  and  thus  distinguish  it  from  A(tu^  ,u^^) 
which,  as  indicated  by  Eq.  (8),  has  the  same  dimensions  as  1/X(u)).  We 
also  note  that  only  the  modulus  of  the  coupling  coefficient  is  shewn  in 
Tabla.  I . For  comparison  purposes  we  have  also  indicated  the  theoretical 
value  of  the  coupling  coefficient  based  on  a nonlinear  theory  which  was 
carried  cut  under  a research  grant  from  the  Texas  Atomic  Energy  Research 
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Foundation,  The  good  agreement  between  theory  and  experiment  gives  us 
confidence  that  indeed  coupling  coefficients  characterizing  quadratic 
nonlinear  interactions  may  be  determined  from  experimental  measurements 
of  the  bispectrum. 


6,  SUMMARY 

On  the  basis  of  the  examples  cited  in  the  preceding  section,  it 
does  appear  practical  to  utilize  digitally  implemented  bispectral  analysis 
as  an  aid  in  analyzing  and  interpreting  fluctuation  data  characterized  by 
nonlinear  or  parametric  wave  interactions.  It  should  also  be  stressed  that 
such  concepts  and  techniques  apply  to  all  wave  phenomena.  For  example, 
such  techniques  may  possibly  be  useful  in  analyzing  parametric- radar  data 
corresponding  to  scattering  from  a rough  sea,  Bispectral  analysis  may 
also  provide  additional  signature  information  by  which  various  radar  tar- 
gets may  be  classified. 


This  work  on  bispectral  analysis  of  nonlinear  wave  interactions 
will  be  described  in  greater  detail  in  journal  articles  that  will  be  prepared 
and  written  within  the  next  few  months. 
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Figure  1,  Computer  Generated  Plots  of  an  Auto- Power  Spectrum  (top) 
and  the  Corresponding  Bicoherence  Spectrum  (middle)  for 
a Case  of  Wave- Wave  Interaction.  Shown  at  the  Bottom  Is 
a Contour  Plot  of  the  Bicoherence  Spectrum. 
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Figure  2.  Computer  Generated  Plot  of  an  Auto- Power  Spectrum  (top)  and 
the  Corresponding  Bicoherence  Spectrum  (middle)  for  the  Case 
of  a Turbulent  Power  Spectrum.  Shown  at  the  Bottom  Is  a 
Contour  Plot  of  the  Bicoherence  Spectrum. 
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ACOUSTIC  WAVE  HARHONIC  GEftRATION 


BICOHERENCE  SPECTRUM 

12 


Figure  3.  An  Example  of  Ion-Acoustic  Harmonic  Generation,  Shown  in 
Top  Left  Is  the  Cross-Amplitude  Spectrum  and  on  the  Right 
the  Corresponding  Phase  Spectrum.  At  the  Bottom  Left  Is 
the  Corresponding  Bispectrum  and  to  the  Right  Is  Shown  the 
Corresponding  Contour  Plot. 
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B.  REFRACTIVE  INDICES  OF  , D^  , He , AND  Ne 

UP  TO  2000  ATMOSPHERESt 
Professor  A.  A.  Dougal,  S.  FI,  Baek  and  G.  Schiirger* 


The  refractive  indices  of  , D^  , He,  O^  and  Ne  were  mea- 
sured at  pressures  from  1 to  2060  atm,  25°  C and  6328a.  A formula  is 
presented  which  gives  the  refractive  index  over  a wide  pressure  range 
using  coefficients  derived  from  the  measurements. 

For  optical  experiments  in  super  pressure  gases  or  applications 
like  laser-fusion  of  gaseous  targets,  it  is  important  to  know  the  focal 
length  of  a lens  inside  a pressurized  gas  cell.  The  focal  length  for  a 
given  pressure  (P)  can  be  calculated  if  the  refractive  index  (n)  as  a func- 
tion of  the  pressure  is  known.  Earlier  workers  made  direct  measure- 
ments up  to  2000  atm  for  wavelengths  (X)  in  the  visible  range,  25° C, 
only  for  Ng  [1]  and  Ar  [2].  Results  are  presented  here  for  the  additional 
gases  of  H^  , Dg  , He,  O^  and  Ne.  When  the  electrical  polarizability 
(a)  and  the  particle  density  (N/V)  are  available,  the  refractive  index  is 
calculated  through  the  Lorentz- Lorenz  formula 

n^  - 1 = 4_  TT  a ii  . (1) 


Values  of  the  electrical  polarizability  are  given,  for  example,  in  [3]. 

Yet  these  values  are  derived  from  measurements  in  low  pressure  gases. 
The  results  of  [1]  and  [2]  show  a decrease  of  the  electrical  polarizability 
of  Ng  and  Ar  at  high  pressure.  This  decrease  is  around  1%  and  has 
little  influence  on  the  refractive  index  calculation  through  Eq.  (1).  The 
pressure  dependence  was  unknown  for  the  gases  investigated  here,  and 
could  be  larger.  Furthermore,  no  data  was  published  about  the  particle 
density  of  Og  above  135  atm  [4],  Therefore  direct  measurements  were 
necessary  to  get  reliable  values  of  the  refractive  indices  of  the  investi- 
gated gases  in  the  super  pressure  range. 


tThis  research  was  supported  in  part  by  the  Joint  Services 
Electronics  Program  under  Contract  F44620-76-C-0089,  and  in  part  by 
NSF  Grant  ENG  74-01785, 
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The  experimental  technique  utilized  an  interferometer  located 
inside  a super  pressure  cell  with  illumination  from  a He-Ne  laser  at 
632Sl  as  is  shown  in  Fig.  1.  Two  glass  discs  were  mounted  in  a holder. 
Three  screws  permitted  alignment  of  the  discs  relative  to  each  other. 

The  reflections  of  the  laser  beam  on  the  inner  surfaces  of  the  discs  pro- 
duced a circular  interference  pattern  at  the  photomultiplier  tube.  In 
order  that  the  reflections  from  the  outer  surfaces  of  the  discs  and  from 
the  cell  window  did  not  disturb  the  desired  interference  pattern,  each 
disc  had  a slanted  outside  surface.  The  two  inner  surfaces  were  tilted 
off  the  perpendicular  to  the  cell  window. 

For  each  run  the  cell  was  fully  pressurized,  and  subsequently 
gas  slowly  was  released.  A small  aperture  over  the  photomultiplier  was 
aligned  to  the  center  of  the  interference  pattern.  As  the  pressure 
decreased,  the  output  of  the  photomultiplier  undulated  and  was  recorded. 
An  m-fold  undulation  is  correlated  with  a refractive  index  change  An 
through  An  • 2d  = m\  , the  change  in  the  optical  pathlength  (d  is  the 
distance  between  the  glass  discs).  We  ignore  that  the  reflecting  surfaces 
were  slightly  less  than  perpendicular  to  the  laser  beam  axis.  The  purity 
of  the  gases,  as  stated  by  the  manufacturers,  was  99.999  % for  , He, 
and  Ar;  99.995%  for  Ne;  99.99%  for  Ng ; 99.95%  for  and  99.5%  for  Dg  . 

From  the  m-fold  undulations,  values  of  the  refractive  index  of 
the  different  gases  are  calculated  and  presented  in  Table  1,  and  graphi- 
cally displayed  in  Fig.  2.  The  gas  pressure  was  recorded  using  a bour- 
don gauge.  We  checked  and  recalibrated  the  gauge  by  measuring  the 
refractive  indices  of  N^  and  Ar  and  then  comparing  it  with  the  earlier 
data  of  Michels  and  coworkers  [1,2].  Our  measured  values  for  Ng  and 
Ar  are  also  included  in  Fig.  2.  The  maximal  error  of  the  pressure  reading 
was  ± 6 atm.  An  additional  uncertainty  in  the  calibrating  process  gives 
a total  possible  error  of  ± 12  atm.  The  error  in  n along  the  ordinate  in 
Fig.  2 is  less  than  0.05%.  The  change  in  the  disc  separation  d due  to 
the  hydrostatic  compression  of  the  disc  holder  was  calculated  from  the 
linearized  form  of  the  equation  for  volume  elasticity  Ad  = dAP/(3-B) 
where  B is  the  bulk  modulus.  An  elastic  compression  over  the  whole 
pressure  range  was  assumed. 

The  particle  density  can  be  calculated  from  the  virial  equation 

PV  = kN  T(1  + BP+.  . .)  (2) 

m A 
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1 

r 


[ 


where  Vj, , k,  N;  , T and  B ore  the  molar  volume,  the  Boltzmann  constant, 
Avogadro's  number,  the  absolute  temperature,  and  the  first  virial  coeffi- 
cient, respectively.  Eqs.  (1)  and  (2)  combined,  give 


(A*  + B*P  + 


2 

4TTP  ^ n +2 

3kT  2 

n -1 


(3) 


with  A*  = 1/a  and  B*  = B/a  , 

Using  the  measured  values  of  n for  various  P in  the  right  side 
of  Eq.  (3),  we  computed  by  the  method  of  least  squares  the  coefficients 
of  the  polynomial  in  P on  the  left  side.  A "best  fit"  is  obtained  for  all 
gases  for  a polynomial  A*  + B*P.  The  computed  coefficients  A*  and  B* 
are  listed  in  Table  2.  Our  results  for  N^  and  Ar  are  also  included  in 
the  table. 


From  Eq.  (3),  n can  be  expressed  as 


e 

Now  Eq.  (4)  permits  calculation  of  the  refractive  index  ) at  6328L 

and  25°  C over  the  full  pressure  range.  The  difference  between  n^^j  and 
n from  Table  1 is  shown  in  the  lower  portion  of  Table  2.  Also  shown  is 
the  difference  between  n^^j  and  , where  n^  ^^  are  the  values  taken 
from  the  literature.  They  are  either  directly  measured  or  calculated  by 
the  Lorentz- Lorenz  formula  from  the  measured  isotherms  together  with 
appropriate  values  of  the  electrical  polarizability 

-24  3 -24  3 

(a  = 0.790  X 10  cm  , Cf  = 0.775  x 10  cm  , 


-24  3 

O'.,  =0.21x10  cm  , 

He 


-24  3 

O'..  = 0,39  X 10  cm  ) [9,3] 

He 


The  relatively  good  agreement  for  these  gases  shows  that  the  change  of 
^ with  the  pressure  is  low.  Using  our  measurements  of  n and  the  particle 
densities  from  [7]  end  [8]  we  calculated  o'  for  He  and  Ne  by  Eq,  (1), 

The  results  gave  a decrease  of  o'  of  the  same  order  as  was  reported  pre- 
viously for  Ng  and  Ar, 
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The  absence  of  higher  order  terms  in  the  polynomial  used  in 
Eq.  (4)  and  the  results  for  n^  ^ ^ in  Table  2 give  us  confidence 

that  Eq.  (4)  can  be  used  for  still  higher  pressures  than  2000  atm. 


2. 

3. 

4. 


6. 


7. 


8. 


9. 
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Figure  1.  Diagram  of  Experimental  Arrangement  for  Refractive  Index 
Measurements  of  Compressed  Gases  Up  to  2060  atm. 
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Table  1.  Refractive  Index  n at  €3281  and  25 “C. 


p (atm) 

n (error  < 

0.05%) 

(error  ± 
12  atm) 

Ho 

°2 

He 

°2 

Ne 

1 ^ 

1.00014 

1 

1.00014 

1.00004 

1.00027 

1.00007 

50 

1.0070 

1.0070 

1.0018 

1.0142 

1.0037 

124 

1.0145 

1.0145 

1.0038 

1.0325 

1.0076 

185 

1.0210 

1.0209 

1.0056 

1.0491 

1.0108 

263 

1.0293 

1.0289 

1.0078 

1.0693 

1.0148 

331 

1.0351 

1.0347 

1.0093 

1.0834 

1.0179 

397 

1.0405 

1.0403 

1.0109 

1.0953 

1.0208 

466 

1.0458 

1.0456 

1.0124 

1. 1057 

1.0236 

535 

1.0506 

1.0503 

1.0138 

1. 1143 

1.0261 

602 

1.0551 

1.0548 

1.0151 

1.1216 

1.0285 

684 

1.0605 

1.0600 

1.0166 

1.1296 

1.0312 

756 

1.0647 

1.0642 

1.0179 

1.1357 

1.0335 

823 

1.0683 

1.0679 

1.0191 

1.1409 

1.0355 

888 

1.0718 

1.0714 

1.0201 

1. 1456 

1.0374 

951 

1.0752 

1.0748 

1.0212 

1.1496 

1.0392 

1018 

1.0784 

1.0780 

1.0222 

1.1539 

1.0410 

1098 

1.0823 

1.0819 

1.0235 

1. 1582 

1.0431 

1166 

1.0853 

1.0849 

1.0245 

1. 1617 

1.0447 

1236 

1.0881 

1.0878 

1.0255 

1.1649 

1.0463 

1301 

1.0906 

1.0905 

1.0264 

1. 1678 

1.0477 

1363 

1.0929 

1.0928 

1.0272 

1.1702 

1.0490 

1432 

1.0956 

1.0954 

1.0280 

1.1730 

1.0503 

1492 

1.0978 

1.0974 

1.0287 

1. 1753 

1.0515 

1567 

1.1004 

1.1002 

1.0296 

1.1781 

1.0529 

1636 

1.1027 

1.1024 

1.0304 

1.1803 

1.0542 

1699 

1.1047 

1. 1045 

1.0312 

1. 1825 

1.0553 

1763 

1.1068 

1.1065 

1.0319 

1.1845 

1.0564 

1827 

1.1088 

1. 1085 

1.0326 

1.1863 

1.0574 

1887 

1.1106 

1.1103 

1.0332 

1.1881 

1.0583 

1945 

1. 1123 

1.1118 

1.0338 

1. 1896 

1.0592 

2007 

1. 1142 

1.1137 

1.0344 

1.1913 

1.0602 

2060 

1.1158 

1.1153 

1.0350 

1. 1927 

1.0611 

f Data  frciti  [5] 
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C,  A GENERALIZED  MODEL  FOR  NONCOHERENT  OPTICAL 

CONVOLVERS  AND  CORRELATORS+ 

Professor  Jerome  Knopp 


In  previous  research  discussing  optical  convolution  or  correla- 
tion in  noncoherent  light,  three  basic  types  of  devices^  have  been 
described  for  performing  these  operations.  (See  Fig.  1.)  [1 , 2 , 3 , 4 , 5 , 6] 
In  each  of  these  devices,  the  two  dimensional  intensity  transmittance  of 
two  transparencies  is  convolved  or  correlated.  In  the  discussion  that 
follows,  reference  will  be  made  only  to  convolution  with  the  tacit  under- 
standing that  the  same  ideas  apply  to  correlation.  (A  device  which 
convolves  becomes  a correlator  by  simply  inverting  and  reversing  one  of 
the  transparencies.)  The  specific  operation  of  the  convolvers  shown  has 
been  explained  in  the  literature  cited  and  they  will  not  be  discussed 
further  except  in  a more  general  context  to  be  explained  below. 


In  past  research  experiments  using  a convolver  of  the  "misfocus" 
type  shown  in  Fig.  1(b),  it  was  found  that  moving  transparency  away 
from  the  lens  still  resulted  in  a convolution.  A geometric  analysis  of 
this  situation  showed  that  the  impulse  response  due  to  an  impulse  at 
results  in  a scaled  spatially  invariant  image  of  i no  matter  what  the 
distance  between  and  the  lens.  This  distance  effects  only  the  scal- 
ing of  the  impulse.  There  is  no  reason  to  mount  X>  as  close  as  possible 
to  the  lens  (unless  some  peculiar  scaling  between  and  x^  is  needed) . 
This  result  is  interesting,  and  applying  it  conceptually  leads  to  the 
generalized  version  of  a single  lens  convolver  shown  in  Fig.  2.  The 
resulting  optical  convolution  between  x^  and  x^  is  described  by 


■fThis  research  was  supported  in  part  by  the  Joint  Services 
Electronics  Program  under  Contract  F44620- 76- C-0089  and  in  part  by 
the  M.I.T.  Lincoln  Labs,  where  Dr.  Knopp  is  presently  on  the  staff. 

^ This  discussion  does  not  include  devices  which  use  mechan- 
ical shifters  to  shift  one  transparency  with  respect  to  the  other. 

^ This  analysis  holds,  assuming  the  usual  small  angle  assump- 
tions made  in  analyzing  optical  convolvers. 
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I(x‘ ,y')  = K 


//  - 


dx  dy  ^ 

mi 


where  K is  a constant  that  accounts  for  the  transilluminance  of  the 
incoherent  source.  The  constants  m^  and  m^  are  the  magnification 
factors  of  and  respectively.  These  factors  can  be  found  by  replac- 
ing one  transparency  with  a 6 - function  (mathematical  impulse  or  physi- 
cally a pinhole)  and  observing  the  image  size  of  the  other  transparency 
in  the  plane  of  convolution.  These  magnification  factors  can  also  be 
found  in  terms  of  the  object- image  geometry  given  in  Fig.  2,  It  can  be 
shown  that 


^2“°1 


m^  = h ~ °2  . (3) 

By  appropriately  adjusting  the  distances  shown  in  Fig.  2 and 
the  focal  length  of  the  lens,  it  is  possible  to  consider  any  one  of  the 
basic  three  types  of  convolvers  as  a special  case  of  the  generalized 
convolver  shown  in  Fig.  2.  This  conceptually  relates  the  convolvers 
shown  in  Fig.  1 and  in  practice  results  in  a more  flexible  system. 
Furthermore,  it  should  be  pointed  out  that  Eq.  (l)  can  be  extended  to 
any  multi-lens  system  where  a convolution  will  always  occur  at  planes 
of  misfocus. 

The  preceding  analysis  did  not  account  for  the  finite  size  of  the 
system  elements.  This  results  in  vignetting  and  bandwidth  limitations 
due  to  diffraction  effects.  These  effects  are  discussed  in  the  recent 
literature  [7] . 
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(Figure  Ic,  Convolution  Using  Shadow  Projection,) 
Figure  1,  Three  Basic  Types  of  Optical  Convolvers, 
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D.  INFRARED  NONLINEAR  OPTICS  IN  MOLECULAR  GASES  t 

Professor  Michael  F,  Becker,  Min  Ho  Kang,  and  Kang  Min  Chung 

We  have  demonstrated  optimized  third  harmonic  generation  in 
SFg  at  10.6  microns.  Phase  matching  and  optimum  focusing  were  stud- 
ied. Measurements  of  the  third  harmonic  susceptibility  show  that  there 
is  significant  population  redistribution  within  the  vibrational  energy  lad- 
der. The  susceptibility  v;as  saturated  at  higher  laser  intensities  due  to 
power  broadening  and  AC  Stark  shifting  of  the  resonant  energy  levels. 

These  effects  limit  the  third  harmonic  conversion  efficiency  to  less  than 
10“®  for  40  ns  CO^  TEA  laser  pulses. 

Optical  third  harmonic  generation  (THG)  has  been  successful 
[1,2,3]  in  the  visible  and  ultraviolet.  THG  in  molecular  gases  at  infrared 
wavelengths  had  not  been  demonstrated  [4,5,6]  until  late  1975.  This 
delay  is  due  to  the  smaller  vibrational-rotational  transition  matrix  ele- 
ments (as  compared  to  the  electronic  matrix  elements)  and  the  longer 
wavelength  itself.  However,  the  choice  [7]  of  a molecular  gas  such 
as  SFg  can  alleviate  the  first  obstacle  by  utilizing  a triple  resonance 
enhancement  in  the  TH  susceptibility  at  CO^  TEA  laser  frequencies.  We 
generated  a maximum  TH  power  of  10  mW  with  an  input  of  2 MW  on  the 
P(24)  line  of  the  10.6  micron  band.  Two  significant  factors  are  identi- 
fied which  limit  the  potential  conversion  of  triple  resonance  enhanced 
THG.  First,  absorption  at  the  fundamental  frequency  reduces  the  effec- 
tiveness of  phase  matching  to  an  enhancement  of  only  five  times  the 
non-phase  matched  conversion  and  prevents  any  improvement  by  further 
increases  in  the  SFg  pressure.  Second,  saturation  of  the  TH  suscepti- 
bility causes  the  obser\'ed  TH  power  to  fall  below  the  expected  cubic 
dependence  on  incident  laser  power.  i 

Absorption  measurements  on  SFg  at  the  P(20)  line  of  the  10.6 
micron  band  have  been  reported  for  power  densities  lower  than  2 MW/cm^ 
by  Armstrong  and  Gaddy.  [8]  We  extended  these  measurements  to  the 
range  of  40  MW/cm®  to  300  MW/cm®,  representative  of  our  third  harmon- 
ic generation  experiments.  The  saturated  absorption  coefficient  at  200 
MW/cm®  is  1/100  of  the  unsaturated  value  of  . 45cm~^  torr~^  and  agrees 


■•■This  research  was  supported  in  approximately  equal  parts  by 
the  Joint  Services  Electronics  Program  under  Contract  F44620-76-C-0089; 
the  National  Science  Foundation;  and  the  Texas  Atomic  Energy  Research 
Foundation. 
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with  the  theory.  [8,9]  This  residual  absorption  has  a deleterious  effect 
on  the  TH  conversion.  Figure  1 shows  the  calculated  behavior  of  TH 
power  as  a function  of  pressure  for  SFg  and  for  an  ideal  lossless  gas 
similar  to  SFg  in  ail  other  aspects.  The  points  represent  the  experimen- 
tal data  for  SFg  . It  is  evident  that  the  ultimate  conversion  in  a lossy 
medium  is  severely  limited  regardless  of  phase  matching  conditions. 

Figure  2 shows  the  calculated  TH  susceptibility  versus  frequency. 
Again,  the  points  represent  experimental  data.  Effects  of  saturation  of 
the  susceptibility  have  been  removed  by  extrapolating  the  cubic  slope  of 
the  TH  power  versus  input  power  data.  While  both  experiment  and  theory 
show  sharp  resonance  effects,  the  experimental  data  falls  below  the 
theoretical  curves.  This  is  due  to  the  severe  population  redistribution 
within  the  vibrational  quantum  system.  Our  theory  assumes  that  the 
molecules  remain  in  their  room  temperature  energy  distribution  and  that 
only  vibrational  ground  state  molecules  contribute  to  the  susceptibility, 
although  ^nis  will  not  be  true  for  power  densities  over  40  MW/cm^  . The 
strongly  saturated  absorption  is  clear  evidence  of  a population  redistribu- 
tion. In  fact,  at  the  absorption  line  center  (947.9  cm“^,  near  the  P(16) 
laser  line)  , the  deviation  of  the  experimental  susceptibility  is  strongest. 
Consider,  however,  the  effective  susceptibility  of  a molecule  in  the 
= 3 vibrational  level.  The  susceptibility  for  this  molecule  is  about 

10“^^  M*/V®  in  good  agreement  with  the  experiment.  A more  accurate 
comparison,  however,  is  very  difficult  since  we  must  know  the  exact 
population  distribution  within  the  complex  vibrational  energy  ladder  as 
a function  of  time  and  laser  intensity.  This  remains  as  future  work. 

Evidence  of  saturation  of  the  nonlinear  susceptibility  through 
power  broadening  [10]  and  AC  stark  shift  is  seen  in  plots  of  TH  intensity 
versus  fundamental  intensity.  Representative  curves  for  several  different 
laser  frequencies  are  shown  in  Fig.  3.  The  P(24)  line  gives  the  greatest 
ultimate  TH  signal  with  P(22),  P(26),  and  P(20)  being  next  highest,  in 
that  order.  If  one  resonant  transition  is  significantly  power  broadened, 
the  slope  is  reduced  from  3 to  2 , and  likewise  to  1 if  two  levels  become 
broadened,  and  to  0 if  all  three  levels  are  broadened.  Figure  3 suggests 
that  a single  level  is  power  broadened  for  the  P(22)  line,  and  that  two 
levels  are  power  broadened  for  the  P(20)  line  at  about  the  same  laser 
intensity.  No  power  broadening  is  observed  for  the  P(24)  line.  Quanti- 
tative analysis  of  the  power  broadening  is  in  progress  and  preliminary 
results  look  reasonable  if  AC  stark  shifts  are  Included.  It  appears 
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that  before  all  three  levels  become  broadened,  catastrophic  disruption 
of  the  SFg  occurs.  This  is  accompanied  by  broadband  visible  and  infra- 
red emission.  At  1 . 8 torr,  the  absolute  molecular  dissociation  probabil- 
ity [ll]  is  0.2%.  At  this  energy  density  (8  j/cm^)  a majority  of  the 
molecules  will  reside  in  the  continuum  of  upper  energy  levels.  [12]  This 
will  further  decrease  the  TH  susceptibility  and  the  generated  TH  power. 

In  conclusion,  we  found  that  the  TH  conversion  efficiency  for 
SFg  is  limited  to  less  than  10~®  by  the  effects  of  residual  saturated  absorp- 
tion, saturation  of  the  TH  susceptibility,  and  by  photo-dissociation  and 
breakdown.  The  second  of  these  effects  is  not  expected  to  be  strongly 
energy  dependent  while  the  last  two  are  energy  dependent.  Therefore, 
shorter  pulses  may  generate  higher  TH  powers  at  least  for  the  P(24)  line. 
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Figure  1 . 


Pressure  Dependence  of  the  Third  Harmonic  Power  With 
Different  Phase  Matching  Conditions. 

(a)  Perfect  phase  matching  for  a lossless  medium 
(Ak  = 0 , a = 0) . 

(b)  SF,  at  optimum  phase  matching  (Ak  = 0,2Ak,p  , a = a. . 

b SFg  Sh 


(c)  A non-phase  matched  lossless  medium  (Ak  = Ak,,  . a = 

^ 6 

(d)  Non-phase  matched  SFg.  Points  are  experimental  data , 
(A  k = A kj  , ct  = a ) . 

In  all  ca  ses  L = 2 3 . 6 cm , = 2 6 . 4 cm . In  (a)  - (c) , the 

focus  is  at  the  center  of  cell.  In  (d),  the  focus  is  L/8  from 
the  front  of  the  cell.  The  laser  energy  was  0.3  joule. 


t = 40  nsec,  on  the  P(20)  line. 


TH  susceptibility/molecule  (r"Vv^) 
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Figure  2.  Frequency  Dependence  of  Third  Harmonic  Susceptibility  in 

SFg  . Total  TH  Susceptibility  is  for  a Pressure  of  1.78  torr. 

(a)  Theoretical  susceptibility  [6]  with  a linewidth  of  .001  cm“^  , 
(Doppler  linewidth).  It  is  assumed  that  only  the  molecules 
in  the  ground  vibrational  state  contribute  to  the  suscepti- 
bility. 

(b)  Same  as  (a)  but  with  a Gaussian  linewidth  of  10.4  cm“’' 
(corresponding  to  the  entire  vibrational  absorption  band). 

(c)  Measured  susceptibility  for  P(30)  - P(8)  CO^  laser  lines. 
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TOTAL  TH  SUSCEPTIBILITY 


LASER  POWER  DENSITY  (GW/cm^  ) 


Figure  3.  Laser  Power  Density  Dependent  Third  Harmonic  Generation 
in  SFg  for  Various  Laser  Lines  in  the  10.6  Micron  Band. 
Configuration  is  the  Same  as  in  Fig.  1.  The  Cubic  Depen- 
dence of  the  Third  Harmonic  Was  Measured  to  Continue  Down 
to  5 uw;  This  Data  Is  Not  Shown  Here. 
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PICOSECOND  OPTICAL  MEASUREMENTS  OF 
ULTRAFAST  PHYSICAL  PHENOMENA  t 
Professor  M.  F.  Becker  and  R.  W.  Gunn 


In  the  experimental  results  reported  here,  we  have  measured 
the  tine  response  of  the  semiconducting-metallic  phase  transition  in 
VO  films  at  68°C  when  excited  by  picosecond  optical  pulses  at  a wave- 
lerrgth  of  I micron.  A single  pulse  from  a mode-locked  Nd:YAG  laser  is 
used  in  combination  with  a delay  line  apparatus  to  both  excite  the  tran- 
sition and  probe  it  at  a later  time.  In  this  way,  its  response  may  be 
measured  over  the  range  of  10  psec  to  10  nsec.  Details  of  the  laser  and 
delay  line  have  been  reported  previously  L 1 < ^ . 

Final  modifications  were  made  to  the  laser  system  irr  order  to 
secure  the  shortest  and  most  stable  mode-locked  pulses  possible  from 
this  type  of  laser.  In  order  to  eliminate  all  possible  etalon  effects  which 
tend  to  lengthen  the  mode-locked  pulses,  the  ends  of  the  NdrYAG  laser 
crystal  were  repolished  at  angles  of  3°  to  a prism-like  shape.  In  addi- 
tion, the  flowing  dye  cell  was  rebuilt  such  that  the  dye  region  was  di- 
rectly adjacent  to  the  laser  output  mirror  and  the  associated  optics  were 
suitably  wedged  so  as  to  reduce  etalon  effects  further.  In  this  config- 
uration, the  average  selected  pulse  length  was  3 5 psec  FWHM,  as  mea- 
sured by  an  autocorrelation  technique  [1,2]. 

Stability  of  the  laser  was  improved  through  precision  power 
supply  regulation,  dye  pump  vibration  isolation,  and  the  use  of  two 
intra-cavity  mode  central  apertures.  The  latter  modification  was  neces- 
sitated by  thermal  focusing  effects  in  the  Nd;YAG  laser  rod  at  repetition 
rates  of  20  to  50  Hz.  Thermal  instabilities  cause  beam  wandering  unless 
the  mode  axis  is  uniquely  defined  by  two  irises  separated  as  far  as  prac- 
tical. The  laser  now  operates  with  better  than  98%  of  the  pulses  well 
mode-locked.  This  is  close  to  the  statistical  limit  for  passively  mode- 
locked  systems . 

0 

Preliminary  measurements  were  made  on  a 2 500A  polycrystal- 
line VOg  film  deposited  on  a sapphire  substrate.  The  film  showed  a 
typical  change  in  resistivity  of  two  orders  of  magnitude  at  the  transition. 

t This  research  was  supported  entirely  by  the  Joint  Services 
Electronics  Program  under  Contract  F44620-76-C-0089 . 
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The  sample  was  mounted  in  the  delay  line  at  the  intersection  of  the  ini- 
tiating and  probing  beams.  The  transmission  of  the  films  was  monitored 
at  1.06  yi;  the  reflectivity  showed  a much  smaller  change  at  this  wave- 
length. Eligh  quality  areas  of  the  sample  showed  a 78%  decrease  in 
transmission  as  it  was  passed  through  the  transition  temperature  by  ther- 
mal excitation.  In  the  laser  initiated  dynamic  measurements,  the  sub- 
strate temperature  was  held  at  20°C  or  53°C;  in  either  case  the  results 
were  identical  since  the  film  heat  capacity  is  negligible,  compared  to  the 
latent  heat  of  the  transition.  For  laser  energy  deposition  densities  below 
the  latent  heat  of  the  transition,  no  change  in  transmission  or  reflection 
was  observed.  For  energy  densities  higher  than  the  latent  heat,  the 
transmission  was  observed  to  decrease  sharply.  However,  the  trans- 
mission decreased  by  11%,  or  only  15%  of  the  change  for  a thermally 
excited  transition.  Also,  the  decrease  occurred  in  less  than  10  psec,[3] 
the  resolution  time  of  this  measurement  technique. 

At  energy  densities  below  the  latent  heat  of  the  transition,  no 
photoexcited  carrier  absorption  was  observed.  It  was  expected  that  the 
photoexcited  carrier  lifetime  could  be  measured.  Either  the  carrier  life- 
time is  so  short  that  a significant  population  never  exists,  or  the  band 
structure  is  such  that  no  significant  change  in  the  transmission  occurs 
regardless  of  the  excited  carrier  population.  This  matter  will  be  studied 
further  in  order  to  estimate  the  photoexcited  carrier  lifetime.  Specific 
conclusions  about  the  mechanisms  driving  the  phase  transition  cannot 
be  drawn  until  the  fate  of  the  photoexcited  electrons  is  known. 

The  incompleteness  of  the  change  in  transmission  has  been 
observed  for  transitions  initiated  with  2 5 to  500  nsec,  1 micron  pulses 
also  [4],  This  incompleteness  may  be  interpreted  as  an  inhomogeneity 
in  the  volume  which  has  changed  phase,  or  as  the  first,  or  fast,  part  of 
a two  step  transition.  This  point  also  will  be  investigated  further.  New 
films  will  be  used  in  which  the  substrate  film  interaction  has  been  well 
characterized  L5]. 
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F.  CHEMICAL  KINETICS  OF  RARE  GAS-OXIDE  LASERS  + 

Professor  John  W,  Keto  and  Charles  F.  Hart 


1,  INTRODUCTION 

The  auroral  transition  in  atomic  oxygen  (o('So)  -♦  OC  Dg) ) ^t 
5577I  has  been  of  long-standing  interest  to  physicists  and  chemists. 
However,  recent  attention  has  been  focused  on  the  attractive  possibility 
of  its  use  in  a high  efficiency,  high  energy  laser  amplifier  [1].  This 
transition  satisfies  many  of  the  criteria  desirable  for  such  an  amplifier; 


(1)  The  lower  level  'D^j  is  more  rapidly  quenched  than  is  ‘So 

by  a large  variety  of  gases. 

(2)  The  transition  is  stable  against  superfluorescence. 

(3)  'So  has  a relatively  long  lifetime. 

(4)  Conventional  optics  may  be  easily  employed  with  the  5577i 

visible  wavelength. 

In  the  last  few  years,  a variety  of  rare  gas-oxide  mixtures  have 
been  shown  to  lase  [2,3],  yet  adequate  kinetic  models  of  the  specific 
reaction  sequences  necessary  for  a detailed  understanding  and  operation 
of  these  systems  is  lacking.  Insufficient  information  exists  about  the 
mechanisms  of  excitation  transfer  from  the  rare  gases  to  O(’So),  and 
accurate  values  for  many  rate  constants  are  required. 

This  report  describes  spectral  surveys  of  N^O  in  Ar,  in  Ar, 
and  Og  in  Ne,  which  have  been  taken  in  our  laboratory  as  the  preliminary 
step  in  a thorough  kinetic  study  which  will  be  started  shortly.  These 
results  have  Indicated  those  processes  where  an  in-depth  study  is  needed, 
and  have  pointed  to  Improvements  in  the  experiment,  which  are  now  being 
instituted . 


tThis  research  was  supported  entirely  by  the  Joint  Services 
Electronics  Program  under  Contract  F44620-76-C-0089. 
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2,  APPARATUS 

We  are  constructing  a sophisticated  apparatus  for  kinetic  studies 
on  electron- beam  excited  gas  samples.  The  computer- controlled  diagnos- 
tic system  will  monitor  both  time  dependent  fluorescence  and  pulsed  dye 
laser  absorption.  Over  the  past  year,  the  hardware  for  the  photon  count- 
ing and  timing  system  has  been  completed.  Final  programming  and  test- 
ing of  these  systems  for  kinetic  studies  is  near  completion. 

Most  previous  measurements  of  excitation  transfer  in  mixtures 
of  oxides  in  rare  gases  have  been  at  pressures  below  2 50  Torr.  However, 
experiments  on  rare  gases  suggest  greater  efficiencies  of  production  of 
excited  states  at  higher  pressures  [4]  and  we  have  therefore  planned 
kinetic  studies  to  pressures  of  30  atm.  At  high  pressures,  small  impurity 
levels  have  been  found  to  rapidly  quench  excited  states  [5,6].  Our  cur- 
rent cell  uses  standard  ultrahigh- vacuum  techniques.  By  baking  to  400°C 
at  10“^  Torr,  very  high  purities  are  achieved.  In  addition,  the  rare  gases 
used  are  cleaned  in  an  ultra  high- vacuum  titanium  furnace  with  a result- 
ant impurity  level  of  less  than  1 ppm.  With  the  closed  cell  required  for 
these  high  pressures,  difficulties  maintaining  the  dopant  gas  composition 
can  occur,  limiting  the  time  duration  of  studies  on  a single  sample.  We 
are  now  building  several  diagnostic  systems  for  monitoring  gas  composi- 
tion during  the  experiments . 

In  addition  to  the  previously  mentioned  construction,  preliminary 
measurements  on  oxide  mixtures  were  initiated.  For  these  measurements 
the  electron  accelerator  provided  a continuous  current  of  lyA  at  150  KeV. 
1/2  meter  visible  and  VUV  scanning  monochromators,  associated  optics 
and  cooled  photomultipliers  were  used  to  monitor  the  spectra.  In  the 
middle  of  the  past  year,  the  electronics  necessary  for  pulse  counting 
were  completed  and  a significant  improvement  in  signal-to-noise  ratio 
and  sensitivity  over  previous  results  using  phase- lock- loop  techniques 
were  achieved.  The  visible  system  was  absolutely  calibrated  using  an 
Optronic  Laboratories  tungsten  ribbon  lamp.  Spectral  surveys  of  N^O 
in  Ar,  O^  in  Ar  and  O^  in  Ne  were  completed. 
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3.  RESULTS  FOR  MIXTURES  OF  O^  IN  ARGON 

Shown  in  Fig.  1 is  the  spectra  of  the  O('So)  - OCD^)  transition 
at  5577A.  The  shape  of  the  line  core  was  found  to  be  independent  of  argon 
pressure  over  the  range  180  - 8500  Torr,  suggesting  predominately  free- 
free  rather  than  bound-free  transitions  [7],  At  very  high  sensitivities  a 
banded  structure  in  the  wings  of  the  O('So)  line  may  be  an  indication  of 
bound- free  transitions  which  have  been  previously  observed  in  Argon- 
oxide  mixtures  [8], 

Weak  bands  in  the  region  3700  - 4500 A were  noted  and  attributed 
to  the  Og  E -*  X^E,  ground  state  (Schumann- Runge)  transitions.  With 
increasing  Og  density,  there  also  appears  to  be  enhancement  in  intensity 
of  the  1900A  band  which  is  present  in  pure  argon.  The  specific  assign- 
ment of  this  band  in  pure  argon  is  uncertain  [9],  and  hence  it  is  difficult 
to  speculate  on  the  cause  for  enhancement.  This  increase  in  1900  A 
intensity  might  also  be  explained  by  an  as  yet  unidentified  transition. 

The  possibility  of  a second  transition  will  be  determined  by  the  future 
time  dependent  study. 

The  intensity  of  the  5577A  transition  as  a function  of  argon  and 
Og  densities  is  shown  in  Figs.  2 and  3,  and  that  of  the  1250A  band  due 
to  Arg*  transitions,  in  Figs.  4 and  5.  Ar^,*  is  known  to  be  formed  in  the 
three-body  reaction  [5,6], 

Ar*  + 2 Ar  - Ar^*  + Ar  , (1) 

and  we  have  verified  the  linear  dependence  of  the  Arg*  radiation  on  argon 
pressure  predicted  for  the  conditions  of  this  experiment. 

The  O(’So)  intensity  varies  at  low  pressures  as  the  argon  pres- 
sure to  the  3/2  power,  and  at  higher  pressures  as  the  pressure  squared. 

At  the  highest  pressures  of  argon  (11  atm)  it  appears  to  become  constant. 
Experiments  at  higher  pressures  will  be  done  to  verify  this  effect.  These 
results  contrast  with  those  of  previous  authors  who  claim  a linear  increase 
with  added  rare  gas  [10]; however,  these  have  been  in  experiments  at 
low  pressures  (<250T),  or  in  systems  where  the  impurity  levels  are  not 
as  carefully  controlled.  An  additional  explanation  may  be  the  affect  of 
dissociation  products  of  Og  on  the  formation  and  quenching  of  reactions 
leading  to  O('So)  production.  Since  this  is  a closed  cell  experiment, 
care  must  be  taken  to  include  the  affect  of  any  such  products.  A mass 
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spectrometer  and  laser  absorption  system  now  being  installed  will  moni- 
tor impurities  and  dissociation  products. 

The  dependence  of  O('So)  intensity  on  concentration  is  par- 
ticularly pronounced  and  indicates  the  reaction  mechanism  to  be  highly 
sensitive  to  this  parameter.  The  initial  constant  intensity  shown  in 
Fig.  3 is  probably  due  to  a small  {<  1 ppm)  residual  of  Og  in  the  argon. 

The  simplest  kinetic  model  one  can  postulate  to  explain  popula- 
tion of  the  0{'So)  level  is  that  of  collisional  dissociation  by  excitation 
transfer  from  bound,  excited  argon  molecules: 


O {‘So)  + O + 2Ar  . 


(2) 


This  reaction  is  probably  insufficient  to  explain  our  results  since  it 
implies  a linear  dependence  of  O('So)  on  Ar  concentration. 


It  is  of  special  interest  tc  determine  the  quantum  efficiency  for 
energy  transfer  from  Ar^  * to  O('So)  (if  this  is  the  major  transfer  mecha- 
nism to  O('So)) , since  this  wUl  determine  in  large  part,  the  overall  effi- 
ciency of  any  argon-Og  laser  amplifier.  Even  though  it  is  not  practical 
at  this  time  to  absolutely  calibrate  our  VUV  Optical  system,  by  making 
some  reasonable  estimates  of  various  VUV  optical  components  and  utili- 
zing the  absolute  calibration  of  the  visible  system,  we  were  able  to 
estimate  the  absolute  efficiency  of  the  VUV  at  1250^.  By  assuming  that 
reaction  (2)  is  the  major  quencher  of  Ar^  * we  have  determined  that  at 
150  PSIG  Ar  and  an  number  density  of  approximately  1 x 10^® /cm  , 
that  the  quantum  efficiency  for  excitation  transfer  from  Ar^  * to  O('So) 
is  40.3%  ± 17.7%.  Here  we  also  assume  de-excitation  of  O('So) 
proceeds  only  by  quenching  and  collisional  radiation.  Hence  our 
estimate  is  a lower  limit,  in  that  both  O3  and  ground  state  oxygen,  which 
are  possible  reaction  products  for  this  experiment,  have  larger  deactiva- 
tion rates  of  O('So)  than  does  [1]. 
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4.  MIXTURES  OF  N^O  IN  ARGON 

This  series  of  experiments  gave  results  comparable  to  that  of  Og 
in  Argon  indicating  that  N^O  very  rapidly  dissociates  into  N^  and  Og. 

Small  amounts  of  NO  are  also  produced  since  the  NO  Y bands  were  weakly 
present. 

We  attempted  to  measure  the  rate  of  N^O  dissociation  by  . 

monitoring  the  time  dependence  of  spectral  features  in  fresh  mixtures 
for  times  following  initial  excitation  by  the  electron  beam.  Since  intensi-  ; 

ties  did  not  vary  after  a few  seconds  from  initial  excitation,  dissociation 
to  an  equilibrium  configuration  is  rapid.  Mass  spectrometer  analysis  of 
gas  composition  is  necessary  for  a proper  understanding  of  the  kinetics  | 

of  these  mixtures  in  future  experiments.  I 

5.  MIXTURES  OF  Og  IN  NEON  | 

No  O('So)  radiation  was  observed  in  these  experiments,  indica-  l 

ting  that  the  energy  difference  between  Ne^*  and  O^  is  too  large  to  allow  | 

for  rapid  collisional  excitation  transfer.  Weak  Schumann- Runge  bands  ■ 

were  observed.  f 
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Figure  2.  O('Sg)  Absolute,  Integrated  Intensity  is  Shown  As  A Function 
of  Argon  Concentration  at  Fixed  Og  Concentration.  The  Error 
Bars  Represent  Both  Uncertainty  in  Absolute  Calibration  and 
in  Statistical  Deviation. 
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Figure  3. 


O('Sq)  Absolute,  Integrated  Intensity  As  A Function  of 


Concentration  at  150  PSIG  Argon.  The  Vertical  Error  Bars 
Represent  Uncertainty  in  Absolute  Calibration  and  Statistical 
Deviation  While  Horizontal  Error  Bars  Represent  Cumulative 
Error  Due  to  Repeated  Doping. 
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ELECTRON  BEAM  EXCITATION  STUDIES  OF 
POTENTIAL  DYE  VAPOR  PHASE  LASERS  t 
Professor!,  W.  Keto,  G.  Morowsky,  and  F.  K.  Tittel 


Currently  the  tunable  dye  laser  is  used  nearly  exclusively  for 
high  tunable  power  at  visible  and  near-ultraviolet  wavelengths  for  photo- 
chemistry, isotope-separation,  and  other  scientific  and  industrial  appli- 
cations. This  has  led  to  the  exploration  of  techniques  for  increasing  the 
power  output  and  efficiency  of  dye  lasers.  Currently,  dyes  are  dissolved 
in  a solvent  and  indirectly  pumped  by  flashlamps  or  lasers.  One  prospect 
for  increased  efficiency  and  power  is  direct  excitation  of  a dye  in  a vapor 
phase  with  an  electron  beam  or  discharge.  It  is  the  objective  of  this  work 
to  investigate  several  electron  beam  pumped  dye  vapor  systems  to  estab- 
lish optimum  conditions  for  electron  beam  excited  vapor-phase  dye  lasers. 

In  initial  work,  we  observed  for  the  first  time  strong  fluorescence 
from  an  electron  beam  excited  dye  vapor  [1].  These  experiments  used  a 
Pulsrad  110  electron-beam  accelerator  producing  a 0.7  5 Mev,  20kA,  20  nSec 
electron  beam.  The  experiments  studied  POPOP  vapor  and  established  the 
feasibility  of  using  various  buffer  gases  for  enhancement  of  the  energy 
transfer  to  the  dye  vapor.  Over  the  last  contract  period  an  improved  ver- 
sion of  the  laser  cell  was  designed  and  constructed.  This  cell  must  have 
optimum  access  for  the  electron  beam  and  laser  beam,  withstand  high 
buffer  gas  pressures  (up  to  10  atm) , and  be  able  to  withstand  the  tempera- 
tures of  450°C  necessary  for  dye  vaporization.  A diagram  of  the  laser 
cell  and  electron  beam  arrangement  is  shown  in  Fig.  1.  The  reaction  cell 
is  equipped  with  two  high  pressure  and  temperature  resistant  sapphire 
windows  and  gas  handling  connections  for  mixing  of  different  buffer  gases. 
An  oven  with  provision  for  optical  access  and  temperature  monitoring  can 
be  placed  over  the  dye  cell.  The  electron  beam  is  coupled  into  the  exci- 
tation region  transverse  to  the  optic  axis  through  a 2 5^m  thick  Havai^ 
foil  window  of  8 x 1 cm  cross  sectional  area.  A Rogowski  coil  surrounds 
the  graphite  cathode  for  monitoring  the  electron  beam  flux.  Technical 
details  of  the  cathode-anode  construction  are  show  in  Fig.  2,  From  the 

tThis  research  was  supported  in  part  by  the  Joint  Services 
Electronics  Program  under  Contract  F44620-76-C-0089 , the  Energy  Research 
and  Development  Administration , The  Robert  A.  Welch  Foundation,  and 
the  National  Science  Foundation, 
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burn  pattern  on  the  Havar  foil,  it  could  be  established  that  a knife-edge 
shaped  cathode  placed  12  mm  from  the  anode  window  produced  a very 
narrow  excitation  profile.  In  order  to  optimize  the  performance  of  the 
electron  beam  accelerator  it  was  necessary  to  establish  a compromise  of 
relatively  long  life-time  for  the  anode  foil,  effective  electron  beam  trans- 
mission through  the  anode  foil  and  optimum  impedance  matching  of  the 
diode  to  the  Blumlein  discharge  network  of  the  electron  beam  accelerator. 
With  an  anode-cathode  spacing  of  12  mm,  reliable  results  were  obtained 
for  at  least  20  shots  of  the  accelerator.  The  optical  emission  of  the  dye 
vapors  was  monitored  simultaneously  by  a grating  spectrograph  and  a 
fast-risetime  Valvo  UVHC  20  photodiode  whose  output  was  observed  on  a 
Tektronix  7904  oscilloscope.  In  order  to  differentiate  between  the  pump 
radiation  of  the  buffer  gas  and  the  dye  fluorescence  suitable  cut-off  filters 
were  used . 

We  have  studied  in  detail  the  excitation  process  of  three  selected 
dyes,  POPOP,  p-terphenyl,  and  the  xanthene  dye  N-92  by  monitoring  the 
electron  beam  excited  fluorescence  of  (i)  the  pure  dye  vapor,  (ii)  a dye 
vapor  with  a rare  gas  buffer  such  as  xenon,  and  (iii)  a dye  vapor  with  a 
buffer  gas  which  can  produce  u.v.  laser  emission  such  as  Kr  +-  NF 3 [2], 

Xe  + NF3  [3],  and  Ar  + N^  [4].  The  following  sections  describe  results 
for  each  dye . 


1.  ZANTHENE  DYE  N-92 

Zanthene  N-92  is  selected  as  a potential  lasing  dye  because 
its  absorption  band  perfectly  matches  the  emission  range  of  the  Ar  ^ Ng 
and  Xe  + NF3  systems.  Electron-beam  excitation  of  N-92  dye  vapor 
revealed  intense  fluorescence  whose  dependence  on  vapor  pressure  as 
derived  from  temperature  measurements  is  plotted  in  Fig.  3.  A slight 
saturation  is  evident  at  temperatures  approaching  380°  C.  At  this  temp- 
erature, breakdown  of  the  molecular  structure  seems  to  occur.  This  tran- 
sition is  irreversible,  since  the  fluorescence  does  not  reaapear  at  lower 
temperatures.  This  observation  is  in  agreement  with  a proposed  eJ^Dlana- 
tion  for  the  observed  discontinuity  in  the  vapor  pressure  curve  of  this  dye 
by  Schafer  and  Steyer  [5].  Furthermore,  Fig.  5 shows  that  the  N-92 
fluorescence  increases  more  than  10  fold  in  the  presence  of  Xe  buffer 
gas.  From  Fig.  5 we  may  also  conclude  that  collisional  energy  transfer 
from  Xe  to  the  dye  is  more  efficient  than  radiation  transfer  from  known 
rare  gas-halide  laser  mixtures,  such  as  XeF^  or  Ar-Ng**  . The  XeF  laser 
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mixture  was  prepared  according  to  Ref.  3 by  using  a mixture  of  Ar,  Xe , and 
NF3  in  the  ratio  of  250:25:1  at  a total  pressure  at  room  temperature  of 
1 atm.  The  presence  of  XeF  obviously  enhances  the  N-92  fluorescence 
to  a higher  degree  than  expected  from  the  Xe  content  of  the  mixture  as 
compared  with  the  addition  of  pure  Xe.  This  may  be  partially  due  to  the 
large  amount  of  Ar  present,  which  has  been  proven  to  be  nearly  as  effec- 
tive in  energy  transfer  as  Xe.  Compared  to  XeF,  Ar-N^  mixtures  are 
less  effective  in  generating  the  N-92  fluorescence  intensity  as  a function 
of  Xe  pressure.  As  shown  in  Fig.  4,  the  dye  intensity  increases  with 
increasing  xenon  pressure  reaching  a maximum  at  p - 2 atm.  There  are 
several  possible  kinetic  models  which  explain  this  behavior:  and  we  are 
currently  trying  to  determine  experimentally  which  one  is  most  probable. 


2.  POPOP 


The  most  promising  dye  in  the  vapor  phase  investigated  by  us 
so  far  is  the  dye  POPOP  [p-phenylene-bis-(5. -phenyl-2  . -Oxazole)].  It 
is  photochemically  stable  at  elevated  temperatures  which  is  a prerequisite 
for  long-term  laser  operation  and  its  fluorescence  peak  for  a dye  vapor- 
buffer  gas  mixture  is  centered  around  390  nm.  The  fluorescence  intensity, 
observed  within  the  spectral  range  of  30  nm  centered  at  360  nm,  is  plotted 
as  a function  of  vapor  pressure  for  pure  POPOP  vapor  and  various  buffer 
gases  in  Fig.  5,  apart  from  the  linear  increase  of  the  POPOP  fluorescence 
with  POPOP  vapor  pressure  up  to  30  torr.  It  is  interesting  to  note  that 
the  KrF  supported  dye  fluorescence  breaks  down  completely  at  higher 
vapor  pressure.  This  behavior  is  probably  a consequence  of  the  small 
amount  of  NF3  present  in  the  system  (He,  Kr,  and  NFa  were  mixed  in  a 
ratio  of  500:50:1  according  to  Ref.  2 with  a total  pressure  of  1 atm)  , which 
can  no  longer  compete  with  the  high  pressure  of  POPOP.  In  addition  to 
shielding  of  the  resonant  energy  transfer  in  the  rare  gas-halide  system  by 
the  presence  of  appreciable  amounts  of  dye  vapor,  the  addition  of  He 
finally  quenches  the  dye  fluorescence  as  shown  by  the  reduced  fluorescence 
intensity  of  the  system  POPOP  + He  (cf.  Fig.  8 and  Fig.  3 of  Ref.  1)  . 
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3.  P-TERPHENYL 

p-terphenyl  has  a fluorescence  range  which  covers  the  poten- 
tially useful  300  to  360  nm  range  of  the  UV  spectrum  and  is  convenient 
to  use  in  the  vapor  phase  due  to  its  high  vapor  pressure.  It  was  found 
(Fig.  6)  that  the  weak  fluorescence  of  pure  p-terphenyl  vapor,  compared 
to  the  dyes  N-92  ana  POPOP,  could  not  be  enhanced  to  the  same  degree 
of  brightness  by  the  addition  of  a KrF  mixture  and  other  buffer  gases. 

In  fact,  as  soon  as  the  dye  vapor  pressure  exceeded  40  torr  corresponding 

to  270°C  the  intense  KrF  radiation  was  quenched  completely,  as  observed 

previously  for  the  POPCPrAr-N^  mixture  [1]  and  for  POPOPrKrF  Even  the 
presence  of  pure  Xe  gas,  which  exhibits  the  most  effective  buffer  gas 
properties  in  the  case  of  the  two  other  vapor  phase  dyes,  led  to  saturation 
of  the  p-terphenyl  fluorescence,  although  the  de-excitation  energy  o 
Xe,  dimer  (7.2  eV)  matches  better  the  required  singlet  state  of  p-terphenyl 
compared  to  the  other  dyes  investigated.  This  might  be  a consequence  of 
an  increased  intersystem  crossing  rate  k,T  in  p-terphenyl  due  to  the  large 
heavy-atom  effect  of  Xe. 

The  successful  generation  of  intense  dye  vapor  fluorescence 
by  relativistic  electron  bombardment  via  a suitable  buffer  gas  demonstrates 
the  feasibility  of  electron  beam  pumping  of  complex  organic  dyes  ® 
vapor  phase.  From  the  spectral  profile  of  the  fluorescence  of  the  three  dyes 
studied  so  far  with  an  ultraviolet  absorption  range  from  220  to  ^90  nm,  e 
tails  of  the  energy  transfer  mechanism  for  the  xenon  buffer  gas  e in 

ferred.  A close  energy  match  of  the  dye  vapor  absorption  peak  to  the  de- 
excitation energy  of  the  X4  excimer(7.2eV)  ^oes  not  appear  to  be  ^ 
sential,  since  this  would,  in  fact,  suggest  that  p-terphenyl  would  fluoresce 
stronger  than  POPOP  (singlet-singlet  transition  energy  3.8  eV)  or  N-92 
(Singlet-Singlet  transition  energy  3.1  eV) . However,  the  strong 
intensity  obtained  for  electron  beam  excited  N-92  indicates  that  the  elec 
tron  stopping  power  of  xenon  (the  rare  gas  used  in  these  experiments  with 
the  highest  atomic  number)  makes  it  an  attractive  buffer  gas.  The  feas- 
ibility of  enhanced  excitation  of  the  dye  vapors  through  radiation  transfer 
from  rare  gas  halide  excimer  systems  was  examined  in  some  detail.  How- 
ever, quenching  of  the  excimer  fluorescence  by  the  dye  vapor  seems  to 
preclude  any  significant  amount  of  radiation  transfer. 
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Of  all  the  mixtures  studied,  only  N-92  and  Ar-Ng  mixtures 
seem  to  show  partial  radiative  energy  transfer.  This  may  be  due  to  the 
fact  that  the  buffer  gas  mixture  contains  a higher  concentration  of  the 
minor  component  (Ng)  as  compared  to  the  fraction  of  NF3  in  the  rare 
gas  mixtures. 

In  summary,  the  optimum  conditions  for  spontaneous  emission 
of  electron-beam  excited  dye  vapor-buffer  gas  mixtures  in  our  experimental 
arrangement  have  been  established.  The  presence  of  xenon  gas  at  a 
pressure  of  about  2 atm  appears  to  produce  the  best  fluorescence  enhance- 
ment. So  far,  internal  radiation  transfer  processes  have  not  been  effective 
in  the  excitation  of  vapor  phase  dye  lasers,  since  the  presence  of  even 
small  amounts  of  dye  vapor  inhibits  the  excitation  mechanisms  of  the 
efficient  rare  gas  halide  excimer  systems.  Work  is  now  in  progress  to 
establish  the  feasibility  of  electron  beam-pumped  vapor  phase  dye  lasers 
based  on  these  experimental  observations. 
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N-92  Fluorescence  Intensity  Versus  Vapor  Pressure  for  Pure 
N-92  Vapor  and  Several  Dye  Vapor-Buffer  Gas  Mixtures,  The 
Ar-N  Mixture  Ratios  Indicated  Refer  to  Pressure  Ratios  in  Torr, 
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Figure  4.  Saturation  Behavior  of  N-92  Fluorescence  In  the  Presence  of 
Xe  Buffer  Gas. 


Figure  5.  POPOP  Fluorescence  Intensity  Versus  Vapor  Pressure  for  Pure 
POPOP  Vapor  and  Several  Dye  Vapor- Buffer  Gas  Mixtures. 
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Figure  6.  Normalized  Absorption  and  Fluorescence  of  P-Terphenyl  in 

Cyclohexane  Solution®®  • with  Spectral  Position  of  KrF  Laser 
Line  Indicated. 
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H.  COLLISION-INDUCED  SCATTERING 

OF  LASER  LIGHT  t 

Professor  Lothar  Frommhold , Mike  Proffitt, 

Mark  Barrington,  and  T.  D.  Raymond 

The  major  effort  of  the  laser  light  scattering  work  is  directed 
toward  the  measurement  of  absolute  scattering  intensities  of  collisional 
pairs  of  atoms  of  the  heavier  rare  gases.  By  direct  comparison  of  our 
collision-induced  Raman  spectra  of  the  rare  gases  with  various  well- 
studied  spectra  of  nitrogen  [1]  and  hydrogen  [2] , we  were  able  to  obtain 
a number  of  calibration  points.  For  example,  at  a frequency  shift  of 
12  cm  and  at  a temperature  of  24°C,  each  pair  of  argon  atoms  scatters 
an  amount  of  laser  light  per  1 cm  ^ bandwidth  equivalent  to  2.58  x 10 
nitrogen  molecules  in  the  pure  rotational  line  J 2 -»  0 at  11.94  cm  , if 
irradiated  with  the  same  power.  This  information  can  readily  be  used  to 
obtain  scattering  cross  sections,  the  photoattenuation  coefficients  versus 
pressure,  or  all  similar  basic  light  scattering  data  desired  [1].  This  cali- 
bration in  argon  could  be  shown  to  be  valid  over  a gas  density  range  from 
essentially  zero  to  at  least  15  amagat.  We  believe  that  such  measurements 
have  not  been  reported  before.  Presently,  attempts  are  being  made  to  fur- 
ther improve  the  accuracy  of  these  numbers. 

We  feel  that  these  data  are  significant  for  the  following  reasons; 

(1)  In  this  way,  the  magnitude  of  the  anisotropy  of  the  polari- 
zability tensors  of  the  rare  gas  atomic  pairs  can  be  determined, 
which  are  not  sufficiently  well  known  at  present.  (This 
aspect  of  the  work  is  supported  by  the  Robert  A.  Welch 
Foundation) . 

(2)  Such  measurements  can  be  readily  used  to  predict  accurately 
the  attenuation  of  laser  beams  in  dense  gases,  as  far  as  it 
is  due  to  collision-induced  scattering  of  pairs  of  atoms. 


tThis  research  was  supported  in  part  by  the  Joint  Services 
Electronics  Program  under  Contract  F44620-76-C-0089 , and  in  part  by 
The  Robert  A.  Welch  Research  Foundation  for  Basic  Research  in  Chemistry. 
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It  was  pointed  out  by  Gerardo  [3]  that  for  every  light  scattering  process 
there  is  a corresponding  stimulated  scattering  process  which  might  well 
determine  an  upper  level  of  the  power  density  for  gas  lasers,  as  well  as 
for  the  propagation  of  high-powered  laser  pulses  through  dense  gases.  A 
thorough  understanding  of  the  collision-induced  light  scattering  processes 
of  atomic  pairs  hence  appears  to  be  of  interest. 

At  the  highest  gas  densities  (for  argon  above  about  12  amagats) 
our  spectra  showed  at  the  lowest  frequencies  (~3  cm  an  enhancement, 
which  was  not  expected  on  the  basis  of  the  theories  of  collision-induced 
spectra  of  binary  pairs.  We  were  now  able  to  show  that  this  new  feature 
of  the  spectrum  does  have  a simple  explanation:  it  is  produced  by  the 
far  wings  of  the  so-called  intra-collisional  process  predicted  by  von 
Kranendonk  [4]  some  time  ago.  In  other  words,  since  subsequent  collisions 
of  the  same  atom  produce  another  time-modulation  of  the  dipole  moment 
induced  by  the  laser  field,  there  is  a pressure-dependent  broadening 
mechanism  superimposed  with  the  collision-induced  scattering  of  isolated 
pairs.  This  intercollisional  scattering  gives  rise  to  a Lorentzian-like 
profile  of  a width  proportional  to  the  pressure.  The  far  wings  of  this 
Lorentzian  begin  to  show  up  in  our  experiments  at  high  pressures. 
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ELECTRON  SCATTERING  FROM  ATOMS  AND  MOLECULES^ 
Professor  Manfred  Eink  and  M.  Kelley 


During  the  last  year  research  was  carried  on  in  three  different 
areas:  (1)  determination  of  charge  density  distributions  in  atoms  and 
molecules  from  our  electron  diffraction  data;  (2)  the  thermodynamics  and 
composition  of  nonvolatile  compounds  like  inorganic  salts  and  metals; 
and  (3)  intramolecular  energy  transfer  among  various  vibrational  modes 
observed  by  superelastic  electron  scattering.  The  above  projects  are 
listed  in  older  of  their  degree  of  completion. 

With  our  electron  diffraction  unit  we  have  added  to  the  N^  data 
new  CO  and  O^  results.  The  density  maps  have  been  constructed  with 
the  technique  described  in  [1].  However,  this  approach  is  plagued  with 
several  approximations.  We  have  significantly  improved  the  evaluation, 
and  the  new  maps  are  presenting  distribution  functions  from  which  we 
expect  to  derive  several  moments  which  can  be  compared  with  optical 
measurements.  In  the  course  of  the  density  studies  it  became  apparent 
that  our  diffraction  unit  is  also  capable  of  recording  the  data  precise 
enough  that  the  molecular  structure  parameters  can  be  derived  with 
unprecedented  accuracy.  In  order  to  learn  the  limit  of  counting  machines 
as  a new  and  more  advanced  structure  tool,  we  engaged  in  a detailed 
study  on  the  data  correlation  as  it  applies  to  our  results.  The  results  are 
very  important  to  structural  chemistry,  as  well  as  for  the  structure  deter- 
minations of  amorphous  thin  films.  The  results  are  as  follows: 


r(C-F) 
i(C-F) 
r(F-  • -F) 
i(F*  --F) 


Table  1.  Structure  Analysis  of  CF^. 


Data  Without  Noise  Data  With  Noise  Previous  Results 
1.3325  ± 0.0003  1.3318  + 0.0006  1.323  + 0.005 

0.0437  ± 0.0015  0.0466  + 0.0024  0.049  + 0.005 

2.1830  + 0.0002  2. 1830  ± 0.0004  2.160  + 0.005 

0.0586  ± 0.0010  0.0595  + 0.0017  0.059  + 0.005 


Standard  0.0018  0.0031 

Deviation  of  Fit 

fThis  research  was  supported  in  part  by  the  Joint  Services 
Electronics  Program  under  Contract  F44620-76-C-0089  and  in  part  by  the 
Robert  A.  Welch  Foundation. 
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It  Is  obvious  from  Table  1 that  our  data  is  about  one  order  of  magnitude 
more  precise  than  previous  results.  However,  even  we  had  to  add  arti- 
ficial noise  to  cover  up  data  correlation.  Since  the  completion  of  these 
series,  we  have  improved  the  diffraction  unit  by  adding  an  optical  encoder 
to  record  the  angles.  We  hope  that  the  next  data  series  will  prove  that 
our  efforts  have  been  worthwhile,  and  a genuine,  uncorrelated  set  can 
be  recorded, 

A second  group  of  questions  has  been  answered  by  our  last  year's 
efforts.  Every  scattering  record  claiming  precision  in  the  .1%  range  has 
to  be  discussed  from  the  point  of  view  of  the  influence  of  the  finite  scat- 
tering geometry.  In  an  actual  experiment,  the  counted  events  originate 
from  an  extended  scattering  volume,  defined  by  the  volume  common  to  the 
electron  and  gas  beam,  and  the  detector  which  subtends  a definite 
acceptance  cone.  We  have  calculated  in  a computer  experiment  the 
extent  to  which  a model  cross  section  will  be  distorted  due  to  the  averag- 
ing. We  selected  the  range  of  parameters  as  they  apply  to  our  unit.  The 
results  however,  proved  to  be  of  such  general  nature,  that  requests  of 
our  peers  indicate  that  it  is  desirable  to  publish  the  method  and  the  con- 
clusions, and  to  make  the  numerical  results  available  through  a technical 
report. 


High  energy  electron  scattering  research  has  been  funded  for 
many  years  by  JSEP.  However,  our  results  have  found  such  enthusiastic 
response  that  we  feel  that  there  is  a very  good  chance  now  that  the  pro- 
gram might  be  picked  up  by  NSF.  To  this  end,  we  have  written  a two- 
year  $100,000  proposal.  If  the  funding  is  approved,  we  intend  to  switch 
the  JSEP  money  to  a new  experiment  which  we  are  presently  developing, 
and  which  is  outlined  below. 

Due  to  the  infrared  lasers  (COj  , HF  and  others)  an  old  topic 
became  very  interesting  again;  the  behavior  of  molecules  after  they  are 
excited  into  high-lying  vibrational  states.  The  question  of  mode  coupling, 
even  the  meaning  of  a mode,  is  no  longer  so  easy  to  determine  compared 
to  those  modes  lying  deep  in  the  potential  well  of  the  electronic  ground 
state.  It  is  our  intent  to  excite  several  molecules  with  various  lasers 
and  follow  the  internal  relaxation  mechanisms  with  low  energy  electrons. 
These  electrons  act  as  perturbers  and  cause  the  molecule  to  relax.  The 
excess  energy  is  transmitted  to  the  electron,  which  will  arrive  at  the 
detector  with  more  energy  than  when  it  started  from  the  electron 
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gun.  This  is  analogous  to  the  anti-Stokes  lines  in  Raman  scattering. 

The  phenomena  has  been  studied  for  an  electronic  transition  of  Na  by 
Hertel  and  Stoll  [2],  Similarly,  as  in  our  electron  diffraction  success, 
we  intend  to  develop  anew  technology  (the  superelastic  scattering  from 
vibrationally  hot  molecules),  and  thus  truly  contribute  to  the  deeper  under- 
standing of  the  energy  transfer  within  molecules.  We  might  be  too  late 
to  contribute  to  the  successful  completion  of  laser  isotope  separation  or 
laser  plasma  heating,  but  we  will  help  understand  why  certain  processes 
work  and  why  others  will  not.  This  is  not  just  an  esoteric  idea,  but  we 
have  worked  in  low-energy  electron  scattering  for  some  time,  and  we 
think  the  prospect  has  a great  enough  chance  of  success  that  the  challenge 
should  be  accepted  and  brought  to  a productive  stage. 

To  this  end  we  obtained  through  the  government  surplus  program, 
a $30,000  stainless  steel  vacuum  chamber  complete  with  the  10"  diffusion 
pump  and  30  cfm  rotary  pump.  We  have  built  a triple  set  of  Helmholtz 
coils  to  compensate  the  magnetic  field  of  the  earth  to  less  than  1 m Gauss. 
At  present  we  are  building  an  electron  gun  which  hopefully,  due  to  its  new 
electron  optical  system,  has  no  anomalous  energy  spread,  and  should 
therefore  deliver  a strong  monoenergetic  electron  beam  after  the  exit  slit 
of  the  spherical  monochromizer. 

The  third  project  we  are  working  on  is  funded  by  the  Robert  A. 
Welch  Foundation.  The  questions  we  try  to  answer  are  directed 
toward  the  molecular  structure  of  nonvolatile  compounds.  This  project 
was  funded  in  the  first  two  years  by  JSEP,  until  the  Robert  A.  Welch 
Foundation  supplied  continuing  funds.  We  have  at  present  focused 
on  rock  salts  and  are  continuing  to  measure, with  our  quadrupole  mass 
spectrometer,  the  reaction  constant  of  dimerization.  The  nozzle  tempera- 
ture is  constant  to  .5°C  and  the  mass  spectra  are  calibrated  through  a 
background  gas  mixture  of  Ne,  Ar,  Kr  and  Xe. 
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J.  LOW  ENERGY  MOLECULAR  BEAM  MACHINE  t 

Professor  W.  W.  Robertson,  Professor  L.  Frommhold  , 
and  Georges  Jamieson 


The  molecular  beam  machine  had  been  disassembled  to  permit 
accurate  alignment  of  stagnation  chambers  with  the  two  molecular  beam 
axis.  The  machine  has  now  been  reassembled. 

The  chopper  thiming  signal  electronics  have  been  modified  to 
increase  speed  and  improve  reliability.  The  light  bulb  has  been  replaced 
by  an  infrared  LED,  and  the  slow  photodiode  has  been  replaced  with  a 
much  faster  Motorola  MRD  510.  The  risetime  of  the  chopper  signal  is  now 
about  15  nsec,  a considerable  improvement  over  the  previous  2 y,  sec  rise- 
time. This  will  permit  accurate  time  of  flight  analysis  of  the  molecular 
beams  and  also  of  the  scattered  molecules.  The  molecular  beam  chopping 
disk  currently  consists  of  a 2"  stainless  steel  disk  with  ten  slots  cut  in 
it.  The  slots  are  about  equal  in  width  to  the  adjacent  solid  part  of  the 
disk,  which  means  that  time  of  flight  data  can  be  extracted  only  by  exam- 
ining the  leading  or  trailing  edge  of  the  beam  pulse  as  the  beam  is  chopped. 
It  is  more  desirable  to  have  a narrow  slot  which  produces  a short  beam  pulse 
so  that  risetime  measurements  do  not  have  to  be  made.  The  chopper  wheel 
will  be  changed  in  the  near  future  to  include  a slot  comparable  in  size  to 
the  molecular  beam  width.  The  current  chopper  disk  will  be  used  until  the 
preliminary  detector  evaluation  studies  have  been  finished. 

The  molecular  beam  detector  is  currently  being  evaluated.  The 
detector  consists  of  a Pierce  gun  ionizer  which  ionized  the  beam  particles, 
lenses  which  focus  the  ions,  and  a particle  multiplier  to  detect  the  ions. 

The  Pierce  gun  is  functioning  properly,  with  the  exception  that  the  1"  long 
filament  sags  when  it  is  heated.  The  filament  must  be  positioned  exactly 
correctly  or  it  will  touch  the  electron  gun  grid  making  the  detector  unuseable. 
When  the  present  filament  burns  out,  the  filament  holder  will  be  replaced 
with  one  that  allows  the  filament  to  be  spring  loaded  (any  failure  which 
requires  the  opening  of  the  ultra  high  vacuum  section  of  the  beam  machine 
usually  results  in  a delay  of  one  day  before  the  mai;hine  is  ready  to  produce 
a beam,  thus,  reliability  of  parts  is  extremely  important). 


t This  research  was  supported  entirely  by  the  loint  Services 
Electronics  Program  under  Contract  r44620-76-C-0089, 
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A separate,  already  existing  chamber  (not  connected  to  the  beam 
machine)  is  being  renovated  and  \will  be  used  to  conduct  tests  on  the  detec- 
tor, along  with  the  addition  of  a quadrapole  mass  filter  to  the  detector.  This 
chamber  has  windows  in  it,  and  is  already  equipped  with  a QMS.  The 
pump  down  time  is  much  faster  than  for  the  beam  machine,  and  it  is  thus 
hoped  the  beam  machine  can  be  used  at  the  same  time  that  modifications 
are  being  separately  tested  in  the  smaller  test  chamber. 

It  is  believed  that  the  current  detector  will  be  good  enough  to 
allow  scattering  data  to  be  taken  'or  some  pairs  of  gases.  It  will  be 
necessary,  however,  to  make  some  changes  to  the  detector  if  full  use  is 
to  be  made  of  the  beam  machine's  capabilities. 
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A.  MODELING  SUBSURFACE  RESISTIVITIES  WITH 

MAGNETOTELLURIC  DATA  (MT  DATA)  t 
Professor  H.  W.  Smith,  Professor  F,  X,  Bostick,  Jr., 
J , E.  Boehl,  S.  H,  Cheng,  and  D.  E.  Wight 


1.  MAGNETOTELLURIC  AND  DIPOLE-DIPOLE  SOUNDINGS  IN  NORTHERN 
WISCONSIN 

This  project,  which  has  been  concerned  with  the  determination 
of  the  deep  electrical  conductivity  structure  of  the  earth's  crust  using 
dipole-dipole  and  magnetotelluric  soundings,  has  been  completed  during 
the  year.  The  techniques  developed  are  considered  to  be  significant  to 
the  investigation  of  the  possibility  of  subsurface  communication  and  to 
the  site  selection  and  antenna  pattern  determination  for  U.S.  Navy  Pro- 
jects Sanguine  and  Seafarer. 

This  work  has  been  reported  at  meetings  and  conferences,  in 
one  Masters  Thesis,  as  part  of  the  1975-1976  Significant  Accomplishments 
of  the  Electronics  Research  Center,  and  is  currently  being  revised  for 
publication  in  the  Geophysical  Monograph  Series  of  the  American  Geo- 
physical Union  [1,  2,  3,  4l. 


2.  A LAND  TO  SEAFLOOR  ELECTROMAGNETIC  PROPAGATION  EXPERIMENT 
IN  THE  OLYMPIC  PENINSULA 

During  the  summer  of  1976,  scientists  from  Pacific  Sierra 
Research  Corporation,  the  above  group  from  The  University  of  Texas  at 
Austin,  and  Scripps  Institute  of  Oceanography  performed  an  electromag- 
netic propagation  experiment  on  the  Olympic  Peninsula  in  the  northwest 
corner  of  the  state  of  Washington.  The  main  purpose  of  the  experiment 
was  to  send  ultra-low-frequency  (0.1  to  10  Hz)  signals  from  land  to  deep 
receivers  (to  8,000  feet)  on  the  seafloor.  In  addition,  a number  of  sub- 


t This  research  was  supported  in  part  by  the  Joint  Services 
Electronics  Program  under  Contract  F44620-76-C-0089 , and  in  part  by 
the  Office  of  Naval  Research  Contract  N00014-67-A-0126-00'04;  Office 
of  Naval  Research  Contract  N00014-76-C-0620  sponsored  by  ARPA 
Order  No.  3189;  and  National  Science  Foundation  Grant  GA  38827. 
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sidiarv  land-based  measurements  and  magnetotelluric  (MT)  soundings 
were  carried  out.  The  operation  of  a grounded  horizontal  electric  dipole 
transmitter  as  well  as  all  land-based  measurements  was  made  by  The 
University  of  Texas  team.  All  sea-based  measurements  were  made  by  a 
team  headed  by  Professor  C.  S.  Cox  of  The  Scripps  Institute  of  Ocean- 
ography and  overall  management  of  the  project  was  the  responsibility  of 
Dr.  E.  C,  Field  and  Mr.  John  Farquhar  of  Pacific  Sierra  Research  Corpo- 
ration. 


Square-wave  signals,  using  frequencies  from  0.1  to  10  Hz 
and  a peak  moment  of  1.6  x 10^  amper-meters , were  injected  into  the 
earth  by  the  transmitter.  Strong  signals  with  a 20-30  dB  signal  to  noise 
ratio  for  the  fundamental  were  received  on  a 100  m long  electrode  pair 
emplaced  on  the  seafloor  at  a depth  of  about  1,000  feet  at  a range  of 
22  km  from  the  transmitter.  It  is  concluded  that  these  signals  were  from 
an  over-down  mode.  However,  measurements  on  the  seafloor  at  ranges 
from  110  to  13  5 km  at  a depth  of  8,000  feet  using  receiving  antennas 
540  m to  1,000  m in  length  failed  to  detect  any  signal.  It  is  concluded 
that  the  ocean  screened  out  the  over-down  mode  and  that  any  signal  de- 
tected v/ould  have  had  to  propagate  through  the  crust.  Calculations  from 
these  results  indicate  that  a field  would  have  been  detected  at  135  km 
and  8,000  feet  depth  if  the  effective  conductivity  of  the  propagation  path 
were  less  than  3 x 10  ^ mho/m  . A report  [4]  on  the  land-to-seafloor  por- 
tion of  the  project  gives  details  of  the  measurements  and  conclusions. 

The  land-based  measurements  consisted  of  Magnetotelluric 
soundings  at  four  sites  and  DC  dipole-dipole  resistivity  soundings  at 
fifteen  sites  ranging  from  about  4 km  to  the  seaward  side  (northwest) 
to  120  km  to  the  landward  sMe  (east  and  southeast)  of  the  transmitting 
antenna.  On  the  basis  of  tne  earth  model  indicated  by  the  resistivity 
results  from  the  two  sets  of  measurements,  it  would  appear  the  ULF  elec- 
tromagnetic signals  could  be  propagated  reasonably  well  in  the  lithosphere 
to  distances  in  excess  of  100  km  along  a path  extending  to  the  landward 
(east)  of  the  transmitting  antenna.  These  measurements  also  show  a sig- 
nificant increase  in  the  effective  conductivity  at  depth  at  sites  to  the 
seaward  side  of  the  transmitter.  This  suggests  that  either  the  effective 
conductivity  beneath  the  ocean  is  somewhat  greater  than  that  beneath 
the  land  sites,  or  that  there  is  a conductive  obstacle  in  the  crust  between 
the  transmitter  and  the  sea-based  receivers.  The  latter  possibility  could 
be  caused  by  the  presence  of  a subduction  zone,  albeit  inactive,  near 
the  present  shore  line.  It  is  also  a possible  explanation  for  the  negative 
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results  at  the  longer  ranges  for  the  seafloor  propagation  experiment.  These 
results  are  contained  in  a recent  report  [5]  describing  the  land-ba:^ed  por- 
tion of  the  experiment. 


3.  DEVELOPMENT  OF  MAGNETOTELLUR  10  ANALYSIS  AND  INTERPRETATION 
METHODS 

It  has  been  recognized  for  some  time  that  on-site,  real-time 
analysis  and  interpretation  of  Magnetotelluric  (MT)  and  Audio-Magneto- 
telluric(  AMT)  data  is  both  a highly  desirable  and  a reasonable  goal.  To 
be  practical,  however,  the  very  large  data  accumulation  and  processing 
reguirements  of  conventional  methods  must  be  simplified  wherever  possible. 
Substantial  progress  toward  this  goal  has  been  made  by  a combination  of 
hardware,  in  the  case  of  the  AMT  system,  and  by  the  development  of  simp- 
lified analysis  and  inversion  techniques.  Two  new  techniques  which  have 
been  used  on  an  experimental  basis  are  now  fully  operational  and  have  been 
reported  during  this  contract  period. 


The  first  is  a simple,  approximate  method  of  MT  analysis  and 
inversion  and  is  based  on  the  asymptotic  behavior  of  apparent  resistivity 
curves  at  low  frequency.  It  utilizes  the  fact  that  the  asymptote  for  a 
single  section  of  uniform  resistivity  over  a perfectly  conducting  substrate 
depends  only  on  the  depth  to  the  substrate,  and  is  independent  of  the  con- 
ductivity of  the  resistive  section  even  if  that  conductivity  is  a function  of 
depth.  Also,  the  asymptote  for  the  same  resistive  section  over  a perfectly 
insulating  substrate  depends  only  on  the  integral  of  the  conductivity  versus 
depth  profile  to  the  depth  of  the  substrate.  The  simple  method  has  been 
disclosed  to  a number  of  researchers  who  have  used  it  successfully  in  a 
wide  variety  of  situations.  It  has  been  presented  at  a conference  and 
has  been  published  in  the  Conference  Proceedings  [6]. 

The  second  technique  uses  the  Hilbert  Transform  to  incorporate 
phase  as  well  as  magnitude  information  in  the  final  determination  of  the 
apparent  resistivity  versus  frequency  curve,  and  has  greatly  reduced  the 
amount  of  experimental  scatter  which  can  often  be  excessive  at  noisy 
sites.  The  process  has  been  termed  "phase  smoothing"  and  leads  to  what 
is  essentially  an  unbiased  estimate.  The  process  has  been  presented  at 
a conference  and  is  to  be  included  in  an  Appendix  to  the  AGU  Monograph 
paper  cited  in  Section  1 above. 
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"Ocular  Damage  Due  to  Light  Generated  Retinal  Temperature  Rise," 

Professor  A.  J.  Welch,  Principal  Investigator,  April  1,  1973  through 
March  31,  1976. 

Department  of  Health,  Education  and  Welfare  PHS  Grant  HL  18247-01, 
"Induced  AV  Nodal  Conduction  Changes  in  Normal  Hearts,"  Professor 
P,  C.  Richardson,  M.D.  , Principal  Investigator,  June  1,  1975  through 

May  31,  1976. 
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(RESEARCH  GRANTS  AND  CONTRACTS  FROM  FEDERAL  SOURCES) 


U.S.  Coast  Guard,  DOT-CG-63441-A , "A  Study  of  Broad  Band  High  Level 
Audio  Systems,"  Professor  E,  L,  Hixson,  Principal  Investigator,  Septem- 
ber 28,  1976  through  September  27,  1977. 

U.S,  Department  of  Commerce,  National  Bureau  of  Standards,  "Second 
International  Conference  on  Software  Engineering,"  Professor  R.  T.  Yeh, 
Principal  Investigator,  September  30,  1976  through  March  31,  1977 


OTHER  THAN  FEDERAL  FUNDS 

Bureau  of  Engineering  Research,  "Studies  of  Concurrent  Systems  Using 
Petri  Nets  and  Related  Methods,"  Professor  Tilak  Agerwala , Principal 
Investigator,  September  1,  1975  through  August  31,  1976. 

Bureau  of  Engineering  Research,  "Correctness  and  Fault  Tolerance  in 
Concurrent  Computer  Systems,"  Professor  Tilak  Agerwala,  Principal  Inves- 
tigator, September  1,  1976  through  August  31,  1977. 

Bureau  of  Engineering  Research  Award,  "Multilayer  Dielectric  Devices 
for  Integrated  Optics ,"  Professor  A.  B.  Buckman,  Principal  Investigator, 
September  1,  1976  through  August  31,  1977, 

Bureau  of  Engineering  Research  Contract  (Tracer,  Inc.),  "Research  Relating 
to  the  Design  and  Development  of  a Feasibility  Model  of  a Bone  Conduction 
Calibrator,"  Professor  E.  L,  Hixson,  Principal  Investigator,  September  5, 

1975  through  March  6,  1976. 

Bureau  of  Engineering  Research  Initiation  Award,  "Estimation  and  Analysis 
of  Nonlinear  Stochastic  Systems , " Professor  S.  I.  Marcus,  Principal  In- 
vestigator, September  1,  1975  through  August  31,  1976. 

Bureau  of  Engineering  Research  Grant,  "Analysis  of  Linear  Systems  with 
Multiplicative  Poisson  Noise , " Professor  S.  I.  Marcus,  Principal  Inves- 
tigator, September  1,  1976  through  August  31,  1977. 

Bureau  of  Engineering  Research  Award,  "Research  on  Solid  State  Interface 
Reactions;  Surface  Profiling  of  Impurity  and  Defect  States  at  Semiconductor 
Interfaces,"  Professor  R.  M,  Walser,  Principal  Investigator,  September  1, 

1976  through  August  31,  1977, 


168 


(RESEARCH  GRANTS  AND  CONTRACTS  FROM  FEDERAL  SOURCES) 


Bureau  of  Engineering  Research,  "New  Faculty  Research  Initiation  Award," 
Professor  G.  L.  Wise,  Principal  Investigator,  September  1,  1975  through 
July  31,  1977. 

International  Research  and  Exchanges  Board  Award  to  Professor  B.  F. 
Womaclc  for  support  of  Visiting  Soviet  Scholar,  August  1,  1976  through 
May  31,  1977. 

Japanese  Ministry  of  Education  Grant,  "Experiments  in  Multimodal  Search 
Techniques,"  Professors  M.  Oda  and  B.  F.  Womaclc,  October  1,  197  5 
through  September  30,  1976. 

Research  Corp.  , #275-40-4062,  "Study  of  the  Chemical  Kinetics  Important 
to  the  Production  and  Loss  of  Excited  Molecules  in  Electron  Beam  Excited 
Gas  Mixtures,"  Professor  J.  W.  Keto,  Principal  Investigator,  November  1, 
1976  through  November  1,  1977. 

The  Robert  A.  Welch  Foundation  Grant  F-534,  "Electron  Scattering  from 
Alkali  Halides  Vapors,"  Professor  Manfred  Fink,  Principal  Investigator, 
May  1,  1973  through  May  31,  1979. 

The  Robert  A.  Welch  Foundation  Grant  F-600,  "Spectroscopy  of  van-der- 
Waals  Molecules  and  of  Biomolecular  Pairs,  " Professor  Lothar  Frommhold, 
Principal  Investigator,  June  1,  1975  through  May  31,  1978. 

Sandia  Laboratories  Contract  No.  02-7963,  "Magnetotelluric  Survey  at 
Kilauca  IKi  Lava  Lake,"  Professors  F.X.  Boxtick  and  H.  W Smith, 

Principal  Investigators,  May  5,  1976  through  September  15,  1976. 

Sandia  Laboratories  Grant  04-77731,  "Chemical  Kinetics  and  Ratiative 
Processes  Important  to  the  O^So  Laser,"  Professor  J . W.  Keto,  Principal 
Investigator,  February  6,  1976  through  February  7,  1978. 

Texas  Atomic  Energy  Research  Foundation  Contract  for  "Controlled  Thermo- 
nuclear Fusion  Research  and  Engineering,"  Professors  M.  F.  Becker, 

J.  Keto,  and  E.  J.  Powers,  Principal  Investigators. 

URI  Grant  R-268,  "The  Use  of  Harmonic  Analysis  in  Estimation  and  Control 
Theory,"  Professor  S.  I.  Marcus,  Principal  Investigator , September  1,  1975 
through  August  31,  1976. 
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